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Abstract 



We consider a number of aspects of the problem of defining time observables in quantum 
theory. Time observables are interesting quantities in quantum theory because they of- 
ten cannot be associated with self-adjoint operators. Their definition therefore touches on 
foundational issues in quantum theory. 

Various operational approached to defining time observables have been proposed in the 
past. Two of the most common are those based on pulsed measurements in the form of strings 
of projection operators and continuous measurements in the form of complex potentials. One 
of the major achievements of this thesis is to prove that these two operational approaches 
are equivalent. 

However operational approaches are somewhat unsatisfying by themselves. To provide a 
definition of time observables which is not linked to a particular measurement scheme we 
employ the decoherent, or consistent, histories approach to quantum theory. We focus on the 
arrival time, one particular example of a time observable, and we use the relationship between 
pulsed and continuous measurements to relate the decoherent histories approach to one based 
on complex potentials. This lets us compute the arrival time probability distribution in 
decoherent histories and we show that it agrees with semiclassical expectations in the right 
limit. We do this both for a free particle and for a particle coupled to an environment. 

Finally, we consider how the results discussed in this thesis relate to those derived by 
coupling a particle to a model clock. We show that for a general class of clock models the 
probabilities thus measured can be simply related to the ideal ones computed via decoherent 
histories. 
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1 Introduction 



The time is out of joint; O eurscd sprite, 
That ever I was born to set it right! 

Shakespeare, Hamlet 

1.1 A Brief History of Time in Quantum Theory 

The treatment of time observables is an important loose end in quantum theory [81, 80]. 
Whilst there is no indication of any disagreement between the predictions of quantum theory 
and the outcomes of experiment, time observables constitute a class of quantities that are 
observables in classical mechanics, but for which no satisfactory treatment exists in standard 
quantum theory. 

As well as being an important topic in its own right, a consistent treatment of time 
observables might also shed light on problems in quantum gravity and quantum cosmology. 
One of the most tricky problems faced by any attempt to quantise gravity is that of the 
problem of time [64] . In essence the symmetries of general relativity prevent the identification 
of any variable to play the role of time in quantum gravity. Whilst there is no generally agreed 
way of tackling this problem one school of thought holds that any approach to quantising 
gravity must proceed by treating time and space on an equal footing. A theory treating time 
and space in a truly symmetric way would have as observables probabilities for entering given 
regions of spacetime rather than, for example, statements about the position of a particle 
at a fixed moment of time. Since these are not the kind of probabilities one normally deals 
with in quantum theory it is instructive to see if the same exercise can be performed in 
non-relativistic quantum theory. That is, can one assign probabilities to questions of the 
form, "What is the probability of finding a particle in a given region of space. A, in a given 
interval of time, [to,^i]?" 

Asking questions of this sort quickly leads us to another problem in non-relativistic quan- 
tum theory, which is this: "What exactly is the status of the variable t that appears in 
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Schrodinger's equation?" It should be clear immediately that it is not related to any oper- 
ator on the Hilbert space of the system. Indeed it is not really a quantum quantity at all 
and is probably best thought of as referring to some external classical clock. This hybrid 
approach to dynamics, having a quantum wavefunction depend on an external classical clock 
variable, could be seen as troubling. An operationalist might reply that the co-ordinate x 
also refers to an external classical world, but this line of reasoning is over simplistic for three 
reasons. The first is that both position and its conjugate variable, momentum, appear as 
operators in quantum mechanics, whereas energy is represented by an operator in quantum 
mechanics, but time is not. There is therefore a lack of symmetry in the treatment. The 
second reason is that quantum mechanics may be written in a representation-free manner 
in terms of abstract vectors in Hilbert space. Here there need be no mention of position 
or momentum, and yet the Schrodinger equation for evolution in this abstract space still 
includes the time coordinate. The final reason is the most profound, and again touches on 
the lack of symmetry in the treatment of space and time coordinates. In non-relativistic 
quantum mechanics, the object. 



gives the probability of finding the particle in the interval [x, x + dx] at the time t. What 
is the corresponding probability for finding the particle at the position x in the interval 
[t,t + dt]'? It is not given by 



since finding the particle at a position x at the times ti or t2 are not exclusive alternatives, 
see Fig. (1.1). Yet in the standard interpretations of quantum mechanics the wavefunction is 
supposed to contain all the information about a given system [67]. If the square norm of a 
general wavefunction is not a probability distribution on a set of times does this mean there 
can be no satisfactory treatment of time observables in standard quantum theory? 

To get a better idea of where the difficulty lies it is instructive to look more closely at the 
classical case. Here the solution of a particular problem may be expressed in terms of an 
equation of motion for the particle, of the form. 



where Xq represents the initial condition. This is the solution to the question: "Where is 







X = f{t,Xo,Po) 



(1.3) 
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t t 




(a) Probabilities for Space-like Surfaces (b) Probabilities for Time-like Surfaces 

Figure 1.1: (a) Since every trajectory crosses a surface of constant time only once, passing 
through disjoint intervals A or B at a given time represent exclusive alternatives. 
(b) In contrast trajectories may cross a surface of constant position many times. 
Therefore passing through disjoint intervals A or B are not exclusive alternatives. 
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the particle at time i?" Under general conditions one may invert this equation of motion to 
give, 



which is the solution to the question: "When is the particle at position x?" Crucial to this 
analysis is the concept of a trajectory and that the particle will with certainty be found 
somewhere on this trajectory, at any time. This concept is of course lacking in quantum 
theory. 

In classical mechanics observables are given by functions on phase space, and so Eq.(1.4) 
is equivalent to some phase space function. 



Crucially, unlike classical observables like the energy or angular momentum, construction 
of Eq.(1.5) requires knowledge of the equations of motion, Eq.(1.3). Strictly speaking then, 
even in classical mechanics questions such as "When does the particle reach the origin?" have 
a different character to questions like "What is the energy of the particle?" . The observable 
coresponding to the former is a function on phase space which is constructed using the 
equations of motion. In a sense, therefore, the surprise is not that quantum theory struggles 
to deal with time observables, but that classical mechanics deals with them so easily. This 
is essentially a consequence of the realist nature of classical physics. 

Approaches to the problem of defining time observables in quantum theory fall roughly into 
three camps. Although it is not the intention of this thesis to provide a complete overview 
of the literature on time observables (see Ref. [77] for a detailed review in the case of arrival 
time), the following characterisations will be useful in understanding how the approaches 
taken in this work relate to previous studies. 

Trajectory Before Quantisation 

These are approaches which in some way borrow the classical structure of a trajectory and 
use it in the quantum analysis. One way to do this is to begin with the classical description, 
in terms of Eq.(1.3) and use this to construct the classical phase space observable Eq.(1.5), 
as outlined above. One then tries to quantise this quantity to obtain the operator. 



t = F{x,Xo,po) 



(1.4) 



t = t{x,p). 



(1.5) 



(1.6) 
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This category of approaches includes the Aharanov-Bohm time operator [3] . Apart from the 
obvious issues surrounding whether Eq.(1.6) is well defined for a given classical observable, 
there is something troubling about this analysis. Since the quantum system we are trying to 
analyse will not be explicable in terms a particle following a single trajectory, it is not clear 
whether using the classical trajectory to obtain the operator Eq.(1.6) is a valid procedure. 
Even if the resulting operator is well defined the relation between this observable and the 
results of measurement remains unclear. Since this is essentially the question of how to un- 
derstand the quantisation of a given classical system we are unlikely to make much headway 
with this issue here. It does, however, motivate trying to find an alternative way to proceed. 

Other approaches that fall into this category include the analysis of Kijowski and its vari- 
ants, see eg Refs.[68, 61]. Here one begins by deciding what properties one would expect 
the quantum probability distribution to possess, such as time translation covariance or pos- 
itivity, based on the properties of the classical observable. One then uses this "wish list" 
of properties to determine the correct operator, with the eventual hope that a sufficiently 
detailed list of properties might uniquely determine this operator. It is not clear that this 
approach should work in all cases however: it is easy to think of quantum observables, such 
as angular momentum or the energy levels of a harmonic oscillator, that have a very different 
character from their classical counterparts. 

The general attitude of this thesis is that whilst any of the analyses above might arrive at 
the correct definition of a time observable in a particular case, they cannot be relied upon 
to do so in all cases. 

Trajectory After Quantisation 

These approaches include those based on the de Broglie-Bohm (dBB) interpretation of quan- 
tum theory [63], analyses based on the WKB approximation [40] and also the decoherent 
histories approach presented in this thesis. Here one first quantizes the system and then 
looks for structures to play the role of the classical trajectories used to define the classical 
observables. For example for gaussian states the center of the wavepacket will follow simple 
curve in configuration space. Provided one does not perform measurements that disturb this 
state this limited notion of trajectory can be used to find approximate arrival and dwell times 
for the particle, limited by the width of the wavepacket. Of course these states form only a 
hmited class of possible states, and it is a general feature of this category of approaches that 
they may fail for highly non-classical states. 
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In the DH approach, described in more detail below, the aim is to assign probabilities to 
coarse grained histories of a quantum system. Provided these probabilities can be assigned 
consistently one can then use these histories in much the same way as the classical trajectories 
to discus time obscrvablcs. 

In the dBB interpretation, the wavefunction is supplemented by a real particle follow- 
ing a trajectory computed from the wavefunction. Since a quantum particle always comes 
equipped with a unique trajectory computation of time observables is, at least conceptually, 
as trivial as in the classical case. The price one pays, however, is that the behavior of these 
trajectories generally depends on whether or not one choses to carry out a measurement on 
the system. Thus although in dBB probabilities for time observables can be obtained with- 
out the need for an observer, the probabilities obtained in this way do depend on whether 
or not there is an observer. We will not have much to say in this thesis about dBB, for more 
information see Ref.[63]. 

Trajectory Free Approaches 

These approaches may also be called "operational." Here one takes a completely different 
tack and instead looks for other ways to define time observables classically, and then tries 
to repeat this analysis quantum mechanically. 

These approaches include studies based on model clocks, see Ref . [88] and Chapter 8 of this 
thesis. Here the idea is to couple the system we wish to measure to some additional degrees 
of freedom which function as a clock. In this context a clock is simply a quantum system 
coupled to our system of interest in such a way that the final state of the clock is correlated 
with the time observable we wish to measure. More details about this scheme can be found 
in Chapter 8. 

In classical mechanics many time observables can be defined by imposing absorbing bound- 
ary conditions on the system. See Chapter 5 for a full discussion in the arrival time case. 
The analogue of these absorbing boundary conditions in quantum theory in the inclusion of 
a complex absorbing potential, see Ref. [76] and also Chapter 5 of this thesis. 

Typically the probability distributions obtained from these approaches are complicated 
functions of the properties of both the particle and measuring apparatus. There may or may 
not exist limits in which dependence on the measurement parameters drops out, depending 
on the set up. 
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As well as the technical differences between these approaches the different approaches 
display very different attitudes towards the observables we are trying to define. In particular 
the first and third approaches assume that an observable such as the arrival time of a particle 
is always a well defined quantity and that for any system one may obtain the probability 
distribution of that observable. The second class of approaches, in contrast, are based on 
recovering the classical notion of trajectory after quantisation, and thus they are generally 
only well defined for semiclassical systems. Proponents of these approaches would claim 
that quantities such as the arrival time distribution are only defined for suitably classical 
systems^. 

An obvious question is whether there is anything linking these three different approaches 
to time observables. This is an important question, firstly since although any particular 
approach might claim to have the solution to these problems the validity of these approaches 
themselves is questionable. If we could demonstrate that for a particular time observable 
the same results are obtained via a number of approaches this would give us confidence that 
this result is correct. Secondly, the different approaches discussed above have their roots in 
different interpretations of quantum theory. Although we will try and avoid philosophical 
questions about the foundations of quantum theory in this thesis if we were to show that, 
for example, expressions for time observables in de Broglie-Bohm theory were identical to 
those obtained via decoherent histories this might shed interesting light on the relationship 
between different interpretations of quantum theory ^. 

1.2 The Quantum Zeno Effect 

There are a number of interesting quantum effects relevant to the definition of time observ- 
ables, the most important of which is probably the quantum Zeno effect. The effect was first 
noticed by AUcock in his seminar scries of papers on the arrival time problem [6], but the 
first serious study was undertaken by Misra and Sudershan in Ref . [74] . 

The quantum Zeno effect is often expressed as a "freezing" of evolution due to repeated 

^In much the same way as it is claimed that time itself is only emergent, from an ultimately timeless theory 
of quantum gravity [13]. 

^Not to mention the fact that although any interpretation of quantum theory is constrained by the require- 
ment that it must agree with the predictions of standard Copenhagen quantum theory, since Copenhagen 
quantum theory struggles to incorporate many time observables there is a possibility that different inter- 
pretations of quantum theory may be free to make different predictions for time observables. There may 
therefore be experimental tests that could be done to distinguish between them. Exciting as this line of 
thought is, we shall not dwell on it in this thesis. 
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measurement. Much has been learned about the phenomenon since the early investigations 
however and it is now understood that there is more to the effect than simply preventing 
evolution. Repeated projections onto a subspace allow for unitary evolution within that 
subspace, the effect of the repeated measurement being simply equivalent to imposing "hard- 
wall" boundary conditions, preventing the particle from leaving the subspace [27]. It is also 
known that there also exists an "anti-Zeno" effect, that is, a speeding up of evolution due 
to repeated measurement [26], although we shall have nothing to say about the anti-Zeno 
effect in this thesis. 

The quantum Zeno effect plays an important role in study of time observables because, 
at least for the "trajectory after quantisation" and "trajectory free" approaches outlined in 
the previous section, the definition of these observables involves statements about the state 
of the system at multiple times. This means one has to consider sequential measurements 
on the system of interest, if these are performed too frequently the quantum Zeno effect 
will result. However the situation is complicated by the fact that most discussions of the 
quantum Zeno effect take place in a setting appropriate for quantum optics or quantum 
information, ie a system with a small number of discrete energy levels, not in the context 
of particles in ID. In the rest of this section we therefore outline the standard discussion of 
the quantum Zeno effect and the way in which this discussion might be modified to make it 
relevant to time observables. Our discussion will be brief, since the details of the Zeno limit 
in the arrival time problem etc are worked out more thoroughly in Chapter 4. 

We begin by introducing the Zeno effect in the way it appears in the original discussion of 
Misra and Sudarshan, Ref. [74]. The idea is to utilise some form of product formula to show 
that the repeated process of evolution followed by measurement is equivalent to restricted 
evolution in the subspace defined by the measurement. We consider a system with initial 
state which evolves under a Hamiltonian H for a total time T, interspersed with N 
projective measurements P. The state at time T is given by, 

= [ex^{-iHT/N)Pf (1.7) 

One then considers the limit N ^ oo, and one can show, under certain conditions on H and 
P, 

V^{T) = lim Vn{T) = exp{-iHzT)P, (1.8) 
where the Zeno hamiltonian Hz — PHP. This shows that evolution under repeated mea- 
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surement leads to unitary evolution in the subspace defined by the measurements. 

If we consider the evolution in Eq.(1.7), with P = 6{x) and H the free hamiltonian for 
a particle in ID, we have a simple model of the type of measurement set up we might use 
to determine the time at which the particle arrived at a; = 0. Without concerning ourselves 
too much with the details, since these are covered in Section 1.4 and the rest of this thesis, 
it is clear that the Zeno effect sets a lower limit on the resolution with which we can define 
an arrival time using this set up. 

Although this analysis is mathematically precise this precision is obtained at the expense 
of providing any great physical insight. A key question concerns the time scale on which the 
Zeno effect begins to occur. That is, what is the minimum spacing between measurements, 
T/N, that I can take without causing significant reflection from the measuring aparatus. 
This is important since we can hardly expect to achieve iV — )■ oo in the real word. One way 
to approach this is to look for the time scale on which the state may be said to leave the 
subspace defined by P. The idea is that projections more frequent than this will "interrupt" 
the state in the process of exiting the region, so that this time scale provides an upper bound 
on the separation between projections necessary to give rise to the Zeno effect. We follow 
the analysis of Peres [87] and consider an initial state, with evolution for a time t followed by 
a measurement which projects back onto the original state. The survival probability, S{t), 
is given by, 

S{t) = \\{ilj\exp{-iHt)\ilj)\\'^ 

= l-t\AHf + 0{t^)^l-{t/T,y + 0{t^) (1.9) 

where we have defined the "Zeno time"^, 

= {AHr = |^) - (^| H . (1.10) 

We see from this expression that the short time behaviour is quadratic in t, thus we confirm 
that in the limit of infinitely frequent measurements the state does not evolve. Secondly, 
we see that there is a characteristic time scale, tz — {AH)~^, on which the evolution takes 
place. 

Although we have arrived here at what is normally called the "Zeno" time, it is not clear 
that this really is the time scale on which the projections give rise to reflection and thus the 

^Note that there are many different notions of the time it takes a system to make a transition. See [92] for 
a more detailed discussion. 
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Zeno effect. For a minimum uncertainty gaussian wavepacket with spatial width a peaked 
around some momentum po Eq.(l.lO) gives = mcr/ |po|- This is indeed the time scale on 
which the state crosses the origin, however this is an essentially classical timescale and it 
seems unlikely that this is relevant for reflection, which is a quantum effect. It is also worth 
pointing out that if we really are limited in the accuracy of our description by the time 
scale 1 / AH then it is impossible to formulate a description of crossing time probabilities in 
the normal sense. This is because our minimum temporal resolution is the same order as 
the time taken for our wavepacket to cross the origin, so the probabilities for crossing will 
essentially be 1 for one time interval and for the rest. 

Projecting back onto the original state captures the notion of the time on which the state 
changes, and this is certainly the same as leaving some subspace of the full Hilbert space, but 
we are interested in a particular subspace, that of states with support only in ,x > and this 
analysis docs not capture that. There is another problem, which is that Eq.(1.9) is a Taylor 
expansion of the initial state about t = 0. It is known that for general wavefunctions the 
behaviour may not be analytic at t = 0. Put another way, the Zeno time is expressed in terms 
of moments of the wavefunction, and these may not exist for general wavefunctions. This is 
particularly apparent when one tries to extend this analysis to more general measurements. 

One of the major achievements in this thesis is to show, in Chapter 4, that the minimum 
temporal resolution is in fact set by 1/ where E is the energy of the system, and not by 
Eq.(l.lO). 

1.3 The Backflow Effect 

Another quantum effect relevant to the understanding of time observables is the backflow 
effect [12, 25]. Consider an initial state \ip) consisting entirely of negative momenta evolved 
freely for a time t and compute the standard Schrodinger probability current at the origin, 
defined in such a way that classically it should be positive. This can be written in operator 
form as, 

J{t) = ^ (VI e^'^'mx) + mP}e-'"' IV') ■ (1-11) 

(We will discuss the current in much more detail in later sections.) Now the operator 
{—){p5{x) + 5{x)p) is not positive, even when restricted to act on states with negative mo- 
mentum. This means although Eq.(l.ll) should be positive classically, it need not be quan- 
tum mechanically. What is remarkable is that very little is known about this effect. One 
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thing that is known is that the amount of probabihty that can flow backwards is bounded, 
in the sense that, 

r dtJ{t)>\ (1.12) 

where A is some dimensionless constant computed numerically to be A ~ —0.0384517 [25]. 
The existence of such a bound is somewhat unexpected, Bracken and Melloy in Ref. [12] con- 
jecture that it is a "new dimensionless quantum number." Certainly there is no dependence 
on % in this constant and as such it has no obvious classical limit. 

One interesting set of questions concern the state for which this bound is obtained, does 
it have an analytical expression? Does it have an obvious physical interpretation? Another 
question concerns the typicality, or otherwise, of backflow. Do relatively standard wavefunc- 
tions such as gaussian wavepackets or superpositions of these give backflow, or do we require 
more pathalogical states? If we can indeed produce backflow with more familiar states, is it 
possible to get close to the maximum amount of backflow? 

Some work towards addressing these questions is in progress [102], but much remains 
unclear. A better understanding of the backflow effect would shed much light on the problem 
of deflning time observables in quantum theory. This is because, together with the Zeno 
effect, it is these quantum effects that provide the fundamental limitation on the accuracy 
with which time observables can be defined. 



1.4 The Arrival Time Problem in Quantum Mechanics 
1.4.1 General Theory 

Much of this thesis will be concerned, directly or otherwise, with the question of analysing 
the arrival time problem in quantum theory, pictured in Fig. (1.2). It therefore seems useful 
to give a brief account of the problem here and also discuss some of the features any solution 
must posses. 

The classical "arrival time problem" is the following, "Given a free particle which, at i = 
has position and momentum pq, what is the time at which the particle crosses the origin?" 
The solution, of course, is that the particle crosses at a time. 



In order to compare the definition of time observables in classical and quantum physics it is 
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t=0 



Figure 1.2: The Arrival Time Problem: What is the probability that an incoming wavepacket 
crosses the origin during the time interval [ti,t2\ ? 

useful to work in terms of classical phase space distributions. For a free classical particle in 
ID the phase space probability distribution Wt{q,p) obeys the following, 



dw p dw 
dt m dq 

The solutions of this equation are of the form, 



(1.14) 



Wr{q,p) =Wo{q-pT/m,p), (1.15) 

so that evolution is just a linear transformation in phase space. 

We can write the probability of crossing at a time r in terms of this distribution as, 

dpdq5{T + mq/p)wo{q,p) = dpdq Wr{q,p) (1-16) 

For states consisting of entirely negative momenta this may also be written as, 

nc(r) = y dpdq0{q)wriq,p), (1.17) 

which can be seen by using Eq.(1.14) in Eq.(1.17) and then integrating by parts. 
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Eqs.(1.16), (1.17) express the arrival time distribution in two distinct ways. Eq.(1.16) tells 
us the arrival time distribution is the expectation value of a certain phase space distribution, 
whilst Eq.(1.17) tells us the arrival time distribution can be defined in terms of the "survival 
probability." More precisely the arrival time probability is just the rate at which probability 
leaves the region x > 0, ie the fiux across a; = 0. We will see below that these two expressions 
suggest different ways of obtaining the quantum arrival time distribution. 

The quantum arrival time problem is essentially the following question, "For a given 
wavefunction what is the probability distribution on the set of times at which the system 
may be said to arrive at a; = 0?" For most situations of interest we can put the following 
extra constraints on the wavefunction; at some initial time t = the wavefunction consists 
only of negative momenta and the particle is concentrated in a; > 0. 

The challenge is to understand the correct quantum analogue, n(r), of Eqs.(1.16), (1.17), 
if indeed it exists. We need to do this subject to the following constraints, 

• n(r) must be a valid probability distribution, ie (irn(r) = 1 and n(r) > 0. 

• We must recover Eqs.(1.16), (1.17) in some classical limit. 

There are other constraints we could impose, indeed in Ref.[68] is was shown that de- 
manding that n(r) satisfy some set of conditions abstracted from the classical case is, al- 
most, enough to fix it identically. This is an example of a "trajectory before quantisation" 
approach, see Section 1.1. 

In some ways the most natural way to quantise nc(r) is via Eq.(1.16). This suggests the 
(formal) quantum expression. 



is some quantisation of Eq.(1.13). We hit a problem, however, which is that it has proven 
surprisingly difficult to find a way of defining f as a self-adjoint operator. This means 
Eq.(1.18) is merely formal and cannot be used to define n(T). Work has been done to try 
and get round this problem [77, 3], but we will not consider this further in this thesis. At 
any rate, the criticisms inherent to any approach of this type, as laid out in Section 1.1 apply 
here. 




(1.18) 



where the arrival time operator. 




(1.19) 
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1.4.2 Standard Forms of the Arrival Time Distribution 



We review here some of the standard, mainly semiclassical, formulae for arrival time. It is 
useful to collect these expressions in one place as we will refer to them repeatedly throughout 
this thesis. None of these semiclassical expressions are entirely satisfactory, but any more 
fundamental approach to defining arrival times must reproduce these expressions for some 
suitable class of states and time intervals. 

We consider a free particle described by an initial wavepacket with entirely negative mo- 
menta concentrated in x > 0. A widely discussed candidate for the distribution p{ti, ^2) is the 
natural quantum version of Eq.(1.17), the integrated probability current density [77, 78, 49], 

pitiM)^ I dtJ{t). (1.20) 
Defining P{t) — 9{xt), and introducing the Wigner function [10], defined for a state p{x, y) 

by 

Wt{p, Q) = ^J die-'^^p{q + | 9 - |), (1-21) 
J{t) can be written in any of the following forms: 



dt 

hi) 
2m 



{'4jt \ {pS{x) + S{x)p ) \tpt) 



dpdq^-^^Wtip,q) (1.22) 
m 

This probability is normalised to 1 when ti — and t2 00, and has the correct semiclassical 
limit [49]. This expression arises from considering the survival probability, that is, the 
probability that the particle is still in a; > at some time t, 

p{x>Q,t) = {ilJt\9{x)\ilJt). (1.23) 

The rate at which probability leaves a; > is then given by Eq.(1.22) and this is a candidate 
for the arrival time probability distribution. There is a problem with Eq.(1.22) however, 
which is that it need not be positive, even for states consisting entirely of negative momenta. 
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This is the backflow effect, discussed in Section 1.4. This means we cannot regard Eq.(1.20) 
as the correct arrival time distribution. 

The backflow effect arises because of interference between portions of the state that have 
crossed the origin, and that part which has yet to cross. In a limited sense then, this problem 
is similar to that of defining an arrival time distribution for a classical particle in a potential, 
where the particle may cross the origin and then re-cross at some later time. In the classical 
problem one imposes absorbing boundary conditions to remove that part of the state that 
has crossed the origin. The analogue in the quantum case is the use of complex absorbing 
potentials. This is an example of a "trajectory free" approach, as outlined in Section 1.1. 
Under suitable conditions this approach does indeed yield the the arrival time distribution 
as Eq.(1.20). See Chapter 5 for a detailed discussion of this. 

Note however that wavefunctions consisting of a single gaussian wavepacket never display 
backflow. This means that Eq.(1.20) may indeed be valid for these states. These wavefunc- 
tions are also WKB states and the WKB interpretation, an example of a "trajectory after 
quantisation" approach as laid out in Section 1.1, does indeed give Eq.(1.20) as the correct 
arrival time distribution, at least for these states. 

We mention here briefly that the arrival time analysis of Kijowski [68] gives, 

Ti{t) = -mp\"'mW\^t) (1-24) 

m 

which can be obtained by an operator rc-ordcring of J(t), but which has the virtue of being 
positive. We will not have anything to say about this distribution in this thesis, since it 
seems to be unrelated either to the histories analysis. Chapter 5, or to the behavior of ideal 
clocks. Chapter 8. Other authors have, however, shown that the distribution Eq.(1.24) does 
arise naturally in analyses of the time of arrival based on POVMs, see Ref. [24] for details. 

For arrival time probabilities defined by measurements, considered later in this thesis, one 
might expect a very different result in the regime of strong measurements, since most of the 
incoming wavepacket will be reflected at x = 0. This is the essentially the Zeno effect [74]. 
It was found in a complex potential model that the arrival time distribution in this regime 
is the kinetic energy density 

\l{t) ^ C {iPt\pS{x)P\iJt) . (1.25) 

where C is a constant which depends strongly on the underlying measurement model [23, 41]. 
(See Ref. [82] for a discussion of kinetic energy density.) However because the majority of 
the incoming wavepacket is reflected, it is natural to normalise this distribution by dividing 
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by the probability that the particle is ever detected, that is, 



U{t) 




{i't\PS{x)p\ipt) 



(1.26) 



where (p) is the average momentum of the initial state. This normahsed probabihty distri- 
bution does not depend on the details of the detector. This suggests that the form Eq.(1.26) 
may be the generic result in this regime, although a general argument for this is yet to be 
found. 

It is found in practice that measurement models for arrival times lead to distributions 
depending on both the initial state of the particle and the details of the clock or measuring 
device, typically of the form 



where 11 (t) is one of the ideal distributions discussed above and the response function R{t, s) 
is some function of the clock variables. (In some cases this expression will be a convolution). 
However, it is of interest to coarse grain by considering probabilities p{ti,t2) for arrival times 
lying in some interval [ti,t2]- The resolution function R will have some resolution time scale 
associated with it, and if the interval ^2 — is much larger than this time scale, we expect 
the dependence on R to drop out, so that 



This is the sense in which many different models are in agreement with semi-classical formulae 
at coarse-grained scales. 

In addition to a lack of positivity, there is a more fundamental problem with Eq.(1.20), 
which is that probabilities in quantum theory should be expressible in the form [84] 



Here Pa is a projection operator, or more generally a POVM, associated with the outcome 
a. Eq.(1.20) cannot be expressed in this form and we therefore conclude that it is not a 
fundamental expression in quantum theory. 




(1.27) 




(1.28) 



p(a) =TY(P„p). 



(1.29) 
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The conclusion then is that whilst there is no axiomatic approach that yields Eq.(1.20) as 
the correct arrival time distribution, there is extremely good evidence for it from semiclassical 
and other approaches. One aim of this thesis is to show how Eq.(1.20) may be derived from 
an underlying axiomatic approach to quantum theory, thus providing the justification for 
these semiclassical and operational approaches. 

1.5 The Decoherent Histories Approach to Quantum 
Theory 

The following is a brief review of the aspects of the decoherent histories approach to quantum 
theory relevant for this thesis. More extensive discussions can be found Refs.[30, 31, 37, 85, 
41, 22]. 

Alternatives at fixed moments of time in quantum theory are represented by a set of 
projection operators {Pa}, satisfying the conditions 

J]P, = 1, (1.30) 

a 

PaPb = SabPa, (1-31) 

where wc take a to run over some finite range. In the decoherent histories approach to quan- 
tum theory histories are represented by class operators which are time-ordered strings of 
projections, 

Ca^PaAtn)---PaAtl): (1-32) 

or sums of such strings^ [65]. Here the projections are in the Heisenberg picture and a 
denotes the string (ai, • • • a„). All class operators satisfy the condition 

J^Ca^l. (1-33) 

a 

Note that in the same way that Eqs.(1.30), (1.31) do not define a unique decomposition 
of the state of a system at a moment of time, Eqs.(1.32), (1.33) do not define a unique 



The proper framework in which expressions such as Eq.(1.32) are to be understood is the temporal logic 
framework of Isham and coworkers [65]. Note also that class operators may be defined involving an 
infinite number of times, provided appropriate care is taken over the definition of the infinite products 
involved [66]. 
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decomposition of the possible set of histories of the system. It is important to decide upon 
the set {Ca} to be used at the beginning of the analysis, and to avoid working with class 
operators from more than one set^ [38]. 

Probabilities are assigned to histories via the formula 

p{a) = Tr {C^pCi) . (1.35) 

These probabilities are clearly real and positive, however probabilities assigned in this way 
do not necessarily obey the probability sum rules, because of quantum interference. We 
therefore introduce the decoherence functional 

D{a,P)^TT(c^pCiy (1.36) 

which may be thought of as a measure of interference between pairs of histories. We require 
that sets of histories satisfy the condition of decoherence, which is 

D{a, 13)^0, a^p. (1.37) 

Decoherence implies the weaker condition of consistency, which is that 

ReD{a,a) = 0, (1.38) 

and this is equivalent to the requirement that the above probabilities satisfy the probability 
sum rules. In most situations of physical interest both the real and imaginary parts of D{a, (3) 
vanish for a ^ (3, and so we shall require Eq.(1.37) throughout this thesis. This condition 



^Strictly speaking, one cannot generate logical contradictions by considering class operators from more 
than one set but doing so can nevertheless lead to confusion. This is similar to the choice of a family of 
projectors at a single time, Eq.(1.31). Consider a single spin-i particle. Suppose one finds that, 

prob(s^ =t) = i prob(s^ =i) = i ,^ 
prob{s^ =t) = i, prob(s^ =i) = I 

where prob(s2 =t) is the probability that the spin will be measured to be "up" in the ^-direction. One 
cannot conclude from this that prob(s^ =t a-nd Sx =t) = 3: for example. 

In the same way, consider two different decoherent sets of histories represented by class operators 
{Ca,C/s...} and {Cq.', C/3'...}. Now suppose that using the first set we can conclude that p{a) = 1 
and using the second that p{a') = 1. There is no difficulty here. The two sets of histories provide 
complimentary, not contradictory descriptions of the system. Note also that Eq.(1.35) implies that the 
probability of a given history is independent of the set of histories in which it appears, and so we could 
never derive the genuinely contradictory result that p{a) is 1 in one set of histories and in another. 
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is related to the existence of records [31, 43]. Decoherence is only approximate in general 
which raises the question of how to measure approximate decoherence. The decoherence 
functional satisfies the inequality [22] 



D(a,/3)|2<p(a)p(/3). 



(1.39) 



This suggests that a sensible measure of approximate decoherence is 



(1.40) 



We note briefly that when there is decoherence Eq.(1.33) and Eq.(1.37) together imply 
that the probabilities p{a) are given by the simpler expressions 



Decoherence ensures that q{a) is real and positive, even though it is not in general. 

A few comments are in order about the place of decohcrcnt histories in the foundations 
of quantum theory. Although some early workers tried to present dccoherent histories as an 
"interpretation" of quantum theory, see in particular [37, 85], the more modern view is rather 
that the histories approach is a valuable tool in understanding how to apply quantum theory 
to closed systems, but not an "interpretation" in the usual sense. The histories approach 
does not try to explain the contextual, non-local nature of quantum theory in the way that, 
for example, the de Broglie-Bohm approach does, nor does it try to introduce new physics 
to solve the problems inherent in the quantum description of measurement, in the way that 
collapse theories do. Rather the histories approach is one in which there is no contridiction 
between the essential quantum nature of reality and the existence of a classical realm of our 
experience. The histories approach sees classical physics, together with classical logic, as 
emergent from the true quantum description of our world. One can of course carry out any 
particular measurement at any time and thereby obtain information about the properties of a 
system at one or more times, however the histories approach gives the conditions under which 
one may reason classically with these properties. In particular, if a set of histories satisfies 
the condition of decoherence, then one may treat the system as if^ it possesses the properties 

^Whist again we stress that it is beyond the scope of this thesis to enter into a detailed discussion of the 
philosophical underpinnings of quantum mechanics in general and the decoherent histories approach in 
particular, the choice of words is important here. Claiming that it is consistent to act as if a system 
possesses one of a number of properties is not the same as claiming that it really does possess them. 



g(a) = Tr(C„p). 



(1.41) 
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associated with these histories, regardless of whether one carries out a measurement to check 
this or not. 

For more details of the way in which decoherent histories may be thought of as an "inter- 
pretation" of quantum theory see [38, 39]. 



1.6 The Decoherent Histories Approach to the Arrival 
Time Problem: Introduction and History 



Now that we have introduced the general framework of decoherent histories, we turn to the 
question of applying it to the specific case of time observables. The aim is to derive the 
correct class operators for questions such as "Did the particle cross the origin in the interval 
[ti, ^2]?" Since the solution to this problem is the subject of Chapter 5 of this thesis, we will 
not present the correct class operators here, but rather we will discuss in general terms how 
they may be defined. This is partly to motivate the use of the histories approach, since it 
will quickly become apparent that time observables are easily framed in this language. We 
also wish, however, to discuss the history of the decoherent histories approach to the arrival 
time problem and in particular some of the misleading results that have been obtained in the 
past by considering class operators defined via path integrals. We will see that the reason 
for these previous misunderstandings was a lack of appreciation of the role of the Zeno effect 
in the definition of class operators for time observables and it is a major goal of this thesis 
to correct these issues. 

Let us begin with the general issue of defining class operators for entering some region 
of spacetime A. Some early papers on the decoherent histories approach gave the following 
analysis of the problem [59, 72, 54, 97, 98]. Consider the propagator between two general 
spacetime points, considered as a path integral. Partition the set of paths summed over 
according to whether they ever enter the region A or not. 



that do not enter A. Here is the restricted propagator, defined as a sum over all paths 
which never enter the region A. The path integral over paths that do enter A would appear 




Dxe 



iS 



(1.42) 



where denotes the sum over all paths that enter A, and J^^ denotes the sum over all paths 
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to be the class operator we are looking for. We may therefore define the class operators for 
entering and not entering A as, 

Cncituto) = 9r{tl,to) (1.43) 

Ccituto) = e-^(^i-*°) -5.(^1,^0). (1-44) 

These definitions were used by earlier worker on decoherent histories, in particular in the 
series of papers by Yamada and Takagi [98]. 

However there is a problem with these definitions, which is that they suffer from the Zeno 
effect. To see why this is the case consider the specific example of the arrival time problem. 
Let P = 6{x) denote the projection onto the positive x-axis. The restricted propagator in 
this case may be written as, 

gr{T, 0) = exp {-iPHPr) P (1.45) 

so is unitary in the Hilbert space of states with support only in x > [95, 27]. Compare 
Eq.(1.45) with Eq.(1.8). 

Differently put, an incoming wave packet evolving according to the restricted propagator 
undergoes total reflection, so never crosses x — Q. This means the probability of not crossing 
the origin is 1, whatever the initial state. In the work of Yamada and Takagi [98] and also 
in later work by Wallden [95] this manifests itself as the inability to assign probabilities to 
anything other than states symmetric about the origin, for which the arrival probabilities 
vanish. Clearly these results are unphysical. 

A moment's reflection reveals that this result is not that suprising. We are demanding 
that the class operator for not crossing is a sum over paths that spend exactly zero time in 
X < 0. Our intuition would suggest that this restriction is too sharp and that we should 
instead sum over all paths spending a time less than some time 5t m. x <Q. 

These heuristic notions are very difficult to make precise in the path integral formulation 
of decoherent histories however. Instead, it is more helpful to attempt to deflne C„c in the 
manner of Eq.(1.32). A natural starting point is 

= Pe-'^'P ■ ■ ■ e-'"'P, (1.46) 

where there are N projectors and r = A^e. We recover the restricted propagator if we take 
the limit e ^ 0, compare Eq.(1.46) with Eq.(1.7). In light of this, one idea is to simply 
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decline to take the limit e ^ and define the class operator for not crossing to be Eq.(1.46) 
for some e. Clearly if e is large enough the system will be monitored sufficiently infrequently 
to let the wave packet cross x = without too much reflection. Eq.(1.46) is an important 
expression for the rest of this thesis, however as it stands this proposal is somewhat vague. 
It is a major goal of this thesis to turn this heuristic proposal into a concrete way of defining 
class operators. The key question concerns the timescale e. These issues will form the subject 
of Chapter 4. 

A second option is to take the limit of e — t- in Eq.(1.46), but "soften" the projections 
to POVMs, that is, to replace P = 9{x) with a function which is approximately 1 for large 
positive X, approximately for large negative x, and with a smooth transition in between. 
Although we have introduced the decoherent histories approach to quantum theory in Section 
1.5 above in terms of projection operators, the discussion is readily generalised to include 
POVMs ^. It would be interesting to pursue this idea, but we shall not do so in this thesis^. 

1.7 Summary and Overview of the Rest of This Thesis 

The major goal of this thesis is to give a decoherent histories analysis of the arrival time 
problem, that will point to the way to a description of general time observablcs in quantum 
theory. We wish firstly to know what the correct quantum analogues of Eqs.(1.16), (1.17) 
are, and how they may be obtained from first principles. We also wish to understand how 
the scmiclassical result, Eq.(1.20), emerges in some limit from this quantum result. 

The key question concerns the definition of the appropriate class operators for these ob- 
servablcs. In Section 1.5 we saw that the Zeno effect limits the usefulness of the path integral 
definitions of the class operators, so our challenge is to work with class operators defined in 
the manner of Eq.(1.46). These class operators represent evolution in the presence of pulsed 
measurements. One of the major achievements of this thesis will come in Chapter 4, where 
we will show that such evolution is, under appropriate conditions, equivalent to evolution 
under continuous measurement in the form of a complex step potential. Before we can do 
this we will need some results about propagation in the presence of these step potentials and 

^Thc; iciason we have not discussed POVMs is tlic following. Consider two distinct histories, {a, /?}, defined 
by strings of projection operators. Now turn the projection operators into POVMs by smearing them 
with a gaussian of finite width and call these new histories {q,/3}. These smeared histories are no longer 
distinct and so even if D{a, /?) = we will have D{a, ^) ^ 0. This makes it very hard to decide between 
histories that interfere, and histories that simply overlap. 

^In fact in some ways we do exactly this in later sections, since the complex potential introduced in Chapters 
4 and 5 is very much like a POVM. Our motivation for introducing it there is rather different, however. 



31 



Chapters 2 and 3 will be devoted to explaining the necessary path integral methods. 

Complex potentials also arise in other standard approaches to the arrival time problem 
in classical and quantum mechanics. These potentials are generally of the form V{x) = 
—iVox{x), where xi^) is the characteristic function of some region of configuration space. If 
we are to work with these complex potentials we must be able to deal with the propagator 
between two arbitary points in configuration space in the presence of these simple complex 
potentials. In Chapter 2 we introduce the Path Decomposition Expansion (PDX), a useful 
technique for evaluating path integrals with piecewise defined potentials. In Chapter 3 we 
then use this technique to evaluate the propagator for the step and Dirac delta function 
potentials. As well as introducing the full expression for the PDX, which is exact, we also 
introduce a useful semi-classical approximation, valid when the height of the potential is 
small compared with the energy of the incoming state. 

In Chapter 4 we ask the following question: Are the schemes of pulsed and continuous 
measurement, represented by complex potentials and strings of projection operators respec- 
tively, in any sense equivalent? We will find that the two approaches are indeed related, 
in the sense that the propagators for evolution under the two measurement schemes are 
equivalent under certain conditions. 

In Chapter 5 of this thesis we use the results of Chapters 2, 3 and 4 to attack the arrival time 
problem. Once we have proven that evolutions under pulsed and continuous measurements 
arc equivalent wc use this to replace the class operator, Eq.(1.46) with evolution in the 
presence of a complex potential. In Chapter 5 this will allow us to undertake a thorough 
examination of the arrival time problem and to obtain the arrival time distribution n(t). 

However, despite the apparent complexity of the approach, the final results will be re- 
markably simple and will suggest that there may be a less rigorous but more intuitive way 
to arrive at them. In Chapter 6 wc will therefore present a simplified derivation of the class 
operators for the arrival time problem, using only a semiclassical approximation. We will 
also show how these class operators may be extended in a straightforward way to situations 
where there are states incident on the origin with both positive and negative momenta. 

With the main goal of this thesis achieved, in Chapters 7 and 8 we turn to some related 
questions that lie somewhat outside the main scope of this work. In Chapter 7 wc look 
at arrival times for open quantum systems, trying to understand the way in which the 
current, Eq.(1.20), emerges as the classical arrival time distribution, both in general terms, 
and by extending the decoherent histories analysis to this situation. In Chapter 8 we turn 
to the issue of arrival and dwell times defined via ideal clocks, finding agreement with the 
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decoherent histories analysis in the appropriate hmit, but also achieving an understanding 
of the emergence of the current as the correct arrival time distribution for a very general 
class of model clocks. 

We summarize the results of this thesis in Chapter 9, and outline some possible areas for 
further study. 
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2 The Path Decomposition Expansion 



In this chapter we describe some useful path integral techniques. We will make frequent 
use of the results in this chapter throughout the rest of this thesis. Since this chapter is 
essentially technical background, the majority of the material presented here is not original 
research. In Section 2.2 we introduce the Path Decomposition Expansion, then in Section 2.3 
we use this to compute the scattering states for a complex step potential. In Section 2.4 we 
then introduce a useful semiclassical approximation, valid when the height of the potential 
is small compared with the energy of the incoming state. 

2.2 The Path Decomposition Expansion (PDX) 

In this section we introduce the Path Decomposition Expansion (PDX), a useful technique 
for evaluating path integrals with piecewise defined potentials. Although for the most part 
the application of these results will be to the simple situation where we have a particle 
incoming on an imaginary potential of step function form, in this chapter we will deal with 
real potentials and we will begin by assuming the more general form V{x) — Vo9{—x)f{x). 




There is surely no greater wisdom than 
to mark well the beginnings and endings 
of things. 



Francis Bacon, Essays 1625, Of Delays 



2.1 Introduction 
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We wish to evaluate the propagator 



g{xi,T\xo,Q) = {xi\exp{-iHoT - iVo9{-x)f{x)T) \xo), (2.1) 
for arbitrary xi and Xq > 0. This may be calculated using a sum over paths, 

g{xi,T\xo,0) = J Vxexp^iS), (2.2) 

where 

']mx''-Voe{-x)f{x)) , (2.3) 



2 

and the sum is over all paths x(t) from a;(0) = Xq to x{t) — Xi. 

To deal with the step function form of the potential we need to split off the sections of the 
paths lying entirely in ,t > or x < 0. The way to do this is to use the path decomposition 
expansion or PDX [8, 9, 48, 42]. Consider first paths from xq > to Xi < 0. A typical such 
path may cross x = many times, but the set of paths may be partitioned by their first or 
last crossing times. We therefore split every path into three parts: (A) a restricted part that 
starts at Xq and does not cross x — 0, but that ends on a; = at time ti, (B) an unrestricted 
part from a; = to x = that may cross x = many times and, (C) a further restricted 
part from ,x = to xi that does not re-cross x — 0, Fig. (2.1). As a consequence of this, it is 
possible to derive the formula, 

i r dgr 
g{xi,T\xo,Q) = dti g{xi,T\Q,ti)-^{x,ti\xo,0)\^^^. (2.4) 

Here, gr{x,t\xo,Q) is the restricted propagator given by a sum over paths of the form (2.2) 
but with all paths restricted to x{t) > 0. It vanishes when either end point is the origin but 
its derivative at x = is non-zero (and in fact the derivative of gr corresponds to a sum over 
all paths in X > which end on x = [48]). 

It is also useful to record a PDX formula involving the last crossing time t2, for xq > 
and xi < 0, 

i r dgr 

g{xi,T\xo,Q) = -— J dt2 -^{xi,T\x,t2)\^^^^ giO,t2\xo,0). (2.5) 
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Figure 2.1: A typical path from xq to xi 



These two formulae may be combined to give a first and last crossing version of the PDX, 

1 Z"*^ dg dg 

g{xi,T\xo,0) = ^J dt2 j dti ^(xi, t|x, ^2) |^=o ^(0,^2|0,ii) -^{x,ti\xo,0)\^^Q. 

(2.6) 

This is clearly very useful for a step potential since the propagator is decomposed in terms 
of propagation in 2; < and in x > 0, essentially reducing the problem to that of computing 
the propagator along x — 0, g{0,t2\0,ti). (Sec Figure 2.1) 

For paths with xq > and xi > 0, Eq.(2.4) is modified by the addition of a term 
gr{xi, t|a;o, 0), corresponding to a sum over paths which never cross x — 0, so we have 

i r dgr 
g{xr,T\xQ,Q) = dti g{xi,T\0,ti)-^{x,ti\xo,0)\^^^^ + gr{xi,t\xo,0). (2.7) 

Again a further decomposition involving the last crossing, as in Eq.(2.6) can also be included. 
In general the usefulness of these expressions relies on the various partial propagators that 
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occur in the PDX formulae being easier to compute than the full propagator. This depends 
crucially on the form of the potential, from now on we will assume that f[x) = 1 in the 
above expressions, and thus that the potential is of simple step function form. 

The various elements of these expressions are easily calculated for a potential of simple 
step function form V{x) — Vo9{—x). The restricted propagator in x > is given by the 
method of images expression 

gr{xi,T\xo, 0) = e{xi)e{xo) {gf{xu t\xo, 0) - gf{-xi, t\xo, 0)) , (2.8) 

where gf denotes the free particle propagator 

gfix„r\xo,0) = [—) exp j . 



(2.9) 



It follows that 

^(^,tiKO)|_, = 2^(0,ti|a;o,0)^(a;o). (2.10) 

The restricted propagator in a; < is given by an expression similar to Eq.(2.8), multiplied 
by exp(iVor). 

Note that this means that Eqs.(2.4), (2.5) and (2.6) can be written as, 

(xi|e-^(^°+^°^(-*))Mxo) = - rciii(xi|e-*(^°+^°^(-*))("-*i)(^(x)pe-^^°*Mxo) (2.11) 

m Jo 

= - r dt2 {xi\ e-^(-f^o+^°)(^-*2)p5(£)e-'^^°+^°^(-*))*^ |a;o)(2.12) 
m Jo 



^ \ f dt2 [\ti{xAe-'^"°+^°^^^-^^^p5(x) 
Jo Jo 

^^-i(Ho+Vo9{-mt2-h)§^^^P^-iHoh 1^^^ (2.13) 

where 6{x) = |0) (0| and |0) denotes a position eigenstate \x) at a; = 0. These operator forms 
of the PDX are the ones we shall use most often. 

The only complicated propagator to calculate is the propagation from the origin to itself 
along the edge of the potential. We will show in the next chapter that in the case V{x) = 
Vo9{—x) this is given by [99], 



37 



However as well as allowing us to obtain exact results for the propagator in certain cases, 
the PDX also suggests some semiclassical approximations we could make to simplify calcu- 
lation for a more general class of potentials. We will discuss this in Section 2.4. 



2.3 Using the PDX: Scattering States for a Complex 
Step Potential 

In this Section wc use the PDX to derive the standard representations of the scattering 
solutions to the Schrodinger equation with the simple complex potential V{x) = —iVo9{—x). 
These are known results but this derivation confirms the validity of the PDX method and 
also allows a certain heuristic path integral approximation to be tested. The results will also 
be useful for the decoherent histories analysis in Chapter 5. 
The transmitted and reflected wave functions are deflned by 

i^ix, r) = 9{-x)iptr{x, r) + 9{x) {ipref{x, t) + 'll)f{x, t)) (2.15) 

Here, ■0^^ is the transmitted wave function and is given by the propagation of the initial state 
using the PDX formulae Eq.(2.4) or Eq.(2.6). For large r, a freely evolving packet moves 
entirely into x < so that the free part, r) is zero, leaving just the transmitted and 

reflected parts. Following the above deflnition, the transmitted wave function is given by 

il^tr{x,T) — — - I ds I dv (x| exp (— iifos)p|0) e~^°^ 
^ Jo Jo 

X {0\exp{-iHv)\0) {0\pexp{-iHo{T -V - s))\tlj), (2.16) 

where |0) denotes the position eigenstate \x) at x — 0. Also, we have introduced s — t — ti 
and V — t2 — t-i, and H — Hq — iVQ9{—x) is the total Hamiltonian. The scattering wave 
functions concern the regime of large r, so we let the upper limit of the integration ranges 
extend to oo. 

Writing the initial state as a sum of momentum states and introducing E — p^/2m, 
we have 

i^trix.T) = \ I dp I (is(a;|exp(-ii/os)p|0) e^^^+^^°)^ 
^ J Jo 

X / (0|exp(-iii'v) |0) e^-^XOlp) e-^^^V'(p). (2.17) 
Jo 
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To evaluate the s integral, we use the formula [92], 



^0 



Xs + 



(S)' e^p(^l^l^2;^)' (2-18) 



2s 

from which it follows by differentiation with respect to x and setting X — E + iVq that 

/oo 
(is(a;|exp(-ii/os)p|0) e^^^+'^o^" = mexp (^i|a;|[2m(S + j . (2.19) 

The V integral may be evaluated using the explicit expression for the propagator along the 
edge of the potential, Eq.(2.14), together with the formula. 

We thus obtain the result, 

^t,{x,T)^ J -^exp(^-ix[2m{E + iVo)]'^ -iEr^ ^t,(p), (2.21) 

where 

Mp) = ^ .^^,1, i^ip). (2.22) 

{l + E-2{E + iVo)2) 

Note that in this final result, it is possible to identify the specific effects of the different 
sections of propagation: the propagation along the edge of the potential corresponds to the 
coefficient in the transmission amplitude Eq.(2.22) (which is equal to 1 when Vq = 0), and 
the propagation from final crossing to the final point produces the Vq dependence of the 
exponent. These observations will be useful below. 

The refiected wave function i/jref is defined above using the PDX Eq.(2.7) (rewritten using 
Eq.(2.6). The first term in Eq.(2.7), the crossing part, is the same as the transmitted case, 
Eq.(2.17), except that Vq = in the last segment of propagation, from x = to the final 
point, and also the sign of x is reversed. We must also add the effects of the refiection part 
of the restricted propagator, and this simply subtracts the refiection of the incoming wave 
packet. The refiected wave function is therefore given by, 

/dp 
-= exp (ixp - i^Jr) ipref{p), (2.23) 
y ZTT 
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where 



1- E-'2{E + iVo 



[l + E-"2{E + iVoy^) 



^ip). (2.24) 



We thus see that the PDX very readily gives the standard stationary wave functions [6] , 
without having to use the usual (somewhat cumbersome) technique of matching eigenfunc- 
tions at x = 0. In fact, this procedure is in some sense already encoded in the PDX. 



2.4 A Semiclassical Approximation 

In this section we discuss a useful semiclassical approximation that allows us to obtain 
the approximate propagator for evolution under simple step potentials of the form V{x) = 
Vo9{—x) for Vo real or imaginary, under the condition that |Vo| is small compared with the 
energy of the system. 

Pictorially, the PDX equation Eq.(2.6) says that propagation from xo > to xi < may 
be considered in three parts, initial propagation to x — 0, propagation along x — and 
propagation away from x = 0, see Fig.(2.1). Classically however, assuming we can neglect 
reflection, we expect the propagator to be dominated by the single classical path. This path 
is a simple straight line that crosses x — only once. One might therefore expect that the 
propagator for the final segment from x = to Xi < in Eq.(2.4) is given approximately by 

{xi\ exp {-iHos - Vo9{-x)s) |0) ^ {xi\ exp {-iHos) |0) exp {-Vqs) (2.25) 

and similarly that, 

(0| exp {-iHos - Voe{-x)s) \xo) ^ (0| exp {-iHos) \xo) (2.26) 

It is not entirely clear that this is the case, however. On the one hand, the usual semiclassical 
approximation indicates that paths close to the straight line paths dominate, but on the other 
hand, paths in x < are suppressed, so maybe the wiggly paths that spend less time in a; < 
make a significant contribution. Since this approximation is potentially a useful one, it is 
useful to compare with the exact result for the transmitted wave packet calculated above. 
We therefore evaluate the following approximate expression for the transmitted wave func- 
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tion, 

i;tJx,T) ^ - I ds (xle-'^^'lO) 6'^°' (0\pe-'^°^^-'^m (2.27) 
Jo 

This is the PDX, Eq.(2.4), in operator form with the semiclassical approximation described 
above and we have set s = r — ti. We now take r — )■ cxd in the integration and evaluate. The 
key integral is, 

r ds (a;|e-'^»^|0) e^^^+'^o^* = ( ^.^"^ ' exp ( -ix[2m{E + iVo)]'A (2.28) 
Jo \2[E + iVo)J \ / 

where we have used Eq.(2.18) (and recall that a; < 0). The resulting expression for the 
transmitted wave function is of the form Eq.(2.21), with 

Arip) = -—r-!——Tm (2.29) 

This agrees with the exact expression for the transmission coefficient Eq.(2.22) only when 
Vq = 0, with the difference of order Vq/E for small Vq. This establishes that the approxima- 
tion is valid for Vq much less than the energy scale of the initial state. 

Crucial to the usefulness of this approximation, however, is that the form of the trans- 
mitted wavefunction is reproduced exactly, the discrepancy being in the prefactor. In the 
situations where we make use of this semiclassical approximation we are always assuming re- 
fection may be neglected, and thus the prefactor is in fact irrelevant, since it can be obtained 
by normalization. This makes this semiclassical approximation a very useful one. 

We make some brief comments here about the significance of the semiclassical approxi- 
mation in the arrival time problem. The classical phase space distribution representing the 
current, 

J(p,x) ^^6{x) (2.30) 
m 

does not have a unique quantisation, because of the ambiguity in the operator ordering of p 
and S{x). The standard Schrodinger current comes from the symmetric choice 

i = ^ mx) + 5ix)p) . (2.31) 
Now in the PDX one can write the propagator from a point > to a point Xi < in 
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either of the forms 

(xi|e-*^^|xo) = - / dt {xi\e-'"^^-'^6{x)pe-'''' \xo) (2.32) 
m Jo 

= - / dt(a;i|e-^^('^-*)p(5(£)e-^^*|xo). (2.33) 
m Jo 

These two integrals are equal, but the integrands are not. They differ in in terms of whether 
it is the incoming or outgoing paths which are restricted. However the semiclassical approxi- 
mation says that the integrands, thought of as path integrals, are dominated by paths which 
cross only once and are therefore in effect "restricted" on both sides. This means we have, 

{xi\ e-^^(^-*)5(i;)pe-'^* |a;o) ^ e-'^(^-*)p5(x)e-^^* |a;o) . 

Because the times and positions are arbitrary this approximation means that there is essen- 
tially no difference between the terms pS{x) and 5{x)p. This means the semiclassically any 
result obtained by using the first crossing PDX, say, could equally well be obtained using 
the last crossing formula. 
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3 The Propagator for the Step and 
Delta Function Potentials, using the 
PDX 

The journey of a thousand miles begins 
with one step. 

Lao Tzu 

3.1 Introduction 

One of the standard approaches to defining the time of entry to a given region of configuration 
space in quantum theory is to introduce a complex potential localised in the region of interest, 
and to define the arrival time probabilities in terms of the rate at which the wavefunction 
is absorbed by the potential [6]. In the simplest case of the arrival time problem in ID 
this means using a complex potential of the form V{x) = —iVo6{—x). If we are to extract 
the arrival time probability from this set up we must be able to solve for the propagator 
in the presence of simple complex potentials. In fact, close inspection reveals that the 
propagator for a particle in the presence of a complex potential can generally be obtained 
by first considering the equivalent real potential, and then analytically continuing the result. 
This means our task is to consider propagators for a particle in the presence of a real step 
potential. Delta function potentials may be used in a similar way. 

Calculations involving step and delta function potentials also occur in many other branches 
of physics. Step potentials can be used to represent "hard wall" boundary conditions, and 
are also involved in tunneling calculations. Delta function potentials can be used to model 
point interactions, especially in the low energy limit where the details of the process are 
largely independent of the form of the scattering potential [5] . 
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We wish to compute the following propagator, 

r'x{T)=xi 



g{xi,T\xo,0) ^ / Pxe^^ (3.1) 

J x{0)=xo 



where 



T / -1 

' mx 







and the potential ^(.t) is cither a step potential \Q(^—x\ or a delta function potential, A b[x\ 
In this Chapter we present a derivation of the propagators in the vicinity of these potentials 
by making use of the path decomposition expansion presented in the previous chapter. The 
problem of calculating the full propagator therefore reduces to that of calculating the partial 
propagator for the interval where the path crosses the origin, (^(0, t2|0, ti). We will show how 
these partial propagators may be derived using the Brownian motion definition of the path 
integral [57, 89]. The full propagators may then be obtained with the help of Eq.(2.6). 

For the step potential the full propagator has been derived in Refs.[15, 11, 16, 17], but 
we shall derive the partial propagator, and then direct the reader to Ref.[15] for details of 
the use of the path decomposition expansion to recover the full propagator. The partial 
propagator we will derive is given by 

9(0.T|0,0)^-,(i^) (3.3) 

This result was used in Chapter 2 to compute the scattering states for a complex step 
potential. 

For the delta function potential we will derive the full propagator and we quote the result 
from Ref.l 



/•oo 

(7(a;i,T|a;o,0) = (7/(a;i,T|a;o,0)-Am / dii e-^^'XIxil + |a;o| + T|0, 0), (3.4) 



where 

m \ 1/2 



^,(^i,T|xo,0)=(^) e-(---)V2^ (3.5) 



is the free propagator. This chapter is based on Ref.[99]. 
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3.2 The Brownian Motion Definition of the Propagator 



We begin with a review of some of the details of the Brownian motion approach to computing 
propagators. For more details see Refs.[57, 89]. The first step is to switch to working with 
the Euclidean propagator g by means of a Wick rotation and we specialise immediately to 
the case of xq = xi — 0. That is, we wish to calculate 

^(0,T|0,0) = / Vx e-^^ (3.6) 

Jx{0)=0 

where Se is the Euclidean action given by 



This propagator may be viewed as a conditional probability density for a random walk on 
the real line. The second step is then to make this integral over paths into a concrete object 
by defining it as the continuum limit of a discrete sum on a lattice. 

To establish conventions and demonstrate the basic ideas we compute the case of a free 
particle, following closely the treatment in Ref.[57]. We consider a rectangular lattice with 
spacing in the time direction of e, spacing in the x direction of r] and we consider propagation 
for a time T = 2en, so we have 2n steps in our paths (the reason for this choice is that it 
simplifies a number of later expressions, and avoids clumsy factors of 1/2). The conditional 
probability ■u(0,T|0,0) to start at (0,0) and end at (0,T) is given by the number of paths 
connecting the start and end points, divided by the total number of possible paths. The 
set of paths from (0,0) to (0,T) is bounded by the extremal paths that take n steps to the 
left/right, followed by n steps to the right/left, see Fig. (3.1). Since a path must have the 
same number of steps to the right as to the left to end up back at a; = we find. 



The Euclidean propagator g is then defined as the continuum limit oi u/2ri where we take 
e, ?7 — )■ 0, n — )■ oo, keeping e/rj^ = m and T = 2en fixed. That is. 




(3.7) 






(3.9) 
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which is the expected result for the Euchdean free propagator. 




Figure 3.1: typical path from x = to x — 0. The dashed line shows the area within which 
all paths must remain. 



3.3 The Step Potential 

The propagator along the edge of a step potential is given by Eqs.(3.6) and (3.7), with 
V{x) — \9{—x). We can write this as 

^(0,T|0,0)= J Vxexp(^- dt^^^exp(^-X dt9{-x)^ (3.10) 

which is similar to the free particle case except that paths are weighted by a factor e~^'^ 
where 

T 

dte{-x) (3.11) 

is the length of time spent in x < 0. In the lattice case the corresponding conditional 
probability u{0,T\0,0) is given by a sum of paths, each weighted by a similar factor. This 
may be written as 

1 " 

u{0,T\0,0) = ^ke-'''\ (3.12) 

fe=0 
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where is the number of paths spending a time 2ke < T in the region a: < 0. Expresions 
for these Uk are known and are in fact independent of k [70] . They are equal to the Catalan 
numbers 

n+l\n J 

where 2n is the total number of time steps. We can see this as follows. First, note that 
the number of paths that never enter x < is given by C„, this being one definition of 
the Catalan numbers. Next consider the following mapping on any path spending a time 
2ek<T inx < 0, Fig.(3.2). 

1. Start from t=0, and follow the path until it first crosses x = (if it doesn't cross then 
stop, the path is in the set of non-crossing paths.) 

2. Follow the path until it comes back to x — again, note the step at which this happens. 

3. Swap the section of path before this step with the section after it. 

4. The new path will now spend 2 fewer time steps in x < 0. 

By repeated application of this mapping, any path can be transformed into one which never 
crosses a; = 0. The important point about this mapping however, is that it is bijective [70], 
which proves that the number of paths spending time 2ke in a; < must equal the number 
of paths that never cross, for any value of k. This shows that rik = Cn for all k < n. So we 
now have 



u 



^ n ^ n 

(O'^IO'O) - ^E^"-"'^' = ^E^"'^' (3-14) 

fc=0 k=Q 

^ (3.15) 



22n \^ l_ g-2eA 

This coincides with the free particle case if A = 0. Since we plan to take the continuum limit 
we can Taylor expand the exponential in the denominator to first order in e, and use the 
following useful assymptotic form for C„ [70] 

An 

(J ~ 

to get, 

I _ g-2e(n+l)A 

M(o,r|o,o) 



v^2en3/2 \^ A 
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Figure 3.2: Application of the rule to the path on the left produces the path on the right, with 
2 timesteps less spent in x < 



Now take the continuum limit as in Section 3.2 to obtain, after some simple algebra, 

5(0,r|0,0)=(-) 4^ (3.16) 

which implies 

, ^ / m \i/2 (1 - e-^^^) 

,(0.T|0.0) = -.(-) (3.17) 

This is our first result, the propagator along the edge of a step potential. The full propagator 
from xq to xi may now be obtained by making use of Eq.(2.6) [15]. Note also that this result 
may be analytically continued in a simple way to imaginary potentials. 



3.4 The Delta Function Potential 

We now wish to evaluate Eqs.(3.6) and (3.7) with the potential given by V{x) — X6{x). In 
a similar way as for the step potential, we can write this as, 

^(0,T|0,0)= J Vxexp(^- dt^^^exp(^-\J^ dtS{x{t))^ (3.18) 
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which is again similar to the free particle case except that paths are weighted by a factor 
e~^"' where 



(7= / dtS{x{t)) (3.19) 
Jo 



is the number of times a given path crosses x = 0. We model the delta function as a square 
potential of width r), and height l/r] so that every crossing of x = is weighted by a factor 
of e~^^^'^. Since = m, we can rewrite this as e~^^"*. We take the potential to be located 
on the left of x = so that to start and end at x = involves an even number of crossings. 
Such choices are only made for convenience and have no significance in the continuum limit. 
The conditional probability density u may be partitioned in a similar way to that for the 
step potential. However now the partitioning is with respect to the number of times a path 
crosses the square potential. Since the number of crossings will always be even we partition 
into classes of paths that cross 21 < 2n times, so the conditional probability is 

1 " 

u{0, T|0, 0) = — J] J(n, /)e-2'-^^ (3.20) 
1=0 

Where J{n, I) is the number of paths of 2n steps that cross the delta potential 21 times. It is 
known that these J{n, I) are given by the 2/*^* convolution of the Catalan numbers [96], this 
can be demonstrated by writing a general path in terms of sums over non-crossing paths. 
These convolutions form the diagonal elements in Catalan's triangle [96, 86, 94]. We need 
to know the n^'^ element in the 2/*^* diagonal from the right, which will give us J{n, I). 
From the formula for the elements of Catalan's triangle [96], 

_ {n + k)\{n-k + l) 

k\{n + l)\ ^"^-^^^ 

it follows that 

,x / , ,x (2n)!(2Z + l) f2n\ n\n\(2l + 1) 

J{n,l) = c{n + l,n-l) = . ^ / = 7 ^77 T^TTT 3-22 

^ ^ ^ ^ {n - l)\{n + 1 + 1)\ \n J {n-l)\{n + l + l)\ ^ ' 

where we have extracted the binomial factor for later convenience. We cannot perform the 
summation in Eq.(3.20) directly as we did for the step function potential, so we take the 
continuum limit first to leave ourselves with an integral. In order to do this we need a 
simpler form for the J's. We need the asymptotic form for n — >■ 00, but plotting Eq.(3.22) 
as a function of / shows that the dominant contribution comes from taking / ^ 00 as well. 
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It is possible to derive the following asymptotic form, 



J(n,0~r")-e-^^^ (3.23) 

and as this form shows, the maximum value of J occurs when / is of order ^/n. This is 
a consequence of the Brownian motion origin of the paths, which therefore have Hausdorff 
dimension 2 [1]. The statement that a typical path crosses the origin an infinite number of 
times is a consequence of this fractal nature of a typical path. 

If we use the asymptotic form for the J's, Eq. (3.23) and make the change of variable in 
the partitioning, Eq. (3.20), u — 2r)l we can turn the summation into an integral, 

1 /2r?\ Z'^''" V 
«(0.r|«'.0) = ^( JX Wn''"^-'""''- 

We now take the continuum limit as in Section 3.2 to obtain 

^(0,r|0,0) = (^^j duue-'^e-'"^^ 

POD 

= ^^(0,T|0,0) - Am / rfw e-^™"^^(M,T|0,0) (3.25) 

Jo 

which, after rotating back to real time, yields 

POD 

^(0,T|0,0) = c//(0,T|0,0) - Am / rfw e"^™"c//(M, T|0, 0). (3.26) 

Jo 

We can now use the path decompsition expansion to obtain the propagator for paths which 
start at xq and end at Xi. There are 4 cases, depending on the signs of Xq and Xi, we 
shall present one case here, the others follow in a very similar fashion. First note that PDX 
formulae in Chapter 2 tell us that the free propagator obeys the following. 



gf{xi,T\xo,0) = — / dti gf{xi,T\0,ti) —{x,ti\xo,0) 
m Jq ox 



(3.27) 



x=0 



i , dgt . ^. . 
— / dt2 -^[xi,T\x,t2) 
m /q ox 



(7/(0,^2^0,0) (3.28) 



x=0 



provided e(a;o) = e{—xi). (where gf denotes the free propagator, and t{x) is the signum 
function.) (If the start and end points have the same sign then there are paths between the 
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two that never cross a: = 0, so the expression for the path decomposition expansion has to 
be modified as in Eq.(2.7) [48].) 

Secondly note the following identities, 



gf{a,t2\b,ti) = gf{-a,t2\ - b,ti) = gf{Q,t2\b - a,ti) = gf{Q,t2\a - b,ti). (3.29) 

which are just expressions of the symmetry of the free propagator. 

Consider the case where xq < and Xi > 0, noting that u > we first use Eq. (3.27) to 
attach the leg from (xq, 0), 

poo 

^(0,^2^0,0) = gf{0,t2\xo,0) - \m due-^"'''gf{u,t2\xo,Q) 

Jo 

poo 

= gf{0,t2\xo,0) - Xm due'^'^''gf{0,t2\xo - u,0) (3.30) 

Jo 

where we have used Eq. (3.29) to obtain the second line. Since Xq < we have that 
xq — u < 0, so we can attach a leg to (xi, T) using Eq.(3.28), 

roo 

g{xi,T\xo,0) = gf{xi,T\xo,0) - Xm due-^'^''gf{xi,T\xo - u,0) 

Jo 

poo 

= gf{xi,T\xo,0) - Xm due-^"'''gf{\xi\ + \xo\ + u,T\0,0) (3.31) 

Jo 

as expected. Calculation in the other cases proceeds in a very similar manner, and confirms 
that Eq.(3.31) is valid in all cases. 
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4 On the Relationship between 

Complex Potentials and Strings of 
Projection Operators 



Time's a strange fellow; 

more he gives 

than takes 

(and he takes all)nor any marvel finds 
quite disappearance but some keener 

makes 

losing, gaining 

-lovelif a world ends 
more than all worlds begin to(see?)begin 

E.E. Cummings 



4.1 Introduction 

In Chapter 1 wc argued that in dccohcrcnt histories the class operator describing histories 
which remain in some region defined by a projector P is given by C^^, where, 

C'J^p) = e~'"'Pe~'"'...Pe-'"'\ij) (4.1) 

(Here, there are n projection operators and t — {n + l)e and we use units in which ^ = 1). 
Such an object also describes the evolution of an initial states under pulsed measurements. 

We are interested in the specific case of a free particle with P taken to be the projection 
onto the positive x-axis, P — 9{x). For sufficiently small e, Eq.(4.1) is then a candidate for the 
amplitude to remain in x > during the time interval [0, r]. In this Chapter we examine the 
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object in Eq.(4.1) without particular reference to it's relavence in the arrival time problem. 
Our aim will be to show that evolution under Eq.(4.1) is equivalent to evolution under a 
particular form of continuous measurement. We will then use this result in Chapter 5 to 
address the arrival time problem. 

It is of interest to explore the properties of this Eq.(4.1) for a range of values of the 
time spacing e. It is known that as e — > 0, we approach the Zeno limit, in which the state 
becomes entirely confined to the Hilbert subspace of states with support only in a; > 0, so 
that an incoming wave packet from the right is totally reflected [74, 95, 27]. However, it is of 
greater physical interest to explore the regime of non-zero e, in which the system is monitored 
sufficiently well to get some idea of whether the particle is in x > 0, yet not monitored so 
much that an incoming state is significantly reflected at a; = 0. An important question in 
this regime is to determine the value of e for which reflection becomes significant. For an 
initial state with energy width Ai?, a timescale held to be significant is the Zeno time, 

which is the timescale on which the state becomes significantly different from its initial 
value under unitary evolution [87]. For a wave packet of momentum p and spatial width cr, 
the Zeno time is of order ma/p which is the timescale on which the wave packet crosses the 
origin. This indicates that the Zeno time for wave packets is an essentially classical timescale 
and, in Eq.(4.1), relates only to the rate of removal of probability through projection. By 
contrast, reflection in Eq.(4.1) arises as a result of the increase in uncertainty in momentum 
resulting from position projection, an obviously quantum process, so one would expect it to 
have a different timescale, which could be much shorter than then Zeno time. It would be 
of interest to compute this timescale. One reason it is important is that there appears to be 
interesting physics very close to the Zeno hmit [23, 82, 45]. 

A signiflcant result in this area is due to Echanobe, del Campo and Muga, who claimed 
that for finite e the string of projection operators in Eq.(4.1) is approximately equivalent to 
evolution in the presence of a complex potential [23]. That is, 

^-iHep^-iHe . . . p^-iHe _ {-iHt - Voe{-x)T) (4.3) 
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They assert that this result is vahd if, for a given e, Vq is chosen such that two conditions 



Voe > 1 (4.4) 
Vo > /^H (4.5) 

arc satisfied. This is a very useful result since Eq.(4.1) is not easy to evaluate analytically but 
the Schrodinger equation with a complex step potential in Eq.(4.3) can be solved straightfor- 
wardly. Furthermore, such complex potentials have been studied extensively in the literature 
and can often be linked to particular detection methods [76]. 

Given the connection Eq.(4.3), one can determine the conditions under which reflection 
becomes important. Known results on scattering with the complex potential in Eq.(4.3) 
show that, for an incoming state with energy scale reflection becomes significant when 
Vq> E [Q, 49]. Reflection is avoided, therefore, when Vq <^ E, or equivalently, from Eq.(4.4), 
when 

e » ^ (4.6) 

This is much less than the Zeno time for a state strongly peaked in energy. There is therefore 
an interesting regime, namely 

in which the projections in Eq.(4.1) are sufficiently frequent to have a significant effect on 
the system, yet not so frequent that there is significant refiection. 

Eq.(4.6) is a very useful result, but it has been derived on the basis of the claimed ap- 
proximate relationship Eq.(4.3). The derivation of Eq.(4.3) given by Echanobe et al is very 
plausible (and was also hinted at by Allcock [6]), but it is rather heuristic, and the deduced 
connection Eq.(4.4) between Vo and e is rather loose. There is therefore considerable scope 
for a more detailed and substantial derivation. 

The purpose of this chapter is to give a more substantial derivation of the equivalence 
Eq.(4.3) and to deduce a more precise relationship between Vq and e. We will do this by 
computing the configuration space propagators associated with each side of Eq.(4.3) and 
show that they are approximately equal in certain regimes. The propagator associated 
with the complex potential is in fact known already, so the bulk of the work consists of an 
approximate evaluation of the propagator associated with pulsed measurements, the left- 
hand side of Eq.(4.3). 

This is certainly not a rigorous mathematical proof of Eq.(4.3), involving operator norms. 
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error bounds and the like, but a theoretical physicists style of proof involving the approximate 
evaluation of propagators. A rigorous proof would certainly be of interest to construct and 
the work described in this chapter may give some hints in that direction. 

We begin in Section 4.2 with a brief summary of the derivation of Echanobe et al. with a 
small extension of it that turns out to be important and yields an equality relating Vq and 
e, thereby improving on Eq.(4.4). Wc then in Section 4.3 give a detailed formulation of the 
problem and how we solve it. The key idea is to use the path decomposition expansion in 
which the propagators associated with each side of Eq.(4.3) are factored across the surface 
a; = [8, 9, 48, 42]. The problem of proving the equivalence Eq.(4.3) thereby reduces to 
proving it for propagation between points at x = for different times. Since the propagator 
for the complex potential is known, the main work is to evaluate this propagator for pulsed 
measurements. This is actually rather difficult to do directly, but a good approximate 
analytic expression can be obtained for it by attacking the problem from a number of different 
angles. This is described in Sections 4.4, 4.5, 4.6 and 4.7. In Section 4.8 we give a detailed 
discussion of the timescales involved for the approximations to be valid. We summarize and 
conclude in Section 4.9. This chapter is based on Ref.[51]. 

4.2 Review and Extension of Earlier Work 

We first review and extend the derivation of Echanobe et al [23] . They first note that 





where, recall, P = 9{x) and P = 1 — P. It follows that 




(4.9) 



as long as the parameter 



a — Vo€ 



(4.10) 



is sufficiently large that 



e-" < 1. 




Eq.(4.3) then follows from the approximate equivalence. 



exp {—iHe) exp (— Vo^(— x)e) ^ exp {—iHe — Vo^(— x)e) , 



(4.12) 
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which will hold as long as 

Voe'\{[H,9{-x)])\<-h (4-13) 

where the average is taken in the initial state. Echanobe et al put an upper bound on the 
left-hand side using the Schrodinger- Robertson inequality and Eq.(4.13) may then be written 
in either of the two equivalent forms 

a^AH < Vo (4.14) 



or 

ae « — (4.15) 

which implies that the time between projections is much less than the Zeno time, the typical 
timescale on which the state undergoes significant change [87]. (The conditions Eqs.(4.10), 
(4.11), (4.14) are a more precise version of the originally stated conditions Eqs.(4.4), (4.5)). 

Since the parameter a need only satisfy an inequality, Eq.(4.11), the relationship Eq.(4.10) 
between e and Vq is not uniquely determined. However, it turns out that the above derivation 
can be extended somewhat to give an equality between e and Vq. This turns out to be relevant 
to the more substantial derivation given in the rest of this chapter. 

The above result may be written, 

exp {-iHe - Vo0{-x)e) ^ exp {-iHe) P (4.16) 
under the conditions given above. Now suppose a is an integer and write e = ae', where 

Voe' = 1 (4.17) 
Then, since P = P^, we may approximate Eq.(4.16) by 

exp {-iHe - Voe{-x)e) ^ e''"'' Pe''"'' P ■ ■ ■ e''"'' P 



(e-^^'py 



(4.18) 



as long as the contribution from the commutator terms between e~*^^' and P is sufficiently 
small. There will be of order such terms, hence the error in this approximation is of order 
a^e'AH and from Eq.(4.14), this error is much less than 1. The key point here is that even 
the longer timescale e is still much less than the Zeno time and, since nothing changes on 
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this timescale, there is essentially no difference between Eqs.(4.16) and (4.18). 

Wc therefore see that the desired result Eq.(4.3) actually holds for much smaller time 
steps, defined by the equality Eq.(4.17). This is different to the original claim of Echanobe 
et al, since they require the inequality Eq.(4.4). However, the above argument shows that 
the restriction Eq.(4.4) is in fact stronger than necessary and in the following pages our more 
detailed derivation will show that Eq.(4.3) does indeed hold with a timespacing of order I/Vq. 



4.3 Detailed Formulation of the Problem 



In this chapter, we will prove the relationship Eq.(4.3) in a much more substantial way by 
proving the approximate equivalence of the propagators 



gv{xi,T\xo,0) = {xi\exp{-iHT -Vq9{-x)t)\xo) 
gp{xi,T\xo,0) = {xi\e-'"'"Pe- 



-iHe g-iHe 



Pe 



-iHeo 



Xo) 



(4.19) 
(4.20) 



for some relationship between the parameters e and Vq, to be determined. Note that in 
Eq.(4.20), we have chosen the initial and final time spacings to be eo and e„, with 



T = (n - l)e + eo + 



(4.21) 



This turns out to be necessary for the most general proof of Eq.(4.3). For the special case 
e„ = e = eo, we may also write 



gp{xi,r\xo,0) = (xi|e-^^^P(ne) • • • P(2e)P(e)|xo) 

Each of the above propagators may be represented by a path integral, 

'1 



(4.22) 



gv{xi,T\xo 



, 0) = J Vx exp j dt 



-mx -\-iVQ6{—x) 



^p(xi,t|xo,0) 



Pxexp i / dt-mx 
P \ Jo 2 



i.2 



(4.23) 
(4.24) 



where in both cases the paths are from x{0) = xq to ,x(r) = Xi and in the second case, 
Eq.(4.24), P denotes that the paths are restricted to be in the positive x-axis at times 
t — 60 + {k — l)e, k = 1,2- ■ - n. 
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A closely related object that will be important is the restricted propagator, 



gr{xi,T\xo,0) = 



xit)>0 



Vxexp 



( 




(4.25) 



where again the paths are from x(0) = xq to x{t) = xi but with x{t) > for all times 
in [0, r]. This is clearly equivalent to gp in Eq.(4.24), in the limit n — )■ oo, e — )■ with r 
constant. If we take the same limit in the equivalent expression for gp, Eq.(4.20), one obtains 
the following convenient operator form of the restricted propagator: 



The restricted propagator satisfies the Schrodingcr equation in a; > subject to the boundary 
conditions that it vanishes when cither end of the propagator sits on x = 0. For the free 
particle, considered here, one can easily solve for the restricted propagator using the method 
of images and the result is 



The restricted propagator, in any of the above forms, describes the regime of "Zeno dynam- 
ics" , in which all states are confined entirely to the Hilbert subspace of states with support 
only in X > [27]. 

The restricted propagator plays an important role here since not only does gp — > gr in the 
limit n — >■ oo, e — >■ with r constant, but also gv ^ gr ^ Vq ^ oo. It follows that gv and 
gp become arbitrarily close to each other for sufficiently large Vq and n, since they both tend 
to the same limit. This is the underlying reason why we expect the approximate equivalence 
Eq.(4.3) should hold. 

The propagators gp and gv niay be decomposed using the path decomposition expansion 
(PDX), which we discussed in detail in Chapter 2. There are two cases to consider. 
We consider first the case in which the initial and final points are on the same side of the 



gr{T,0) = exp{-iPHPT)P 



(4.26) 




(4.27) 
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(xo,0) 



Figure 4.1: The path decomposition expansion for gp. Any path from xq > att — to afinal 

point Xi > at t = T which crosses a; = has a first crossing of x = at ti and 
a last crossing at t2- The crossing part of the propagator from {xq.O) to (xi,r) 
may therefore be decomposed into three parts: (A) restricted propagation entirely 
in X > 0, (B) propagation starting and ending on x — with the restriction that 
the paths are in x > at a discrete set of times, and (C) restricted propagation 
entirely in x > 0. The corresponding path decomposition expansion formula is 
given in Eq.(4.28). (The g^ term corresponds to paths which never cross x — Q). 
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surface in x > 0. Recall from Chapter 2 that the corresponding PDX has the form 

1 r Z"*^ dg dg 
g{xi,r\xo,0) = i^^^ ^ "^^Vo '^^i '^=o 1^=° 

+ gr(xi,r\xo,0) (4.28) 

This is depicted in Figure 4.1. 

Eq.(4.28) holds for both gp and gv, but the restricted part gr is the same in each case and 
equal to Eq.(4.27) above, since the restriction on paths defining gp and the presence of the 
complex potential in gy in Eqs.(4.23), (4.24) are both redundant if x{t) > 0. 

It follows from the above that gp and gv could differ only in terms of their propagation 
along X — 0, hence to prove the approximate equivalence of the propagators Eqs.(4.19), 
(4.20), we need to prove that 

gp{0,t2\0,ti) ^ gv{0,t2\0,ti) (4.29) 

Note that, unlike gv{0,t2\0,ti), gp{0,t2\0,ti) is not covariant under time translation, that 
is, 

5p(0,t2|0,ti) 7^5p(0,t2-ti|0,0) (4.30) 

although it has an approximate covariance on timescales much greater than e, as we will see 
below. 

The second case is that in which the initial and final points are on opposite sides of the 
surface, Xi < and Xq > 0. In this case we have for gv, 

gv{x^,r\xo,0) = ^2 [ dt2 e-^^^^"*^) ^(x^, r|x, t^) L^o 

X gv{0,t2\0,ti) ^(x,ti|a:o,0)|^^„ (4.31) 

since the last section of the paths is in a; < where the complex potential acts. For gp, if 
Xi < 0, since the paths must be in x > at the given discrete set of times, the last crossing 
time t2 cannot be less than the last time r — at which the projectors act, hence 



1 r Z"*^ dg 

gpix,,T\xo,0) = ^J^_^'^^^J^ dt, -^{Xi,T\x,t2)\^^^ 



X gp{Q,t2\Q,h) ^{x,h\xo,0)\^^^ (4.32) 
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Again we will have to prove that Eq.(4.29) holds, but these two expressions Eqs.(4.31), (4.32) 
also differ in the form of the t2 integral. However, since Vqc 1, the form of the exponential 
in Eq.(4.31) effectively squeezes t2 to lie approximately within e of r, so we have 

dt2 e-^°("-*2) fti / dt2 (4.33) 

Jr—en 

and Eqs.(4.31) and (4.32) are approximately the same. 

We see that both cases reduce to proving Eq.(4.29). In Chapter 3 we showed that the 
propagator along the boundary of the complex potential is given by. 

The main purpose of the remainder of this chapter is to calculate the propagator with 
projection operators along the boundary, 

gp{0, t|0, 0) = {0\e-'"'-Pe-'"' ■ ■ ■ e-'"'Pe-'"'°\0) (4.35) 

and show that the approximation Eq.(4.29) holds, under conditions to be determined. (Here, 
|0) denotes a position eigenstate |a;) a,t x — 0). 
In what follows, our main result is to show that, to a good approximation, 

5p(0,i|0,0)« (— j fp{t) (4.36) 

where fp{t) is a kind of saw-tooth function - a piecewisc linear function with peaks imme- 
diately followed by troughs at tk = + {k — l)e. A; = 1, 2, • • • . Approaching tk from below, 
there is a peak of value 

and approaching tk from above there is a trough of half that size. That is 

/p(i) = ^1^^ + for tk-.<t<tk, k^2,3... (4.38) 

and 

fp{t) = 1 for 0<t<eo (4.39) 
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The functions fp{t) and fv{t) are shown in Figure 4.2. Wc see that fp{t) oscillates with 
period e about fvif), as long as we choose Vq so that fvif) lies between the peaks and 
troughs of fp{t). That is, for large k, we require that 



111 

2k^ k 



(4.40) 



Since t ~ ke, fvif) will lie approximately midway between the peaks and troughs of fp{t) if 



RecaUing that the propagator is attached through the PDX Eq.(4.28) to an initial state, 
the oscillations, and hence the differences between gp and gfy, will be smoothed out as long 
as e is chosen to be smaller than the timescale of variation of the initial state. With some 
quahfications (discussed further in Section 4.8), this timescale is the Zeno time, 1/AH. The 
desired approximation Eq.(4.29) will therefore hold in a time-averaged sense, and we thus 
have significant agreement with the extended version of the original argument of Echanobe 
et al described in Section 4.2. 

In the following sections, we evaluate Eq.(4.35) and confirm the form Eq.(4.38) of fp{t). 
In Section 4.4 we evaluate Eq.(4.35) exactly for the cases of one, two and three projections. 
We also show why in general the troughs of fp{t) are exactly half the size of the peaks by 
considering the hmit e„ — > of Eq.(4.35). In Section 4.5 we use a lattice method to derive 
the magnitude of the peaks of fp{t) for large k. These results are substantiated in Section 
4.6, where we use an S-matrix expansion to derive some exact results for a certain time- 
averaged version of Eq.(4.35). In Section 4.7, we fill in some of the gaps in these regimes 
and approximations by computing Eq.(4.35) using numerical methods. 

We will find in the numerical and analytic solutions that the function interpolating between 
the troughs and peaks of fp{t) is not in fact a linear function in general. It is interesting 
however, that, because of the slow time variation of the initial state in comparison to e, 
the precise form of this function turns out to be unimportant, and this is what makes the 
problem more tractable than one might expect. All that is important is the location of the 
peaks and troughs of the saw-tooth function and the period of their oscillation. It is this 
separation of timescales that also restores an approximate time translation covariance to gp, 
even though it does not hold exactly in general, Eq.(4.30). 



4 



(4.41) 
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Figure 4.2: A plot of the functions fv{t) (smooth curve) and fp{t) (saw-tooth curve), defined 
for the complex potential by Eq.(4-34) with V^e = 4/3, and for intermittent pro- 
jectors by Eqs.(4-36), (4-38). (The time scaling is such that e = 1 in the plot, so 
that the peaks occur at integer values oft). We see that fp(t) oscillates around 
fv{t) with period e so the two functions are equal in a time-averaged sense, when 
integrated against functions which are slowly varying on a timescale e. 
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4.4 Exact Analytic Results 



In this section we carry out an exact evaluation of Eq.(4.35) for the case of one, two and 
three projections. For simphcity, we take the initial time interval eo to be e and for the 
case of three projections we are able to carry out the calculation only in when the final time 
€n — €. We thus compute the objects 



gp{0,t\0,0) 



{0\e-'"*\0), iiO<t<e 
(0|e-*^(*-')Pe-^^^|0), if e < t < 2e 



(4.42) 



We will show that 



gp{0,t\0,0) 



i^y ifo<t<6 

H^)' ^ if6<t<26 

I (l + f ArcTan[(t - 2e)/t]^) (^)^ if 26 < t < 3e 



(4.43) 



thereby confirming the approximate form of the saw-tooth function, Eq.(4.38), for the first 
few projections. These expressions all have the property that the value of gp{0,t\0,0) drops 
to half its value immediately after a projection and we show why this is true in general at 
the end of this section. 

To compute Eq.(4.42) it turns out to be useful to consider more general objects, in which 
the projections may be P, P or the identity. We therefore consider the object 



gp,c{0,t\0,0) = J^dxi...dxn[—) i 



m 



2TTi(t — ne) 



1/2 



X exp 



imx^ im , , ^ imxi 



2e 



fc=i 



2(t - ne) 



{AAA) 



Here C — {+,—,0...} symbolically stands for the integration ranges of the x^, eg C = 
{+, — ,0...} means < xi < oo, — oo < X2 < 0, — oo < X3 < oo etc. Performing a Wick 
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rotation, e — > —ie and changing variables yields 

dyi...dyn 



/ \ 1/2 r 

9P,c(0.-.*|0,0) = (-) / 



c(ire)"/2(t-ne)'/2 



fe=i 



27r 



7c(e,i) (4.46) 



The cases of no projection and one projection are trivially evaluated, and we easily obtain 
the first two equations in Eq.(4.43). For the case of two projections, we need to evaluate the 
object 

dyidy2 f yf (yi - 2/2)^ yl 



^2, 63) = / exp -^^ - v,,^^ _ ^ (4 47) 

Jo TT^/ele^ V ei 63 €3/ 

(a slight generalization of the Tc defined in Eq.(4.46)). Changing variables to r = yi, 
y = 2/2/2/1 gives, 

^ . ^ r drdy f 2 f 1 iX-y? 

r++(ei, £2, €3) = / — r exp -r <^ — + h — 

^ 1 

271 (6263 + €163(1 - y)2 + 61622/^) 

(4.48) 



Jo 



Iq 2TTa + by + cy 

Where a — 63(61 + 62), b — —26163, c = 61(62 + 63). Now use the change of variables, 

= ^ + (4.49) 

to obtain 

1 



T^^(e„e2,e^ = ^J— J^^u^^^ (4.50) 



Noting that 4ac — 6^ > 0, and using the standard integral [36] 



/ 



du ^'^ ^ = -ArcTan (-) (4.51) 
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we find finally 



r++(ei,e2,e3) = ^ = ( tt + 2ArcTan ( . / )) (4.52) 

47r V(ei + €2 + eg) V VV ^2(61 + £2 + 63)// 

Setting ei = 62 = e and £3 = i — 2e it follows that 

^p(0,i|0,0) = I (1 + ^ArcTan[(i - 26)/^]^) (^)^ (4.53) 

for 2e < t < 3e, so we confirm the third equation in Eq.(4.43). 

It will be useful for the three projection case below to record the following related results. 
A similar analysis to that above yields 

{yi - y2f yl 



r+_(ei,e2,e3) = / dyi dy2 — , exp — ^- 
^0 ^-00 7r^eie2e3 \ ei 



62 63 



TT - 2ArcTan ( J— — — . 1 1 (4.54) 



47rv/(ei + 62 + 63) V VV ^2(61 + 62 + 63) 

In particular we then have 

T++ = T++(6,6,6) = -^ (4.55) 

3v3e 



We are able to evaluate the three projection case only in the situation where all time 
intervals are equal, so we set i — ne = e in the definition of Tc in Eq.(4.46). In this case, Tc 
possesses a number of helpful symmetries. The first is "reflection" symmetry: if we define 
— C by the string obtained by letting (H — >■—),( >■ +), (0 — >■ 0) then 



-c 



Tc (4.58) 



The second symmetry is "time reversal" : if we define C by the string obtained by reversing 
the order of C, then we have 

Tc = Tc (4.59) 
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In addition we have the simple property that 



T...0... = T..+... + r.._... (4.60) 
These properties imply the following for the three projection case. On the one hand we have, 

T+++ = T++0 - = r++o - To+_ + T_+_ (4.61) 

but we also have, 

= r+o+ - = T+0+ - T_+_ (4.62) 
Combining these expressions gives, 

2T+++ = T+0+ + T++0 - To+_ (4.63) 

For each of the objects on the right, we may carry out the full range integral, to leave 
ourselves with an object of the form ^^(ei, 62, €3) computed above. For example we have 
that ^+0+ = 2++(e, 2e, e). Each of these objects may then be evaluated using Eq.(4.52) and 
Eq.(4.54), and combined to give. 

It follows that J 

when t — 4e, thus confirming the fourth equation in Eq.(4.43). 

To end this section, we confirm our claim in Section 4.3 that the function gp{0, t\0, 0) drops 
to half its value immediately after a projection. On the face of it, this involves interpreting 
the expression (0|P which is ambiguous, so must be defined by a limiting procedure (where, 
recall (0| denotes at x — 0). We thus consider the hmit e„ — >■ in the expression 

gp{0,r\0,0) = {0\e-'"'-Pe-'"'P---e-'"'Pe-'^'°\0) 

= J dy {0\e-'^''^Pe-'"'\y){y\P Pe-'"">\0) (4.66) 

We seek to show that the result obtained by this limit is half the result obtained when the 
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final projection is absent. We write 

P=^-{P + P) + i(P-P) (4.67) 

so that 

(0\e-'Her.Pe-'"'\y) = ho\e-'"'"e-'"'\y) + ^(0|e-'^^"(P - P)e-'"'\y) (4.68) 

The first term on the right-hand side yields the sought after factor of a half as e„ — > 0, so it 
remains to show that the second term is zero. We have 

/"OO 

{0\e-'"'-{P - P)e-'^'\y) ^ dx (0|e-^^^"|x) {{x\e-'^'\y) - {-x\e-'"'\y)) (4.69) 

Jo 

where we have used the fact that 

{0\e-'"'- \x) = {0\e-'"'^ I - x) (4.70) 

(for the free particle, considered here). As e„ — )■ 0, (0|e~*''^'^" |a;) — )■ 6{x) and the pair of terms 
in brackets on the right-hand side cancel in Eq.(4.69), so we obtain zero as required. 

4.5 Asymptotic Limit 

We now consider the asymptotic evaluation of gp{0, r|0, 0) for large n in the case e„ = e = eo. 
This will give the values of the peaks of fp{t), Eq.(4.37). 

Recall that in the limit n — )■ oo and e — ?> with r = (n + l)e fixed, gp{xi, r\xQ, 0) tends to 
the restricted propagator gr{xi,T\xQ,0). To obtain the asymptotic form of (?p(0,r|0,0) for 
large n we therefore need to determine the lowest non-trivial correction to this result. It is 
clear from the definition Eq.(4.20) of gp that, close to the limit, we have the general form 

gp{xi,T\xo, 0) = gr{xi, t\xo, 0) + egi{xi, t\xo, 0) + O(e^) (4.71) 

for some function gi{xi, t\xq, 0). Since the restricted propagator vanishes if either xi = or 
xq — 0, the object we need to calculate is 

9.(0,r|0,0)= lim ?l<^^ (4.72) 

e— >0,n— >oo e 
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A standard and convenient way to do this is to rotate to imaginary time f (we use tilde 
to denote Euclidcanizcd time) and then write gp{xi,f\xo,0) as a Euchdean path integral. 
This is then defined in terms of the continuum limit of probabilities of random walks on 
a space time lattice of temporal spacing Af and spatial spacing 77. The details of this 
construction are very conveniently given by Hartle [59] so we will give only the briefest 
account here. Using this language, gp{0,f\0,0) is then the continuum limit of the object 
{2r))~^up{0,f\0,0), where Mp(0,f|0,0) is the probability for a random walk on the lattice 
starting at the origin and ending at the origin at time f, with the restriction that the walker 
is in the positive x-axis at the intermediate times e, 2e, 3e • • • (where e > Af). 

The calculation of defined in this way, for values of e generally greater than the lattice 
spacing Af, is in fact a known problem in combinatorics called the tennis ball problem. It 
appears to have a formal solution, but this solution is too implicit for us to extract a useful 
result [73]. 

Fortunately, however, for the purposes of calculating the limit Eq.(4.72) the results of 
Rcf. [59] are sufficient. For this case, we set e = Af and Mp(0,f|0,0) is then the probability 
for a random walk from the origin to itself, with the restriction that the walker is in the 
positive X-axis at every intermediate step. On a finite lattice Hartle's calculations give the 
result. 



to leading order for small e, 77 [59]. We may use this result to compute the limit Eq.(4.72), 
which, continued back to real time, is 



Through Eq.(4.71), this confirms Eq.(4.37) for large k. 

4.6 A Time Averaged Result 

We have argued that at the peak values of 5fp(0, t|0, 0) with n projections, we have the 
approximate result 




(4.73) 




(4.74) 




(4.75) 
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Wc showed above that this is exact for n = 1,2,3 and true asymptotically for large n. Some 
interesting exact results for any n may be obtained by considering a more general version 
of this object in which the projections are not restricted to act at the given set of times, 
e, 2e, 3e • • • . 

On the one hand, Eq.(4.34) may be expanded as a power series in powers of Vq: 



9v 



^°-^^'''^=&'"t^Kr' (4.76) 



n=0 ^ ' 



On the other hand, the evolution operator with complex potential may be expanded in the 
usual S-matrix expansion, 

°° / -\\n pT pT 

y^-^ / dhT[Vitn)...V{h)] (4.77) 

n=0 



oo 

-iHt ' 



where T denotes time ordering, V{t) = VoP{t) and P(t) = e'^*Pe-^"\ It follows that 
^y(0,T|0,0) = (0|exp(-i(//-il/)r)|0) 

°° / -\\n pr pr 

= E^^^"/ / rfti (0|e-^^^T[P(t.)...P(ti)]|0) (4.78) 

n=0 



Equating powers of Vq in Eqs.(4.76) and (4.78) and writing out the time ordering explicitly, 
we deduce that 

- 2 dt. I I * (0|e--P(t„)...P(tO|0) = (^) (4.79) 

We would like to write this in a form involving P, instead of P. We introduce the reflection 
operator 

R^j,.\.)(-.\ (4.80) 

and note that P = RPR, the Hamiltonian H commutes with P, and P|0) = |0) (where, 
recall, |0) denotes \x) at x = 0). It follows that 

(0|e-^-P(y...P(ii)|0) = (0|e-^-P(i„)...P(ti)|0) (4.81) 
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so that Eq.(4.79) holds with all the P's replaced with P's. Noting that 

- / dtn dtn-1... / dh^l (4.82) 



T 



^0 ^0 



we see that Eq.(4.79) is of the desired general form, Eq.(4.75), but time- averaged over the 
times of the projections. 

The question is now to what extent the time-averaged expression on the left-hand side of 
Eq.(4.79) is close to Eq.(4.75). We may take this further in two different ways. 

First, note that for a real- valued function / of n variables, we have the mean value theorem 

— / dtn dtn-1... / dtifitn, tn-l • • • il) = fiCn, Cn-l " " " Cl) (4-83) 

"T" Jo Jo Jo 

for some set of numbers ^„ > ^n-i > • • • > Ci the interval [0, r]. The integral in Eq.(4.79) 
is easily made into a real integral over a real-valued function by analytic continuation, and 
we therefore deduce the exact result 

(0|e-'-PK„)...Pte)|0)=(^)"^^ (4.84) 

for some set of projection times > ^n-i > • • • > in the interval [0, r]. 

Second, we may expand the integrand in the left-hand side of Eq.(4.79) about the values 
ke, to get some insight into why these particular values have any special significance. We 
have 



(0|e-^^^P(t„)...P(ti)|0) = (0|e-'^^P(n6)...P(6)|0) 



+ E (^^-^^) ^(0|e-^-P(i„)...P(ii)|0)|^ 

k=l 



+ ••• (4.85) 
Inserting in Eq.(4.79), and noting that 

- / dtn / dtn-i... / dti tk = ke (4.86) 
JO Jo Jo 



(where, recall, r = {n-\-l)e) we see the first order term in the expansion Eq.(4.85) averages to 
zero in Eq.(4.79). This is clearly only true for expansion about the special values tk — ke. We 
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therefore obtain the desired result Eq.(4.75) up to second order corrections. This suggests 
that the values = ke are significant because they give the best approximation to the 
average in Eq.(4.79). These results give evidence that Eq.(4.37) holds approximately for all 
n, including the intermediate values not covered in the previous two sections. 

4.7 Numerical Results 

To support the analytic results described in the previous sections we evaluate Eq.(4.35) 
numerically for up to 20 projections and confirm the conjectured form Eq.(4.36), (4.38). It 
is convenient to define a sequence of functions Fn{t, x) defined by 

Fo(i,x) = (x|e-*^*|0) (4.87) 

for < t < e and 

Fn{t,x) = (a;|e-*-^(*-"^) (Pe"'^^)" |0) (4.88) 

for ne < t < {n+ l)e where n — 1, 2, 3 • • • . The desired object Eq.(4.35) (with, for conve- 
nience, eo = e), is then given by 

gp{0,t\0,0) ^ F^{t,0) (4.89) 
The sequence Fn{t,x) may be calculated using the recursion relation, 

POO 

F^{t,x)= / (a:|e-^^(*-"^)|y) F„_i(ne,y) (4.90) 
Jo 

Using a new time coordinate s defined hj t = se, rotating to imaginary time, and defining 
Fn{s,x) = Fn{—it,x), we have 

for n < s < n + 1. 

We have attempted to find an approximate analytic solution to Eq.(4.91) for large n, 
but without success. However, a numerical solution is straightforward and yields all the 
information we require. This was done using a simple mid-point rule. The lattice size was 
chosen to be of order 10~^, and the results were checked for robustness against changes in 
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lattice size. 

The numerical result for fp(t) (defined in Eq.(4.36) in terms of 5'p(0, t|0, 0)) is plotted 
in Figure 4.3, along with our claimed approximate analytic expression for fp{t), Eq.(4.38). 
We see that there is excellent agreement. The values at the peaks and troughs appear to 
agree perfectly. The only small discrepancy is that the interpolating functions between the 
peaks and troughs are not exactly linear. This discrepancy is only noticeable for intermediate 
values of n and in any event since, as argued, the curve is effectively averaged over time in the 
PDX, this discrepancy is insignificant. We therefore find substantial numerical confirmation 
for our our main result, Eqs.(4.36), (4.38). 



/(') 




Figure 4.3: A plot of the numerical calculation of fp{t) (the set of dotted lines), and the 
conjectured analytic expression for fp{t), Eqs.(4-36), (4-38) (the set of straight 
lines). (Again we use a time scaling such that e = 1). There is excellent agree- 
ment between analytic and numerical results at the peaks and troughs of the func- 
tions. The numerical result shows some departure from exact linear behaviour 
between the peaks and troughs but this is insignificant as argued in the main text. 
A plot of fv{t), the smooth curve, is also given for reference. 

The apparently perfect agreement of numerical results with the approximate analytic 
expression Eq.(4.38) at the peaks and troughs is striking. We wonder if the approximate 
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1.0 r 




Figure 4.4: A plot of the function S(t) defined in Eq.(4-92), describing the oscillations of 
gp{0,t\0,0) around gv{0,t\0,0). The bold line represents the numerical calcula- 
tion and the dashed line the analytic result (from Eqs.(4-34), (4-38)). After the 
first few oscillations, it oscillates between ±1/3 (approximately). 



74 



analytic expression is in fact exact at these points, but we have not been able to prove this, 
except for the cases n — 1,2,3 and for large n. 

A useful way of seeing even more precisely the relationship between fp{t) and fv{t) is to 
define the function 

SW = M|-1 (4.92) 

SO that 

fp{t) = {l + S{t))Mt) (4.93) 

This is plotted in Figure 4.4. It is a simple function oscillating around zero between ±1/3 
with period e. In terms of S{t), the relationship between gp and gv then has the particularly 
simple form 

gp{0, t\0, 0) = (1 + S{t))gv{0, t\0, 0) (4.94) 

This relationship is perhaps the most concise summary of the sought-after connection be- 
tween the propagators. 



4.8 Timescales 

We now give a more detailed explanation as to the timescales involved in proving the ap- 
proximate equivalence of gv and gp. We broadly expect that the appropriate timescale is 
the Zeno time of the initial state, tz = 1/AH. However, we have derived a very precise 
connection between evolution in the presence of a complex potential and evolution with pro- 
jection operators so we are in a position to investigate the specific way in which an initial 
state may discriminate between these two types of evolution. 

We have argued that the equivalence boils down to proving the equivalence of propaga- 
tion along the boundary, Eq.(4.29), and we have shown that gp{0,t\0,0) oscillates around 
gv{0,t\0,0) with period e. Suppose we have an initial state ip{x,0). Let us consider the 
change in the wave function arising from replacing gy with gp using the PDX, Eq.(4.28). It 
is ^ 

6ilj{x,,T) = ^ dhJJ dt, '^{xi,r\x,t2)l^^ Sg{0,t2\0,h) ^(0,^i) (4.95) 

where 

Sg{0, t2\0, h) = gviO, h\0, h) - gp{0, t2\0, h). (4.96) 
The important part of this expression is the ti integral, which we expect will be small if the 
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initial state is sufficiently slowly varying in time. 

The results of Section 4.3 and 4.8 show that Sg{0, t\0, 0) oscillates around zero with period 
e. In particular, Eq.(4.94) shows that 



The explicit form of S{t) is given in Figure 4.4, but its important qualitative feature is its 
oscillation with period e, so for simplicity, S{t) may be loosely modelled by the function 
sin(27rt/e). To be definite we take the initial state to be a Gaussian wave packet, so we 
have 



where E — p^/2m, N is a. normalization factor, and we have ignored wave packet spreading 
effects in evolving the state. For such a state the Zeno time is tz — ma /p. In Eq.(4.95) 
the differentiation of ip produces a prefactor which does not contribute to the leading order 
evaluation of the time integral, and if the range of integration is much greater than e, we 
obtain the order of magnitude result 



If £'e <^ 1, the right-hand side is clearly very small if e <^ i^. li Ee > 2, the right-hand 
side if bounded from above by exp(— i|/e^) so again will be small He <^tz- Hence the Zeno 
time of the initial state is indeed the timescale controlling the validity of the approximation, 
as expected. Note also that Eq.(4.99) goes to zero as e — > 0, as it must, since gp and gv 
become exactly equal (and equal to gr) in this limit. 

The only problematic case is that in which the initial state has energy £■ ~ 1/e. In this 
case, the right-hand side of Eq.(4.99) is not necessarily small and the approximation may 
fail. This is not surprising since it is the case in which the oscillations in time of the incoming 
state are comparable to the temporal spacing of the projections, so that the state can "see" 
the difference between the complex potential and projections at a discrete set of times. 

Hence, apart from the above exception, Eq.(4.29) holds for e <^ tz- As outhned in the 
Introduction to this chapter, reflection from the complex potential is negligible ii Vq <^ E 
which is equivalent to e ^ Therefore there is an interesting regime, namely 



Sg{0,t\0,0) = S{t)gv{0,t\0,0) 



(4.97) 




(4.98) 




(4.99) 



1 



1 



(4.100) 



E 



AH 
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in which the approximate equivalence Eq.(4.29) holds, yet there is negligible reflection. This 
regime is important in, for example, study of the arrival time problem using complex poten- 
tials [49], and we will make use of it in Chapter 5 of this thesis. (See also Ref.[92] for an 
interesting discussion of timescales in the Zeno effect.) 

4.9 Summary and Discussion 

This chapter was physically motivated by the desire to understand the effect of periodically 
acting projections onto the positive x-axis for a free particle, Eq.(4.1). A valuable way to 
proceed is to use the conjectured relationship Eq.(4.3) with a complex potential first put 
forward by Echanobe et al. This connection, together with known results on scattering, 
estabhshes the timescale under which significant refiection occurs in Eq.(4.1). We noted 
that the arguments for the relationship Eq.(4.3) are only heuristic and there is scope for a 
more substantial proof. 

We proved Eq.(4.3) by considering the associated propagators Eqs.(4.19), (4.20). We 
noted that an approximate equivalence between these propagators is expected since both 
propagators tend to the restricted propagator gr{xi,T\xQ,0) in the limits e — >■ 0, n — > oo 
and Vo — > oo. The path decomposition expansion reduced the proof of equivalence of these 
propagators to the simpler case of propagation between points lying on the origin, Eq.(4.29). 
The propagator along the origin for the complex potential gv{0,t\0,0) is already known, 
Eq.(4.34), so the bulk of the proof was to derive the analogous result for the propagator 
with projections, (7p(0, r|0, 0). Our main result was to prove that this propagator has the 
approximate form Eqs.(4.36), (4.38), which we proved using a variety of analytic and numer- 
ical methods. In effect, the main achievement of this chapter has therefore been to obtain a 
good approximate analytic expression for the propagator gp that appears in Eq.(4.1). 

We found that gp{0, t\0, 0) oscillates with period e around ^'^(O, t\0, 0) as long as Voe ~ 4/3, 
a result most concisely summarized in Eq.(4.94). The approximate equivalence Eq.(4.29) of 
these propagators then holds in a time-averaged sense as long as the timescale e between 
projections is much smaller than the Zeno time of the initial state, 1/AH (but may fail in 
the special case when the incoming state is peaked about energy E ~ 1/e). These conditions 
agree in essence with those of an extended version of the results of Echanobe ct al, with an 
advantage over their results that a definite relationship between Vq and e is obtained. We 
noted that their restriction Eq.(4.4) relating Vq and e is in fact stronger than required and 
the equality Vbe ~ 4/3 derived here gives the best approximate equivalence between gp and 
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Let us make some comments here on the possible generahty of the connection Eq.(4.3). 
We first note that it may in fact be written 

^-iHep^-iHe . . . p^-^He ^ {-iHr - VoPt) (4.101) 

where P — 1 — P. We have proved Eq.(4.101) for the case in which the projections are 
onto the positive x-axis, but it seems reasonably clear that the relationship will hold for 
projections onto any region [a, b] of the x-axis (as long as it is not too small) and indeed for 
regions of configuration space in a many-dimensional model. Such a potential has been used 
recently in the decoherent histories analysis of quantum cosmological models [47] . Moreover, 
although the proof of Eq.(4.101) given in this chapter relied heavily on the fact that P 
projects onto position, the form Eq.(4.101) and the heuristic argument for it in Section 4.2 
do not rely on the particular form of P. We therefore conjecture that Eq.(4.101) may hold 
for a wider variety of projection operators, not just projectors onto position. This will be 
pursued elsewhere. 

Finally we note that the left hand side of Eq.(4.101) is a class operator in the decoherent 
histories (DH) approach to quantum theory outlined in Chapter 1. The equivalence expressed 
in Eq.(4.101) means we can also use the right hand side as a class operator. The practical 
value of this is that it opens the way to dealing with time observables in DH by replacing 
strings of projectors, which may be hard to work with, with complex potentials which are 
much simpler. However on a more fundamental level, it also suggest new ways to construct 
class operators which would not have been obvious starting from the exposition in Chapter 
1. We will pursue this in the next chapter. 
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5 Arrival Times, Complex Potentials 
and Decoherent Histories 



Tomorrow, and tomorrow, and tomorrow. 
Creeps in this petty pace from day to day 

Shakespeare 



5.1 Introduction 

We saw in the Introduction, Chapter 1, that there are many different approaches to the 
arrival time problem in quantum theory, (see Fig. (5.1) for a recap of conventions). The 



to — 



t. — 




' ► X 

Figure 5.1: The quantum arrival time problem. We prepare an initial state localized entirely 
in X > and consisting entirely of negative momenta. What is the probability 
that the particle crosses the origin during the time interval [^1,^2]? 
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approach we wish to pursue in the next two chapters is based on the decoherent histories 
approach to quantum theory, as outhned in Sections 1.4 and 1.5. We saw in Section 1.5 
that the key objects for this analysis are the class operators for crossing and not crossing 
the origin during a given time interval. The class operator for not crossing was determined 
to be, 

C„,{t) = Pe-'"'P...e-'"'P (5.1) 

where P — 9{x), there are + 1 projections and r = Ne. In order to avoid the Zeno effect 
we need to leave e finite. However in Section 1.5 we concluded that this object was somewhat 
problematic, firstly because we could not determine the required timescale e, and secondly 
because this object is difficult to work with analytically. 

However, armed now with the results of Chapter 4 we are able to make progress. The re- 
sults of Chapter 4 tell us that evolution with the class operator above, Eq.(5.1), is equivalent 
to evolution in the presence of a complex potential, 

C„c(t) = Pe-'"'P...e-'"'P ^ e-'"^-^°^^, (5.2) 

where P — 9{—x), r = Ne and Vq ~ 1/e. Eq.(5.2) holds provided 

(6.3) 

This means we can compute the probabilities and the decoherence functional by replacing 
class operators of the form Eq.(5.1) with evolution under the complex potential, and this 
simplifies the calculations enormously. The equivalence in Eq.(5.2) also tells us that we must 
have Ee » 1, where E is the energy of the incoming state, if we are to avoid refiection. 

Of course the results of Chapter 4 as they stand only give us the class operator for not 
crossing during a given time interval. We also require the class operator for crossing. We 
shall compute this in this chapter, as well as calculating the probabilities and decoherence 
functional for questions such as, "Did the particle cross during the time interval [ti,t2]?" 
Our end result will be a complete analysis of the arrival time problem in the language of 
decoherent histories. Before we do this, however, it is useful to explore some aspects of the 
arrival time distribution defined using complex potentials. This will help us to understand 
the key features of our final result. 

The first aim of this chapter is to look in some detail at the calculation and properties of 
the arrival time distribution defined using a complex potential. In particular, we will use 
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path integral methods, which in some ways are more concise and transparent than previous 
derivations. The second aim of this chapter is to carry out a decoherent histories analysis of 
the arrival time problem, using the equivalence between strings of projection operators and 
complex potentials demonstrated in Chapter 4. 

Our main result is that for intervals of size A ^ I/Vq, these class operators are given 
approximately by 

= e-*^°^ dt tl) (p s{xt) + 6{xt) p) 
Jtk 

= e-'"°'{e{x{h))-e{x{tk+i))) (5.4) 

Significantly, the dependence on the complex potential has dropped out entirely. We will 
show that there is decoherence for an interesting class of states, and, for such states the 
probabilities are then given by 

q{tk,tk+i) = ['^' dt J{t) (5.5) 

This, we will show, coincides with the expected result Eq.(1.8), when Il{t) is integrated over 
a range of time much greater than 1/Vq. Hence there is complete agreement with standard 
results on the arrival time distribution at sufficiently coarse-grained timescales. Furthermore, 
our results shed some light on the problem of backflow - we find that the situations when 
Eq.(5.5) is negative are those in which there is no decoherence in which case probabilities 
cannot be assigned. 

In the next Chapter, we compute the crossing class operators in a simpler but more 
heuristic way, by exploring a semiclassical approximation to Eq.(5.1) keeping e finite. The 
results are essentially the same. 

The rest of this Chapter is arranged as follows: In Section 5.2 we consider the classical 
arrival time defined via a complex potential. This well help us understand how the proba- 
bilities we obtain depend on the choice of Vo, or equivalently e, in Eq.(5.2). In Section 5.3 
we then repeat this analysis for the quantum case, using the PDX to calculate the arrival 
time distribution function. We make heavy use in this section of the results of Chapter 2. In 
Sections 5.4 and 5.5, we use the results of the previous sections to carry out the decoherent 
histories analysis. We summarize and conclude in Section 5.6. This chapter is based on 
Ref.[49]. 
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5.2 The Classical Arrival Time Problem via a Complex 
Potential 



Before looking at the quantum arrival time problem it is enlightening to look at the corre- 
sponding classical arrival time problem defined using an absorbing potential. This gives some 
understanding of the expected form of the result in the quantum case. In particular since 
the complex potential is a simple measurement model we expect the probabilities obtained 
in this way to be of the form Eq.(1.27), which include an "apparatus response function." 
We will see in this section that exactly the same thing happens in the classical case, and in 
fact the response functions are the same in the quantum and classical cases. This means the 
complex potential is acting like a classical measuring device, and all dependence on it can 
be removed by coarse-graining, in the manner of Eq.(1.28). 

We consider a classical phase space distribution wt{p,q), with initial value wo{p,q) con- 
centrated entirely in q > with only negative momenta. This distribution evolves in the 
presence of a complex potential V{q) — — iVo^(— g). The equation of motion is, 

— ^-^—-2V{q)w. 5.6 
at m oq 

This form may be deduced, for example, by computing the evolution equation of the Wigner 
function with this complex potential and dropping the higher order terms (involving powers 
of %). Eq.(5.6) is readily solved and has solution 

'u;^(p, g) = exp ^— 2 j ds V{q — ps/m)^ Wo{p,q — pr/m) (5.7) 

The survival probability, that is the probability of not having been absorbed by the com- 
plex potential by a time r, is 

/oo roo 
dp / dq Wr{p,q). (5.8) 
oo J —oo 

The arrival time distribution is therefore 



dN f 

n(r) = -—^2Vo dp dq Wr{p, q) (5.9) 

"''^ J —oo J —oo 

where we have made use of Eq.(5.6). This expression is conveniently rewritten by noting 
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that, again using Eq.(5.6), Ii{r) obeys the equation 

^ + 2VoU = -2Vo f dp -Wr{p, 0) (5.10) 
dr J_oo m 



This may be solved to yield 

U{t) = -2Vo[ dte-2^°("-*) / dp -wt{p,0) (5.11) 

From Eq.(5.7), we see that 

Wtip,0) = exp ^— 2Vo J ds 0{ps/m)^ Wo{p, —pt/m) (5-12) 

but since the momenta are all negative the exponential factor makes no contribution. We 
thus obtain 



where 



n(r) = 2Vo r dt e-2^°("-*) J{t) (5.13) 
Jo 

J(t) = - / dp -wo{p, -pt/m) (5.14) 

J —oo ^ 



The current J{t) is the usual classical arrival time distribution that we would have expected 
in the absence of the absorbing potential. 

Eq.(5.13) has the same form as the expressions for the quantum arrival time wc expect to 
obtain from a measurement model, Eq.(1.27). The response function is given by the quantity, 

R{t,s) = 2Voexp{-2Vo{T - s)). (5.15) 

In Chapter 1 we noted that R{t, s) is related to a sort of coarse graining in time, we now 
demonstrate this for this specific case. Eq.(5.13) gives the probability n(r)(ir for arriving 
during the infinitesimal time interval [r, r + dr]. Supposcr we consider the probability for 
arriving during a finite time interval, [ti , T2] ■ This is given by 

T2 



p(T2,ri) = / (ir n(r) 

Jti 

= r dr r dt 2Vo 6-''^°^^-'^ J(t) (5.16) 

Jti Jo 
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Rearranging the order of integration and integrating over r yields, 

p(r2, Ti) = r dt r dr 2H e-2^°("-*) J{t) 

Jo Jti 

+ r dt r dr 2Vo e-2^«("-*) J{t) (5.17) 

Jti Jt 

+ /""^ dt (1 - e-'^°(^^-*)) J(i) (5.18) 

Jti 

1 /Vo plays a role as a fundamental short timescale in the problem, so now suppose we assume 
that Ti, T2 and (t2 — Ti) are all much greater than I/Vq. It follows that all the exponential 
terms may be dropped in Eq.(5.18) and we obtain the very simple result, 

/■T2 

p(t2,Ti) ^ dt J{t) (5.19) 

Jti 

That is, all dependence on the resolution function R{t, s) and the complex potential pa- 
rameter Vo completely drops out when we look at probabilities defined on timescales much 
greater than I/Vq. This result is very relevant to the decoherent histories analysis considered 
later where it is natural to look at the arrival time during a finite time interval much greater 
than e. 

5.3 Calculation of the Arrival Time Distribution via a 
Complex Potential 

We now repeat the analysis of the previous section for the quantum case. The analogue of 
Eq.(5.6) is the inclusion of the complex potential, 

V{x) = -iVoe{-x) (5.20) 

in the Schrodinger equation. With such a potential, the state at time r is 

|V(r)) = exp {-iHoT - Vo9{-x)t) It/j) (5.21) 
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where Hq is the free Hamiltonian. In a similar way to the classical case the idea is that 
the part of the wave packet that reaches the origin during the time interval [0, r] should be 
absorbed, so that 

N{t) = {^{rmr)) (5.22) 

is the probability of not crossing x — during the time interval [0,t]. The probability of 
crossing between r and t + dr is then 

n(r) = (5.23) 

Complex potentials such as Eq.(5.20) were originally considered by AUcock in his seminal 
work on arrival time [6] and have subsequently appeared in detector models [44, 18]. (See 
also Ref.[76, 62] for further work with complex potentials). Prom the results of Chapter 4 we 
see that evolution according to Eq.(5.21) is essentially the same as evolution under pulsed 
measurements. 

For large Vq, the wave function defined by the evolution Eq.(5.21) undergoes significant 
refiection, with total refiection in the "Zeno hmit", Vo — > oo [74]. Here, we are interested in 
the opposite case of small Vq, where there is small refiection and Eq.(5.23) can give reasonable 
expressions for the arrival time distribution. A number of different authors [6, 18, 62] indicate 
that the resulting distribution is of the form 

/oo 
dt R{Vo,T-t)J{t) (5.24) 
■oo 

where J{t) is the current, 

J{t) = {M^-^{pS{x) + S{x)p}\^o) 



and. 



R{Vo, t) = 2Voe{t) exp {-2Vot) (5.26) 



It is therefore closely related to the current J{t) but also includes the infiuence of the complex 
potential via the "apparatus resolution function" R{Vo, t). Our aim in this section is to prove 
this claim. 
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With the complex potential Eq.(5.20), the arrival time distribution Eq.(5.23) is given by 



U{t) = 2Vo{A\0{-x)\i>r) 



(5.27) 



where 



exp {—iHt) IV') 

exp {-iH^T Vod{-x)T) 1^) 



(5.28) 



(so we use H — Hq — iVo9{—x) to denote the total non-Hermitian Hamiltonian) . We are 
interested in calculating this expression for the case in which Vq is much smaller than the 
energy scale of the initial state. (The very different limit, oi Vq — > oo, the Zeno limit, has 
been explored elsewhere [45]). 

One way to evaluate Eq.(5.27) is to use the transmitted wave functions, Eq.(2.17). How- 
ever, we give here instead a different and more direct derivation using the PDX. We use the 
first crossing PDX, written in operator form, Eq.(2.13), 



{x\ eicp{-iHT)\ijj) = / dt {x\ exp(-i//(T - t)) S{x)p exp {-iHot) (5.29) 



Now note that the operator S{x) — |0)(0| has the simple property that for any operator A 



(where, recall, |0) denotes the position eigenstate \x) at a; = 0). Inserting Eq.(5.29) in 
Eq.(5.27), together with the property Eq.(5.30) and the change of variables s — t — t, 
s' — T — t' , yields 



We are aiming to show that this coincides with Eq.(5.24) with the current Eq.(5.25), and the 
main challenge is to show how the p5{x)p combination turns into the combination p5{x) + 
6{x)p in the current Eq.(5.25). 




5{x)AS{x) ^ 5{x){0\A\0) 



(5.30) 




X (0| exp {iHW) \x){x\ exp {-iHs) |0) 

X {ipl exp {iHq^t — s'))p 5{x) p exp {—zHq^t — s)) \ip) 



(5.31) 
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Consider first the x integral. Since we are assuming small Vq, we may use the semiclassical 
approximation Eq.(2.25), which reads 

{x\exp{-iHs)\0) « (^-^y^'exp (^i^ - V^os) (5.32) 
The X integral may now be carried out, with the result, 

n(.) = 4 r^s' [is i^x'^ 



Jr. .In \2TriJ Cs-.sMI 



X {ipr\ exp {—iHQs')p S{x) p exp {iHos) \ipr) (5.33) 

where |'0t) denotes the free evolution of the initial state. 
We now carry out one of the time integrals. Note that, 

pT pT pT pT pT pT 

/ ds' ds^ ds' / ds+ ds ds' (5.34) 

Jo Jo Jo J s' Jo J s 

In the first integral, we set u — s', v — s — s', and in the second integral we set w = s, 
V — s' — s. We thus obtain 

dv — J— 



V2 



X 



— {iPt\ exp {—iHQu)p 5{x) p exp {iHo{u + v)) 



12 



+ — ^^{ipr\ exp {-iHo{u + v))p 8{x) p exp {HIqu) \tpj 

{-i)-2 



(5.35) 



The factors of l/(±i)2 in front of each term come from a careful consideration of the square 
root in the free propagator prefactor (and must have this form because II{t) is real). 

We will assume that VqT ^ 1, which means that the integrals are concentrated around 
u — V — 0. This means that we may take the upper limit of the v integral to be oo, and it 
may be carried out, to yield. 



U{t) = 2Vo f du e 
Jo 



■2Vou 



X ^{^|;r_^\p5{x)E{p) + E^{p)6{x)p\^|Jr-u) (5.36) 
zm 
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where the operator E(p) is given by 



r (5.37) 

[2m{Ho + iVop 



For Vq much less than the energy scale, 



E(p) ^ (5.38) 

so S(p) is simply the sign function of the momentum, which is —1 in this case, since the 
initial state consists entirely of negative momenta. Finally, writing u — t — t, we obtain 

n(r) = 21^0/ dte-^''''^^-'^ ^-^{iJt\{pS{x) + 6{x)p)\ilJt) 
Jo 

= 2^0 Tdt e-2^°M J(i) (5.39) 
Jo 

We therefore have precise confirmation of the classical result Eq.(5.13), and also agreement 
with the expected quantum result, Eq.(5.24). 

Some comments are in order concerning the positivity of the result for II{t). The ex- 
pression (5.36) is positive because it was derived from the manifestly positive expression 
Eq.(5.31). (Two approximations were used: the semiclassical approximation Eq.(5.32), and 
the condition Vqt 1, neither of which affect the positivity of the result.) 

However, to obtain the final result Eq.(5.39) we took the hmit Vq — )■ in the current part 
only of Eq.(5.36), leaving behind the Vo-dcpendent term in the exponential part, and the 
resulting expression is not guaranteed to be positive. In Eq.(5.39), J{t) is not always positive 
due to the backfiow effect [78, 12] and integration over time does not necessarily remedy the 
situation. (See the Appendix for a more thorough discussion of this point). The lack of 
positivity for a n(r) obtained in this way is not surprising since taking the limit Vq — )■ in 
one part of the expression Eq.(5.36) only but not the other will not necessarily preserve its 
property of positivity. The violation of positivity is generally small, however, so Eq.(5.39) 
may still be a good approximation to the manifestly positive expression Eq.(5.36). 

It should also be added that it would be misleading to explore the first order corrections 
in Vo in going from Eq.(5.36) to Eq.(5.39), since comparable correction terms have already 
been dropped in using the semiclassical approximation Eq.(5.32). 
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5.4 Decoherent Histories Analysis for a Single Large 
Time Interval 



We turn now to the decoherent histories analysis of the arrival time problem. The important 
question we wish to answer is; What is the probability of arriving during some general time 
interval [^1,^2]? We answer this question in detail in the next section, here we address a 
simpler question; What are the probabilities for crossing or not crossing during the large 
time interval [0, r]? The reason for tackhng this question first is that it gives us useful 
insight into the more general case and also suggests an interesting way of thinking about the 
situations in which we might have decoherence. 

The class operators for not crossing and crossing during this time interval are 

Cnc = ex.v{-iHoT-Voe{-x)T) (5.40) 

Cc = exp i-iHoT) - exp {-iHoT - Vo9{-x)t) (5.41) 

and they satisfy 

C„c + Ce = e-*^°^ (5.42) 
We are interested in the probabilities for not crossing and crossing, 

Pnc{r) = Tr (CncpCl) (5.43) 

Pe(T) = 'lV(CepCt) (5.44) 

and the off-diagonal term of the decoherence functional, 

= IY(C„ePe^^°-)-p„e (5.45) 

These quantities obey the relation 

Pnc +Pc + ^c,nc + ^^nc = 1 (5-46) 

We look for situations where there is decoherence, 

Dc,nc = (5.47) 



89 



(which is usually only approximate), in which case the probabilities then sum to 1, 



Pc{T)+Pnc{r) = l (5.48) 

It is useful to relate some of these expressions to the standard expressions for arrival time 
n(t) defined in Eqs.(5.22), (5.23). To do this, note that pnc is in fact the same as the survival 
probability, N{t) defined in Eq.(5.22), and that Pnc obeys the trivial identity 

Pnc{r) = 1 + ^ (5-49) 

since Pnc(O) = 1. It follows that 

Pnc{r) = 1 - / dt n(t) (5.50) 
Jo 

When there is decoherence, Eq.(5.48) holds and we may deduce that 

p,(r)= / dtU{t) (5.51) 

Hence the decoherent histories analysis is compatible with the standard result, but only 
when there is decoherence. 

We turn now to the calculation of the decoherence functional. Following the general 
pattern described in Eqs.(1.32), (1.33), consider the quantity, 

qUr) = Tr (C„epe^^°-) (5.52) 

From Eq.(5.45), we see that the decoherence functional may be written, 

Dc,nc = Qncir) - Pnc{r) (5.53) 

This means that g„c = Pnc when there is decoherence. Or to put it the other way round, 
decoherence of histories may be checked by comparing g„c with pnc and this is what we now 
do. Recall that pn^ is given by Eqs.(5.49), (5.50) (which hold in the absence of decoherence). 
We may write Qnc in a similar form: 



^0 



qncir) = 1+1 dt^ (5.54) 
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The integrand is similar to Il{t) defined in Eq.(5.27), so we define 

m - (5.55) 

We now have that 

qncir) = 1 - / dt n(t) (5.56) 
Jo 

It then follows that the decoherence functional is 

^c,nc = j^ dt(u(t)-Il(t)') (5.57) 

To compute the decoherence functional we need to calculate n(i), which is given by 

n{t)^Vo{ip\ exp{iHot) 9{-x) ex.p{-iHot - Vo9{-x)t) (5.58) 

This is almost the same as Hit) except that the exponential on the left involves only Hq 
and not the complex potential (and also an overall factor of 2). We therefore follow the 
calculation of n(t) in Section 5.3 with small modifications. With a little care, one may see 
that the final result is the same as that for Il{t), Eq.(5.39), except that 2Vo is replaced with 
Vq, that is, 

n{t) ^Vo f ds e-^°(^-^^ J{s) (5.59) 

^0 

This result holds for timescales greater than I/Vq and under the semiclassical approxima- 
tion (which required S> Vq so is equivalent to the requirement of negligible reflection 
encountered above). Finally, a calculation similar to that of Eqs.(5.16)-(5.19) implies that 

/ dt n(t) ^ dt J{t) (5.60) 
Jo Jo 

as long as Vqt ^ 1. Since this result is independent of Vq, n(t) will satisfy the same relation. 
We thus deduce that 

Dc,nc ~ (5.61) 

hence there is decoherence, under the above conditions. 

This analysis gives us an interesting way of thinking about decoherence. Decoherence is 
achieved if Eq.(5.59) is equal to Eq.(5.39). This will be the case if the current varies on a 
timescale that is long compared with 1/Vo » 1/E, where E is the energy of the wavepacket. 
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The timescale on which the current varies is given, for a gaussian wavepacket, as 1/AH. In 
Chapter 4 we noted that whereas the timescale given by is genuinely quantum, 1/AH is 
in fact a classical timescale. This means we will have decoherence if there is a large separation 
between the quantum and classical timescales associated with the incoming wavepacket, a 
satisfying result. 

5.5 Decoherent Histories Analysis for an Arbitrary Set 
of Time Intervals 

We now turn to the more complicated question of much more refined histories, that may 
cross the origin at any one of a large number of times, during the time interval [0, r]. This 
corresponds more directly to the standard crossing probability, Il{t)dt, the probability that 
the particle crosses during an infinitesimal time interval [t,t + dt]. 

5.5.1 Class Operators 

We have defined class operators Eqs.(5.40), (5.41) describing crossing or not crossing during 
a time interval [0, r]. We now split this time interval into n equal parts of size e, so r = ne 
and we seek class operators describing crossing or not crossing during any one of the n 
intervals. We first note that 



where C„c(e) and Cc(e) are defined as in Eqs.(5.40), (5.41) except that here they are for a 
time interval [0,e]. We now use this to decompose e^'-^"^ into the desired class operators. 
We have 



e-^^°^ = C„c(e) + Ce(e) 



(5.62) 



e 



iHoT 



^^-iHoey-' (ac(6)+C,(6)) 

(e-^«^)""'c„e(6) + e-^^°(--^)Ce(e) 



(5.63) 



Repeating the same steps on the first term, this yields. 



e-^o- = {e-'"°^y-' C„e(2e) + e-^°(--^^)Ce(e)C„e(e) + e-^°(--)C,(e) 



(5.64) 
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Repeating more times eventually yields, 



n-l 



CUr) + e-^^°(^-('=+^)^)Ce(e)C„,(A;e) (5.65) 



jfc=0 



Prom this expression, we see that the class operator for crossing x = for the first time 
during the time interval [/ce, {k + l)e] is given by the summand of the second term, 

C,{{k + l)e, he) = e-^^°(^-("+^)^)Ce(e)C„e(A;e) (5.66) 

We will not in fact work with the class operator Eq.(5.66), since a more useful similar 
but alternative expression can also be found. Taking the continuum limit of Eq.(5.65) (and 
inserting the explicit expression for C„c), we obtain. 



e 



(where, recall, V = Vq9{—x)). This indicates that the class operator for first crossing during 
the infinitesimal time interval [t,t + dt], is 

Ce(t) = e-^^o("-*)ye-^^°*-^* (5.68) 

We do not, however, expect histories characterized by such precise crossing time to be 
decoherent, so it is natural to consider coarser-grained class operators. 



C^^ / dtCc{t) (5.69) 



which represents crossing during one of the N time intervals [tk, tfc+i] of size A, where tk = kA 
with A; = 0, 1 • • • — 1 and r = A'^A. The complete set of class operators Cq, for crossing 
and not crossing is the set of A'" + 1 operators 

and Eq.(5.67) implies that they satisfy 

iV-l 

e-i^or^Cnc+Y^C^ (5.71) 

fe=0 
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To check for decoherence of histories we need to calculate two types of decoherence functional 



Dkk' = Tr (c,V(^^f )^) (5.72) 
L>,,„e = Tr(C,V(Cnc)0 (5-73) 



and this will be carried out below. 



5.5.2 An Important Simplification of the Class Operator 

There turns out to be a very useful simplification in the class operator Eq.(5.68). Consider 
the amplitude 

{x\e'"°^C^{t)\'il)) = yo(a;|e~'^°*^(-x)e-'-^°*-^*|V') (5.74) 

for any x. The right-hand side is very similar to Eq.(5.27), except that there is no complex 
potential in one of the exponential terms and also the "final" state is \x) not (And there 
is also an overall factor of 2 different). Despite these differences, we may once again make 
use of the details of the calculation of Section 5.3, and we deduce from the analogous result 
Eq.(5.39), that 

(x|e*^°^Ce(t)|V') = Vo ! ds e-^°(*-^) ^^{x\ {p 6{xs) + 6{xs) p) (5.75) 

Jo 2m 

Like the derivation of Eq.(5.39), this is valid under the conditions that all energy scales 
are much greater than Vq and all time scales much greater than I/Vq. Now we integrate 
this over time to obtain the coarse-grained crossing time operator, Eq.(5.69), and again use 
approximations of the form Eqs.(5.16)-(5.19) (again using the assumption of timescales much 
greater than I/Vq), to yield, the remarkably simple and appealing form, 

e''''^^Cl= r^'dt ^ {pS{xt) + 6{xt)p) (5.76) 

This may also be written even more simply 

^iHor^k ^ Q^^^t^^^ _ e{x{tk+,)) (5.77) 
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5.5.3 Probabilities for Crossing 

The above expressions for the crossing time class operator are the most important results of 
this Chapter and provide an immediate connection to the standard expression for the arrival 
time distribution. Supposing for the moment that there is decoherence of histories, we may 
assign probabilities to the histories. The probability for crossing during the time interval 

p{tk,h+,) ^Tr {C',p{C^y) (5.78) 

However, as noted in Eqs.(1.32), (1.33) when there is decoherence of histories, this expression 
for the probabilities for histories is equal to the simpler expression, 

q(tk,h^^) = Tr(C,Ve'^"0 

ftk+l 

^ dt J{t), (5.79) 

which is precisely the standard result! The expression for the probabihty q{tk,tk+i) is not 
positive in general (although is real in this it happens), but when there is decoherence, 

it is equal to p(tfc,tfc+i), which is positive. Hence the decoherent histories result coincides 
with the standard result under the somewhat special conditions of decoherence of histories. 



5.5.4 Decoherence of Histories and the Backflow Problem 

There is an interesting connection between decoherence of histories and backflow. To see this, 
consider the following simple case. We consider histories which either cross or do not cross 
the origin during the time interval [^1,^2]- So the crossing and not crossing class operators 
are C and 1 — C, where 

c = e{xi) - e{x2) (5.80) 

where we have adopted the notation Xk — x(tk) (and for convenience we have dropped the 
exponential factor which is just a matter of definition and drops out of all expression of 
interest). The decoherence functional is 

D = (C(l-C)) 

= {C)-{C') (5.81) 
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This may also be written 



D 



(5.82) 



a form we will use below to check decoherence. When there is decoherence, D — and the 
probability for crossing is 



As noted above, (C) is the standard result, Eq.(5.79), for the probability of crossing. 

There is an interesting connection here between backfiow and decoherence. If there is 
decoherence, D is zero so (C) must cancel (C^) in Eq.(5.81), which means that (C) > 0, so 
there is no backfiow. Or we may make a logically equivalent statement: if there is backfiow, 
(C) < 0, then there cannot be decoherence, since \D\ is then greater than the probability 
(C^). Hence, states with backfiow do not permit decoherence of histories. (Note that absence 
of backfiow, (C) > 0, is not itself enough to guarantee decoherence - the stronger condition 
D = must be satisfied). 

This is an important result. The quantity (C) is regarded as the "standard" result for 
the arrival time probability and its possible negativity is disturbing. Here, the decoherent 
histories approach sheds new light on this issue. In the decoherence histories apporach, the 
true probability for crossing is the manifestly positive quantity (C^) and this is equal to 
(C) only when there is decoherence. In particular, when there is significant backfiow, there 
cannot be decoherence, so probabilities cannot be assigned and (C^) is not equal to (C). 

5.5.5 The Decoherence Conditions 

The crossing probabihties described above are only vahd when all components of the deco- 
herence functional, Eqs.(5.72), (5.73), are zero. We therefore address the issue of finding 
those states for which these conditions hold. 

The most important decoherence condition is that Dkk' defined in Eq.(5.72) vanishes, so 
we focus on that. We write the class operator Eq.(5.77) for crossing during the time interval 
[tk,tk+i] as 



For an initial state \ip), the quantity C^lip) is a quantum state representing the property 
of crossing of the origin in the time interval [tk-tk+i]- The decoherence condition Dj^k' = 
is simply the condition that the "crossing states" C^\ip) for different time intervals have 



p{t,,t2) = (C) = (C) 



(5.83) 




e-^°- {9{xk) - 9{xk+i)) 



(5.84) 
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negligible interference. The states C^\ip) consist of an initial state which has been localized 
to a range of time at a; = 0. This is closely related to the interesting question of diffraction 
in time [75] and this connection will be explored in more detail elsewhere [53] . 
The decoherence functional is given by 

= mxk) - e{xk+i)) - e{-xj))) (5.85) 

where without loss of generality we take tj^i < tk- It is a sum of terms each of the form, 

dkj = {e{-Xk)9{x,)) (5.86) 

where tk < tj. Note that 

|4,f <(t^ {e{-xk)e{xj)e{-xk)) (5.87) 

The key thing is that rf^^ has the form of a probability - it is the probability to find the 
particle in a; < at and then in a; > at tj. Semiclassical expectations suggest that this 
is small in general for the states considered here, which are left-moving wave packets, and 
indeed the classical limit of this probability is zero. So this is a useful object to calculate 
in terms of checking decoherence. (Although it may not be small for states with backfiow). 
Note that it also implies that the other parts of the decoherence functional, Eq.(5.85) will 
also be small. In detailed calculations, the upper bound Eq.(5.87) must be compared with 
the probabilities, as in Eq.(1.40). 

5.5.6 Checking the Decoherence Condition for Wavepackets 

Decoherence starts to become lost as the the size A of the time intervals [tk, tk+i] is reduced to 
close to 1/E. This is because the "locahsation in time" produces a corresponding spreading 
in energy, and thus momentum, and for sufficiently small time intervals this is enough to 
give rise to positive momenta with significant probability. We will see in a specific example 
how this works. 

For simplicity we work with the simple wavepacket initial states, more general initial states 
require an environment to produce decoherence and will be considered in Chapter 7. We 
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take the initial state to be, 



where go > 0, po < and |po|c" >> 1- We note that the decoherence functional Eq.(5.82) 
satisfies 

\D\' < 2dl (5.89) 

with given by Eq.(5.87) with k — 1, j — 2. We need some probabilities to compare this 
with. We have that 

|L>|' < (C')((l-C)') (5.90) 

The interesting case is that in which the crossing probability (C^) is somewhat less than 1, 
less than about 1/2, say, in which case the non-crossing probability ((1 — C)^) will be of 
order 1. It is therefore sufficient to compare d"^ with (C^). Now note that 

= {e{x^)e{-x2)e{x^)) + {e{x^)e{-x2)e{-x^)) + {e{x2)e{-xi)). (5.91) 

This means that (C^) is in fact equal to the probability, 

P12 = {e{x{)d{-X2)e{x^)), (5.92) 

up to terms which vanish when D — Q. This is useful since it is now identical in form to the 
expression for and our goal is to show that 

(t<^pi2. (5.93) 
It is useful to work in the Wigner representation [10], defined, for a state p{x,y) by 

W{p, <D = l^jdi e-^^« p{q +\i,q- (5.94) 
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The probabilities pi2 and df^ are then given by 

Pi2 = 2tt J dpdqWp{p,q)Wo{p,q,ti), (5.95) 

dl = 271 J dpdqWD{p,q)Wo{p,q,h). (5.96) 

Here, Wq{p, q, ti) is the Wigner function of the initial state, evolved in time to ti, 

WM exp (- ("-"o-M' _ 2.^(p - p„)=) . (6.97) 

The objects Wp and Wd are the Wigner transforms of the ^-function combinations appearing 
in Eqs.(5.87), (5.92) and are given by 

Wp{p,q) = -^^(g) r dy'-^, (5.98) 



where 



^(p,g)=2g(p+^^-^). (5.100) 



We see that the only difference between the expressions for pi2 and o?^ is in the sign in the 
^-functions. 
The integral 

f(u)^£dy'^ (5.101) 
may be expressed in terms of the Sine integral function Si{x), 

f{u) = ^ - Si{u) (5.102) 

but its properties are not hard to see directly. For large negative u, f{u) ^ n, at u — 0, 
/(O) = 7r/2, and for large positive u, f{u) goes to zero, oscillating around 1/u. (See Figure 
5.2.) 

We now compare the size of pi2 and d'^.MVe assume that the wave packet is spatially broad, 
so a is large and the Wigner function Eq.(5.97) is therefore concentrated strongly about 
P = Po < 0. We therefore integrate out p and set p — po throughout. The most important 
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Figure 5.2: The function f{u). It oscillates around zero for ^ and oscillates around tt 
for M ^ 0. As a function of q, f{u(po, q)) differs from or tt only in a region of 
size l/|po| around q — 0. 

case to check is that in which the wave packet is reasonably evenly divided between x > 
and a; < at time ti so that both pi2 and df^ have a chance of being reasonably large. This 
means that qq + poti/m should be close to zero (to within a few widths a), so for simplicity 
we take it to be exactly zero. 

With these simplifications, we have 

Here, since q < and Pq < 0, we have u{pq, q) > 0. We can evaluate this expression by 
examining the comparative effects of f{u) and the exponential term. From the plot of f{u) 
(see Figure 5.2), we see that it drops to zero a.t u = Uq (which is of order 1) and oscillates 
rapidly around zero for m > i^o, so we expect the integral to be dominated by values of q for 
which < u < uq. The value u — uq corresponds to q — qo where 



|po|A 



1+ "° 



EoA 



1/2 

-11, (5.104) 



100 



where Eq — p'^/2m. In the complex potential calculations, we have assumed that Vq 
and we also assumed that all timescales are much greater than I/Vq, and these together 
imply that EqIs. ^ 1. We may therefore expand the square root to leading order and obtain 

Wc have assumed that the wave packet is sufficiently broad that apo 3> 1, and this means 
that 

\qo\ < a. (5.106) 

This in turn means that f{u) is significantly different from zero only in the range \q\ <^ a, 
and most importantly, in this range, the exponential term in Eq.(5.103) is approximately 
constant. We may therefore evaluate Eq.(5.103) by ignoring the exponential term, integrating 
from to qo and approximating f{u) as 

f(u) ^^-u + 0{u^). (5.107) 

We thus obtain the simple result. 

In the expression for pi2, there is a key difference in that q > which means that u{po, q) 
can be positive or negative. Introducing 

q, = = (5.109) 

(where is the Zeno time) we see that u < ior q < q^ and w > for g > g^- We therefore 
have 

= (2vr3a^)i/2 / ^^^^P(-^) /W^^o,g)) 

^ {2n^l^y/^ l dq exp(^-^^ fiuipo,q)). (5.110) 

Here, f{u{po, q)) ~ tt in the first term, differing from this value only in a region of size l/po 
close to q — 0. In the second term f{u) will tend to be small except for a small region of 
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size 1/po around the origin. 

If (iz ^ o", (that is, /S.^ tz), the second term in pi2 is exponentially suppressed and in the 
first term the integration range is effectively to oo, so we obtain 

1 

Pi2 ~ 2" (5.111) 

This is the expected result, since under the above assumptions on the wave packet, half of 
it will cross a; = if the time interval is sufficiently large. Clearly 3> o?^ in this case so 
there is decoherence. 

If Qz <^ cr, the first term in pi2 is of order qz/<y and the second of order l/(|po|c"), the same 
order of magnitude as d"^. Hence in this case we have decoherence if 



which is equivalent to, 

£;A>1. (5.113) 

Crucially this is equivalent to the condition that there be negligible reflection from the 
complex potential. 

In brief, we therefore get decoherence of histories for a single wave packet under a wide 
variety of circumstances. It is easily seen that similar results also hold for orthogonal su- 
perpositions of gaussian wavepackets, see Fig. (5.3). More general states will require an 
environment to ensure decoherence. This will be addressed in Chapter 7. 

The condition Eq. (5.113) has previously been obtained in studies of the accuracy with 
which a quantum system can be used as a clock [88, 4] (see also Chapter 8.) In these studies 
it was proposed that 1/ E is a fundamental limitation of the accuracy with which times can 
be defined in quantum theory. One weakness of those studies, however, is that it is always 
possible that the limitation is just an artifact of the particular choice of clock model. It is 
therefore very interesting that the same condition arises as a decoherence condition in the 
study of the arrival time problem presented here. The fact that what is, on the face of it, 
a very different analysis of the problem nevertheless arrives at the same limitation on the 
accuracy of the definition of arrival times adds weight to the idea that Eq. (5.113) really is a 
fundamental limit in quantum theory. 
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Figure 5.3: An initial state consisting of a superposition of wave packets may have signifi- 
cant crossings in at least two different time intervals. If the initial packets are 
orthogonal, and the time intervals are sufficiently large (greater than 1/E), the 
packets will remain orthogonal after passing through the time intervals and the 
corresponding histories will be decoherent. 

5.6 Summary and Conclusions 

We considered the arrival time problem using a complex potential to kill paths entering 
X < 0. In Section 5.2 we gave a classical analysis of the problem. We derived a result of the 
expected form exposing the resolution function as an essentially classical effect summarizing 
the role of the complex potential. We also showed that coarse graining over time scales 
much greater than 1 /Vq produces a formula for the arrival time of expected form and which 
is independent of the complex potential. This is an important result for the rest of the paper. 

In Section 5.3, we used the PDX to rederive the standard form of the arrival time distri- 
bution with a complex potential, in the limit of weak potential. The form of this calculation 
turned out to be useful for the subsequent work on the decoherent histories approach. 

In Section 5.4, we considered the decoherent histories analysis for the simple case of a 
particle crossing or not crossing x — during a large time interval [0, r]. We found the 
simple and expected result that the histories are decoherent as long as reflection by the 
complex potential is negligible. The resultant probabilities are consistent with the standard 
result for the arrival time. 

The main part of the decoherent histories analysis was given in Section 5.5, where we 
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first derived the class operators described crossing a; = for an arbitrary set of small time 
intervals. Here we obtained our most important result: the crossing class operator Eq.(5.76) 
for timescales much greater than 1/Vq. This form of the class operator gives an immediate 
connection with the standard result for probabilities when there is decoherence. Indeed, one 
may have guessed the form of the class operator from the standard form of the probabilities, 
and this is pursued in the next chapter. However, it is also gratifying that it can be derived 
in some detail using the complex potential approach used here. 

To assign probabilities, the decoherence functional must be diagonal and we considered 
this condition. We found a variety of states for which there is decoherence, under certain 
more detailed conditions, which we discussed. 

We also noted an interesting and important relationship between decoherence and back- 
flow: If there is decoherence, the probabilities for crossing must be positive so there cannot 
be any backfiow. If there is no decoherence, the integrated current may still be positive, but 
one can say that if there is backfiow there will definitely be no decoherence. This means that 
the decoherent histories approach brings something genuinely new to the arrival time prob- 
lem: it establishes the conditions under which probabilities can be assigned and in particular 
forbids the assignment of probabilities in cases where there is backfiow. 
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6 Arrival Times, Crossing Times and 
a Semiclassical Approximation 



Some historians hold that history ... is 
just one damned thing after another. 



Arnold Toynbee 



6.1 Introduction 

We have seen in the previous chapter that arrival time probabilities can indeed be defined via 
decoherent histories, with the class operators defined in terms of complex potentials. However 
probably the single most important result was that the class operators corresponding to 
arriving in a time interval [ti, ^2] can be written as, 



with P{t) = 0{xt), provided E{t2 — ti) » 1 where E is the energy around which the 
state is strongly peaked. In this limit, all dependence on the complex potential has dropped 
out and the class operator is exceedingly simple. This strongly suggests that there is some 
simple approximation that can lead us to these class operators, without having to resort to 
complex potentials. The calculations in the previous chapter would then be important only 
to establish the conditions under which this approximation holds. One aim of this chapter 
is to present such a derivation. 

Eq.(6.1) was derived for the special case of a wavcpackct incoming from the right, with 
entirely negative momentum. This analysis can be changed in a simple way to allow for 
states consisting entirely of positive momenta, but it is less clear how to handle states with 
momenta of both signs. 




(6.1) 
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This is important for the following reason. Consider the case where the initial state is 
a superposition of two wave packets, identical except that one is the mirror image of the 
other, ie a; — )■ —x,p — )■ —p. The symmetry here implies that the current at the origin is 
always zero, even though we would like to say that both wave packets cross the origin. This 
suggests that crossing probabilities cannot in general be simply related to the current, and 
it is of interest to know how they may be defined instead. 



t 







t2 - 









Figure 6.1: Crossing times. Consider an initial state which is a superposition of a left and a 
right moving wavepacket. If we can define the probability of crossing for the left 
and right moving wavepackets separately, it seems reasonable to expect that we 
can define the crossing time probability for the superposition. 

Useful insight can be gained by looking at the classical case. There if one has momenta 
of both signs one can seperate the classical trajectories into ones that cross the origin from 
the left and ones that cross from the right [77, 19]. The crossing time distribution, Il{t) can 
therefore be written as 

n(t) = n+(t) + n_(t) = j+(o,t)- j_(o,t), (6.2) 

where the current associated with the positive momentum components of the system is given 
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by 

J+(0,t)= / dp^wt{p,0) (6.3) 
Jo m 

and similarly for J_(0,t). Here Wt{p,q) is the classical phase space distribution function. 
This can be written in the alternative form 

U{t) = r dp^-^w,{p, 0) = r dp^wtip, 0) (6.4) 

7-0O J-oo "I 

where e{p) is the signum function. 

In the quantum case, however, we do not have a description in terms of trajectories. This 
means that we cannot justify the procedure used in the classical case. In fact it turns out 
that even if we could, we would not arrive at an acceptable arrival time distribution. This 
is because the quantum analogue of Eq.(6.4) would be 

n,^(t)= / dp^-^Wt{p,0) (6.5) 

J ITT' 

where Wt{p, 0) is the Wigner function of the system. However unhke the classical phase 
space distribution function, the Wigner function need not be positive. Therefore in general 
this expression will not be positive, so cannot be a probability distribution. It is important 
to note that this negativity is not due to any "trajectories moving backwards" since the 
modulus of the momentum has been taken. Rather this negativity is the result of quantum 
interference between wave packets moving in different directions. The second aim of this 
chapter therefore, is to find the correct class operators to describe arrival times where the 
initial state is a superposition of left and right moving wavepackets, to compute the crossing 
probabilities thus obtained, and to check the conditions under which we have decoherence. 
Our aim will be to demonstrate that if we have a left and a right moving wavepacket for 
which the arrival time histories would individually decohere according to the analysis of 
Chapter 5, then the crossing time histories of a superposition of such wavepackets will also 
be decoherent. 

Our analysis follows from the simple observation that although in standard quantum 
theory there is no notion of trajectory, in decoherent histories the concept of a history is 
very similar to that of a trajectory. We will use this insight to show how an analysis similiar 
to the classical case can lead us to an expression for the object representing the probability 
of crossing a surface from either side. 
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The rest of this chapter is arranged as follows. In Section 6.2 we present a semiclassical 
derivation of the crossing class operators previously derived via complex potentials in Chap- 
ter 5. In Section 6.3 we address some issues relating to the differences between arrival times 
and crossing times that are relevant for the rest of this chapter. In Section 6.4 we introduce 
the crossing time class operators, and then simplify them using a similar semiclassical aprox- 
imation to that used in Section 6.2. In Section 6.5 we compute the crossing probabilities 
and decoherence conditions using these crossing class operators. We conclude in Section 6.6. 
The first part of this chapter, Section 6.2, is based on Ref.[50]. 



6.2 Semiclassical Derivation of the Arrival Time Class 
Operators 

In this section we give a semiclassical derivation of Eq.(6.1) following Ref.[50]. We wish to 
see if there is a simpler way of arriving at the main results of Chapter 5, without recourse 
to the machinery of complex potentials. Instead we see if some semiclassical approximation 
can lead to the same conclusions. 

We consider an initial state at t — with p < concentrated in a; > 0. Suppose the state 
crosses the origin during a large time interval [0,t]. It is useful to introduce a discrete set 
of times tk — ke, where k — 0,1 ■ • • n and r = ne. We also define the projection operators 

P = e(x) (6.6) 
P = l-P^e{-x) (6.7) 

Consider first the class operator for remaining in x > (i.e., not crossing x — 0) during 
the time interval [0,t]. We assert that the appropriate class operator is 

Cnc^P{tn)---P{t2)P{tl) (6.8) 

This corresponds to the statement that the particle is in a; > at the discrete set of times 
^1, ^2, ■ ■ ■ tn, but its location is unspecified at intermediate times. We know from our studies 
of the Zeno effect that we must have e » 1/ E to avoid reflection. 

We construct the first crossing class operator by partitioning the histories according to 
whether they are in a; < or a; > at the discrete set of times tk and noting that the class 
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operators must sum to the identity. We write 



1 = 



P{ti) + P{h) 



(6.9) 



Repeating inductively, we obtain 



n 



1 = P{h) + Y,Pitk)P{tk-i)---P{h)P{h) 



k=2 



+ P{tn)---P{h)P{tl). 



(6.10) 



We thus identify the first crossing class operator as 



Ck^P{tk)P{h-l)---P{t2)P{tl) 



(6.11) 



for k > 2, with Ci = P{ti). This clearly describes histories which arc in x > at times 
ti,t2, ■ ■ -tk-i and in ,T < at time t^, so, to within the accuracy set by e, describe a first 
crossing between and t^. The last term in Eq.(6.10) is the non-crossing class operator, 
Cnc- We will actually assume that r is sufficiently large that the wave packet ends up entirely 
in a; < at large times, so Cnd'ip) is essentially zero. This means that we effectively have 



as required. We will generally be interested in class operators describing crossing in intervals 
[ta, ta+i] of size A = me, where m is a positive integer, and these class operators are simply 
obtained by summing. 



Last crossing class operators are similarly constructed but will not be required, as we shall 
see below. 

Consider the strings of identical projection operators P at different times appearing in 
Eq.(6.8) (and similarly in Eq.(6.11)). Given that the final projection P{tn) is onto a; > 
and also that the initial state is localized in ,t > 0, it seems reasonable to suppose that 
the projections at times ti to do not disturb the evolving state too much, under the 



n 




(6.12) 



k=l 




(6.13) 
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condition e > 1/E discussed above. It therefore seems reasonable to make the approximation 

P{tn)---P{t2)P{hm^P{tnm (6.14) 

This is easy to understand in a path integral representation. The right-hand side is in essence 
the amplitude from an initial state concentrated in a: > to a final point in a; > at time tn 
(up to overall unitary factors). The sum over paths will be dominated by the straight line 
path, which lies entirely in x > at all intermediate times. It will therefore be httle affected 
by the insertion of additional projections onto a; > at intermediate time. 

Using this approximation, the crossing class operator Eq.(6.11) operating on the given 
initial state may be approximated as 

Ck « P{tk+i)P{tk) (6.15) 
Rearranging and using the approximation Eq.(6.14) a second time we obtain 

Ck = P{tk) - P{tk+i)P{tk) 

« P{tk) - P{tk+i) (6.16) 

Which is precisely what we set out to show. We have therefore shown that a simple semi- 
classical approximation, Eq.(6.14), can yield the correct class operators without recourse to 
complex potentials. However, it is not clear from the analysis above when this semiclassical 
approximation is valid. For this we do need the results of the previous chapter, which tell us 
that we must have t2 — ti » 1/E. We also note that the semiclassical approximation used 
above means that there is no distinction between first and last crossing, so a last crossing 
class operator would yield the same result, see the remarks at the end of Section 2.4 

6.3 Arrival Times, Crossing Times and First Crossing 
Times 

Given the results of the previous section, it is straightforward to write down the class op- 
erators corresponding to crossing or not crossing from either side of the origin. Before wc 
do this however we wish to pause a moment to address a subtle issue in the arrival time 
problem, which is the difference between crossing times and first crossing times. In classical 
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mechanics, for a particle in a potential, a particle may cross the surface a; = many times. 
One generally defines the arrival time as the first crossing of a; = 0. In fact constructing 
even the classical arrival time distribution for a particle in a confining potential is somewhat 
involved. Generally one imposes absorbing boundary conditions on the system, and the ar- 
rival time distribution can then be taken as the flux at x = for the modified system with 
these boundary conditions. For the case of a free particle, of course, the particle only crosses 
X = once, and so there is no distinction between crossing and first crossing. 

In the quantum case, a similar problem arises for a particle in a potential and presumably 
we could also tackle this problem using the complex potential analysis used in Chapter 5, 
at least for suitably smooth potentials. However even for a free particle in quantum theory, 
the set of possible paths summed over in the path integral definition of the propagator is 
dominated by those that cross a; = an infinite number of times (see Chapter 3.) This 
means that we must distinguish between notions of crossing and first crossing even for a free 
particle in quantum theory. 

In practice, however, this issue is unlikely to cause real problems. Since we know from 
the Zeno effect that we can only deflne arrival times up to some accuracy e >> 1/E, it 
seems reasonable to assume that we will not resolve multiple crossings of the origin by a free 
particle^. 

The consequence of this is that although we presented the derivation of the class operators 
above as if we used a single semiclassical approximation, in fact there are two different 
approximations being used. The first is the approximation, 

C\ = P{tk)P{tk-i) ■ ■ ■ P{t2)P{ti) ~ P{tk)P{tk-i) (6.17) 

As claimed, this is based on the idea that the propagator is dominated by the straight line 
path. Alternatively this may be thought of as expressing the equivalence of first crossing 
and crossing times. In fact this object is an exact class operator, but for a slightly different 
question, which is "What is the probability that the particle crossed the origin from the right 
during the interval [tk-i,tk\'?" See Fig. (6.2) for an illustration of this. 



This assumes that the timescale on which these effects are important is given by 1/E. Since for a 
wavepacket this is the genuine quantum timescale this seems reasonable, but this will be studied further 
elsewhere. 



Ill 



t t 




(a) First Crossing Times (b) Crossing Times 

ure 6.2: (a) First crossing, or arrival time, class operators are sums over paths which 
first cross x = in the interval [tfc_i,tfc]. The extra projections at earlier times 
are included to ensure paths have not previously crossed x = 0. (b) /n contrast 
crossing time class operators are sums over paths that cross x — 0, not necessarily 
for the first time, in the interval [tk-i,tk\. 
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The second approximation we made was the one leading to Eq.(6.29), 

Ck = P{tk+l)P{tk) 

^ P{tk)-P{tk+i). (6.18) 

This is not so much a semiclassical approximation (it is not true in the general classical 
case), rather it depended on us having an entirely a left moving wavepacket. In fact the 
arrival time class operator Eq.(6.18) can be written exactly as, 

P(tk) - P{tk+i) = P(tfc+i)P(tfe) - P{tk+i)P{h), (6.19) 

so that the approximation Eq.(6.18) corresponds to ignoring crossings from the left. 

These remarks will be important when we come to analyse the class operator for crossing 
from either side of the origin. 



6.4 Crossing Class Operators in Decoherent Histories 

We turn now to the questions of defining class operators for crossing the origin when we 
have momenta of both signs. The class operators for crossing the origin from the right, Ccr 
and from the left, Cd are clearly given by, 

Ccr{tk+l,tk) = P{tk+l)P{tk), (6.20) 

Cci{tk+utk) = P{h+i)P{h). (6.21) 

These expressions are simply Eq.(6.15) and its obvious analogue for the case of a right moving 
wavepacket. Our provisional class operators for crossing and not crossing from either side 
during the interval are thus, 

Cc{tk+l,tk) — Ccr{tk+l,tk) + Ccl{tk+l,tk) 

= P{tk+i)P{tk) + P{tk+i)P{tk), (6.22) 

Cnc{tk+l,tk) — i — Cc{tk+l,tk), (6.23) 

Note that, in light of the comments made above, we cannot simplify Eqs.(6.20) and (6.21) 
in the manner of Eq.(6.18). 

In principle Eqs.(6.22) and (6.23) provide all the information we need to define crossing 
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probabilities in quantum theory. In practice however these objects are somewhat awkward 
to work with. What we want is to use the same sort of semiclassical approximations we used 
to simphfy the arrival time class operators to write these crossing class operators in a more 
useful form. 

We begin by defining new class operators which are simply the real part of Eqs.(6.22) and 
(6.23), 

Cc{h,h) = ^{Cc{h,t2) + Cl{h,t2)) (6.24) 
Cnc{tl,t2) = 1-Cc{h,t2). (6.25) 

The use of the Hermitian class operators, Eqs.(6.24, 6.25) is not strictly necessary, since 
decoherence will ensure that, when defined, the probabilities computed from the class opera- 
tors Eqs.(6.22), (6.23) are real. However using the Hermitian class operators makes it easier 
to compare the probabilities with the standard results. 

Now consider the class operator for crossing from the right during the interval [ti,t2]i 
Eq.(6.20). Crossing from the right is, semiclassically, the same as restricting to negative 
momentum and what we want to show is that Eq.(6.20) is equivalent to some simpler class 
operator which includes a restriction on the momentum. We begin by noting that Eq.(6.20) 
can be written as, 

Ccr{ti,t2) = P{t2)P{h)e{p) + p{t2)P{h)e{-p). (6.26) 

We claim that, semiclassically, 

. p{t2)P{h)e{p) ^ 0, 

. p{t2)P{h)e{-p) ^ {P{t,)-p{t2))9{-p). 

The second approximation follows from Section 6.2. The first one is clearly true classically, 
and will be approximately true in the quantum case provided t2 — ti » ^/E where E is the 
energy of the incoming state. These two approximations mean we can write, 

Ccr^{P{t,)-P{t2M-p). (6.27) 

Using the same approximations in Cd, we can write down a semiclassical approximation to 

Cc, 

CcitiM) = (^(^i) - P{t2)m-p) - eip)). (6.28) 
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Finally then we take the real part of this expression to obtain, after some rearranging, 




(6.29) 
(6.30) 



Eqs.(6.29) and (6.30) are the desired class operators for crossing and not crossing the origin 
from either side during the time interval [tk^tk+i\- 

As an aside, it is interesting to note that that Eq.(6.22) may also be obtained from 
Eq.(6.19) by the following argument: Start with the class operator corresponding to the 
integrated current, Eq.(6.19). Notice that the second term on the right of Eq.(6.19) is asso- 
ciated with crossing the origin in the opposite direction to the first term. Flip the sign of 
this second term to correct for this. This leads us to Eq.(6.22). This is exactly the same 
argument as leads to Eq.(6.2) in the classical case. 

Let us state this argument again in a different way. The object, P{ti) — P{t2) is the 
operator representing the integrated current in standard quantum theory, and is thus the 
candidate for the arrival or crossing time distribution. Eq.(6.19) is an identity, but neither 
of the objects on the right hand side has an obvious interpretation in standard quantum 
theory (one would have to talk about measurements to give them meaning.) However in 
decoherent histories we can easily identify these objects as class operators for crossing from 
the right and from the left. We can therefore perform the classical trick of changing the sign 
in front of the second term, to arrive at Eq.(6.22), which is the crossing time class operator. 
Eqs.(6.22) and (6.23) may therefore be derived either from first principles, or by splitting the 
object representing the integrated current, Eq.(6.18) into two objects representing crossing 
from the left and from the right and then making the change of sign in analogy with the 
classical case, Eq.(6.2). 

6.5 Crossing Time Probabilities and Decoherence of 
Crossing Histories 

Now we have the class operators for crossing, we can ask about the probabilities for crossing 
the origin. The easiest way to approach this is to express these probabilities in terms of the 
momentum representation of the density matrix. For comparison we note that the current 
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can be written in this way as 



Jit) 



Tr{P{t)p) = - / dpdp' 



, [P + P') 



p{p,p')e 



(6.31) 



2m 



For a decoherent set of histories, the probabihties are given by, 



P{tl,t2) 




(6.32) 



This means that the probabihties are in fact functions not of the fuU density matrix, but 
rather of the quantity. 



This is somewhat unexpected, but very appeahng. For a superposition of left and right 
moving states, our arrival time probabilities are just the sum of the probability for the arrival 
of the left moving state and that for the right moving one. Interference terms between left 
and right moving states do not appear. Note the close similarity between the way our class 
operators partition the current into left and right moving parts and the classical expression. 



We turn now to the question of for which states and intervals we have decoherence. De- 
coherence between crossing and not crossing during an interval [^1,^2] requires. 



P = o{p)pO{p) + o{-p)pe{-p). 



(6.33) 



Eq.(6.2). 




(6.34) 
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which can also be written as 

D{c, nc) = ^ ' dm ({e{-p)P{t)e{-p) - e{p)P{t)e{p)}p^ 

- J ' dtdsTv {{e{-p)P{t)e{-p) - e{p)P{t)e{p)}p 

X {e{-p)P{s)e{-p) - e{p)P{s)e{p)}) , (6.35) 

D{c, nc) - r dtTr (e{-p)P{t)e{-p)p) - j ' dtdsTi (e{-p)P{t)e{-p)pe{-p)P{s)e{-p)] 

dm {e{p)P{t)e{p)p^ + dtdsTr (^e{p)P{t)9{p)pe{p)P{s)e{p)^ . (6.36) 

The first hne is the decoherence functional for crossing from the right, and the second line is 
the decoherence functional for crossing from the left. This means that a sufficient condition 
for decoherence of the crossing probabilities is that our initial state is a superposition of a left 
moving and a right moving wavepacket, each of which satisfy the conditions for decoherence 
of their corresponding arrival time probabilities. But this is exactly what we set out to prove. 
Detailed calculation of the two terms in Eq.(6.36) has been carried out in Chapter 5, and we 
simply quote the result, that wc have decoherence for gaussian wavcpackcts, and orthogonal 
superpositions of these, provided they are tightly peaked in momentum, and subject to the 
requirement 

{t2 - U)E » 1, (6.37) 

where E is the mean energy of the wavepacket. For more general wavepackets we require 
an environment to ensure decoherence of the arrival time probabilities [100] and we suspect 
that the same will be true here. This issue will be addressed elsewhere. 



6.6 Summary and Discussion 

Our aims in this chapter were twofold. Firstly we wished to understand whether the origin 
of the arrival time class operators might be understood by some intuitive argument. We 
saw that this is indeed the case. Secondly we wished to see how the class operators for the 
arrival time problem might be extended to deal with wavepackets incident from both sides of 
the origin. Again we saw that this is possible and that the probabilities defined in this way 
have the nice feature of being insensitive to the interference terms involving left and right 
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moving wavepackets. What is interesting about the analysis of crossing times is that it is not 
obvious how to do this in the language of complex potentials, yet it seems relatively simple 
in terms of the class operators defined via strings of projection operators. In particular it 
is gratifying that the use of decoherent histories provides a justification for some simple 
semiclassical manipulations, such as the derivation of Eq.(6.22), and this suggests it may be 
usefully employed to address the definition of other time observables. This will be pursued 
elsewhere. 



118 



7 Quantum Arrival Time for Open 
Systems 



Tempora mutantur, et nos mutamur in 
illis. 

Times change, and we change with them. 

WilUam Harrison 

7.1 Introduction 

In the previous two chapters we have discussed how the arrival time distribution may be 
derived in the framework of decoherent histories. We found that the arrival time distribution 
may, under the conditions of dccohcrence, be written as the integral of the current at the 
origin, thus providing justification for the semiclassical expression we wrote down in Chapter 
1, Eq.(1.20). Recall this was, 

p{h,h) = Tr{P{ti)p)-Tr{P{t2)p)^- dtTriPpt) 

dtJ{t), (7.1) 

L 

with 

J{t) = [ dpdqtP}Mw,{p,q), (7.2) 

where P = 0{x), we denote the state at time t by pt and the Wigner function Wt{p,q) is 
defined by Eq.(1.21). Recall also the three problems with Eq.(7.1), 

1. It is not generally positive, even when it would be expected to be so classically. 




119 



2. It was not derived but postulated. 

3. It is not of the general form of the trace of the density matrix times a projection 
operator or POVM, Eq.(1.29). 

Chapters 5 and 6 sought to address these issues by deriving the arrival time distribution 
from decoherent histories and showing how it reduces to Eq.(7.1) in the case where we have 
decoherence. We succeeded in addressing the first and second points, but no mention was 
made of the third. In addition wc generated the new problem that we could only demonstrate 
decoherence of histories for particular initial states. 

In this chapter we look at the definition of arrival times for open quantum systems, that 
is for a particle coupled to some environment. One aim of this is to extend the decoherent 
histories analysis to a general initial state and we discuss this in more detail below. However 
the first aim of this chapter is to see whether the semiclassical expression Eq.(7.1) is any 
less problematic for a particle coupled to an environment. An interesting clue is provided 
by the expression for the current in terms of the Wigner function, Eq.(1.22). Negativity 
of the Wigner function is a neccessary condition for the negativity of the current (in the 
sense of backflow. as discussed in [49]), but it is known that evolution in the presence of 
an environment typically renders the Wigner function positive after a short time [21]. This 
suggests that something like Eq.(7.1) may be an acceptable, if still heuristic, arrival time 
distribution for a system coupled to an environment, at least after some time. The first 
aim of this chapter is to derive the correct analogue of Eq.(7.2) for a system coupled to an 
environment, and to prove that it can indeed be regarded as an arrival time distribution 
after a short time. In particular we show that the arrival time probabilities computed in this 
way may be written as the trace of the density matrix times a POVM. This then addresses 
the first and third points above. 

However although this may allow us to interpret the current as the arrival time distribution, 
we are still left with the task of deriving Eq.(7.1) from some more fundamental quantity. 
This was achieved in Chapters 5 and 6 in the context of the decoherent histories approach 
to quantum theory [30, 31, 37, 85, 41, 22], for the case of a free particle and the derivation 
was shown to hold for states and intervals exhibiting sufficient decoherence. However it was 
also shown that there exist states for which decoherence of histories cannot be obtained 
for reasonable levels of coarse graining, and for which therefore we cannot assign a time of 
arrival probability distribution. States exhibiting backflow are good examples of this. 
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Lack of decoherence for a general state is a situation frequently encountered in the litera- 
ture on decoherent histories. The solution lies in the observation that realistic systems are 
always coupled to an environment and that as such they are most fundamentally described 
by open system dynamics. Because by definition an environment consists of degrees of free- 
dom about which we have no knowledge, and over which we have no control, it is natural 
to coarse-grain over these degrees of freedom. Such coarse-graining generally results in the 
recovery of approximate classical behaviour, and thus decoherence. We therefore anticipate 
that decoherence of histories corresponding to arrival times may be achieved for a generic 
state, provided that particle is coupled to a suitable environment. The second aim of this 
chapter is to examine this scheme. We show that the probabilities calculated in this way are 
approximately decoherent and approximately equal to those computed via the analogue of 
Eq.(7.2), valid in the case of an environment. 

Once we introduce an environment, however, the correspondence is with a classical stochas- 
tic theory, rather than with a deterministic one. The classical trajectories may now cross 
the origin many times due to fluctuations. The arrival time distribution for such a classical 
theory was computed in Ref.[54], and is given by 



where q) is the initial phase space distribution, evolved with a type of restricted prop- 
agator valid in the case of an environment, with boundary conditions, [54] 



In Rcf. [54] Halliwcll and Zafiris presented a decoherent histories analysis of the correspond- 
ing quantum case. Although the conclusion they reach is sensible their analysis actually 
contains a small error. This was the result of a lack of appreciation of the role of the quan- 
tum Zeno effect in these calculations, as outhned in Chapter 1 Section 1.6. We will show 
in this chapter how the analysis can be corrected to give a correct expression for the arrival 
time probability in decoherent histories. We shall be interested in a particular limit of this 
general case where the stochastic trajectories are sharply peaked about the deterministic 
trajectories we would have in the absence of any environment. In this case the restricted 
propagator may be replaced with an unrestricted one, and we will show that quantisation in 
this limit yields an expression of the form Eq.(7.1), but with the Wigner function evolved in 




(7.3) 



wl{p, 0) = 0, for p > 0. 



(7.4) 
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the presence of an environment. 

It will eventually be necessary to specialise to a particular model of system-environment 
coupling in order to obtain quantitative results, but we will begin by considering general 
models which can be written in terms of a master equation of Lindblad form. This kind 
of system has been extensively studied in the decoherence literature and the properties 
of such evolution, such as suppression of interference effects, loss of entanglement and an 
approximate recovery of classical behaviour have been discussed in [21, 56, 55, 33]. 

This chapter is arranged as follows, in Section 7.2 we explore some properties of the arrival 
time distribution for general open systems of Lindblad form, our aim being to derive the 
corrections to the current resulting from the environmental dynamics. In Section 7.3 we 
discuss quantum Brownian motion and specialise the results of Section 7.2 to this case. In 
Section 7.4 we derive some properties of the arrival time distribution for quantum Brownian 
motion, in particular we prove that the current becomes positive after a finite time. In 
Sections 7.5 we briefly recap the expressions for the arrival time distribution that arise in 
decoherent histories, in particular we show that the decoherence condition is closely linked 
to the semiclassical approximation. In Section 7.6 we then examine whether histories corre- 
sponding to arriving in different intervals of time are decoherent. We summarise our results 
in Section 7.7. This chapter is based on Ref.[100]. 

7.2 Arrival Time for Open Quantum Systems 

We consider an open quantum system consisting of a free particle coupled to an environment 
with master equation of the following Lindblad form [69] , 



where Lj are the Lindblad operators and H — Hq + Hi is the free hamiltonian plus a possible 
extra term that depends on the Lj. (This term gives rise to frequency renormalisation in 
oscillating systems.) (In this chapter we explicitly write factors of h to help us see the 
difference between classical and quantum terms in equations.) Specific forms for the Lj may 
be computed for particular models, but for the moment we leave this general. The Lindblad 
equation represents the most general master equation for a Markovian system which preserves 
the properties of the density matrix, in particular positivity. 

We now extend the derivation of Eq.(7.1) to this system. The probabihty of crossing 



dp 
dt 



(7.5) 
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during the interval [^1,^2] is 



t2 



^ "^^^^ ^^^^^^ ^ ^i^)p]pt) dm {[H,, p]pt) 

+ 1 E ^^Tr ([L], P]L,p, + Lt[P, L,]p,) , (7.6) 

where P = 6{x). The first term is the standard current expression, although with the state 
evolved according to Eq.(7.5), and we therefore recover Eq.(7.1) when all the Lj are 0. The 
second and third terms depend on the Lindblad operators, Lj and thus on the form of the 
system-environment coupling. 

To proceed further we specialise to the case where L is a linear combination of x and p, 
L — ax + ibp, where a and b are real constants. The master equation Eq.(7.5) is then 

-Q^ = -^[Ho.pt] -iab[x,{pup}] - -^[x,[x,pt]] - ^[p, (7.7) 

Note that this equation is also identical in form to the exact master equation for a particle 
in a gas environment, given in [20]. 

To derive the arrival time distribution we could simply substitute L = ax + ibp into 
Eq.(7.6), but the algebra is somewhat clumsy, and there arc a number of terms which vanish 
for reasons not immediately apparent from this expression. An equivalent approach is to 
start from Eq.(7.7), multiplying this expression by P and taking the trace. Because the 
second and third terms on the right of this expression have the form [x, ^4] their contribution 
is proportional to 

Tr([£, i]P) = Tr(i[P,£]) = 0, , 

so that only the first and last terms on the right hand side of Eq.(7.7) contribute. The final 
term on the right of Eq.(7.7) gives a contribution, 

-|tV([p, [p,p,]P) = -|tV([p, [P,p]]p,) = -ihb^TT{\p5{x) - 5{x)p\pt), (7.8) 
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so we arrive at the expression, 



2m 



-1 




tl 



dm{\pd{x) + d{x)p]pt) 



dtTr {\pS{x) — S{x)p\pt) . 



(7.9) 



The second term in this expression has a somewhat unusual form and its significance is 
not immediately clear. We will see below that this term is related to diffusion in position. 
There is a connection here to a recent paper by Genkin, Ferro and Lindroth [32], in which 
the authors sought to examine the effects of an environment on the arrival time distribution. 
In that paper the authors used the standard expression for the current valid in the unitary 
case, Eq.(7.2), ignoring the extra terms that arise because of the environment. Although this 
may be a good approximation when we can neglect the effects of dissipation, it is clear that 
there may be significant corrections to the current for strong dissipation. They also pointed 
out that these corrections are equivalent to the presence of extra terms in the continuity 
equation, although they did not compute these explicitly. For the sake of completeness, and 
also because it helps to understand the nature of the extra terms in Eq.(7.9) we will derive 
them here. 

To derive the continuity equation we multiply Eq.(7.7) by 6{x — x) and perform the trace. 
If we neglect the final three terms on the right we arrive at the standard result. 



with J{x,t) defined by an obvious extension of Eq.(7.2). Turning to the extra terms that 
result from the inclusion of the environment, the second and third terms vanish in exactly 
the same way as for the current above, and the correction term is therefore given by. 



dpi 
dt 



{x,x) + -^{x,t) = 0, 



(7.10) 



— Tr([p, \p, pt]]5{x - x)) = 2h%'^-^{x,x), 



(7.11) 



so the continuity equation now reads. 



dpt 
dt 



{x,x) + — {J{x,t) + JD{x,t)) = 0, 



(7.12) 
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where we have identified the diffusive current, 

Jd {x, t) = -2h^b' ^ {x, x). (7. 13) 

This is a specific example of a modification to a conservation equation resulting from open 
system dynamics. For a more general discussion of these issues see [29]. 

The analysis presented above can also be phrased in terms of phase space distributions. 
The Wigner function corresponding to pt is defined as [10] 

Wtip, Q)-^J d^e-i^^ptiq + e/2, q - i/2). (7.14) 
The Wigner transform of the master equation Eq.(7.7) is 

The first term is the standard unitary term, whilst the second term represents dissipation, 
and the third and fourth terms represent diffusion. The arrival time distribution can be 
written in terms of the Wigner function by taking the Wigner transform of Eq.(7.9) 

p{tut2) = J^t J dpdq (^tEls{q)Wtip,q) + ^^(^) ^(^^' «)) • (^-^S) 

Eqs.(7.9) and (7.16) are the sought for generalisation of Eq.(7.1) for the case of a particle 
coupled to an environment. 



7.3 Quantum Brownian Motion 

Quantum Brownian motion [14, 56, 55] is a commonly used form for the environment of an 
open system, partly because it is exactly solvable and partly because in many cases it is a 
good approximation to a realistic environment. For the quantum Brownian motion model 
[56] we have Eq.(7.7) with 




(7.17) 
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Here D — 2m'ykT, where T is the temperature of the environment, and 7 is a phenomeno- 
logical damping constant, see [14]. Eq.(7.7) may be written as 

dpt, . ih f \ , . D . . X 

f d d\ yi^'f f d 9 V 

-7(x-y)(^^-^jM^,y)-^(^^-^j Ptix.y)- (7.18) 

Although we could in principle work with this general case, it is useful to specialise to the 
case of negligible dissipation. There are two reasons for this, the first is that the analysis 
is considerably simplified, and this helps us to see the important effects more clearly. The 
second is that we have a particular aim in mind here, and that is to understand how a sensible 
classical result emerges from the quantum case. The classical case we have in mind is that 
of a heavy particle following an essentially classical, deterministic trajectory, but subject to 
small quantum fluctuations. We therefore restrict our analysis to timescales much shorter 
than the relaxation time 7"^, so we can drop the final two terms in the master equation 
above. 

This master equation may then be solved in terms of the propagator [55] , 
J{x,y,t\xo,yo,Q) 

Taking the same limit in the equation for the Wigner function, Eq.(7.15), gives 

dWt ^ pdWt ^ ^d'Wt 
dt m dq dp^ 



J dxodyoJ {x, y,t\xo, yo, O)po{xo, yo) 



(7.19) 



m 



2nhtJ \2ht 



im 



[{x - xof -{y- yo^] 



Dt 

T 



[{x - yf + {x- y){xo - yo) + (xo - yoY]) ■ (7.20) 
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Evolution of the Wigner function may also expresed in terms of a propagator [55] , 

Wt{p,q) 
K{q,p,t\qo,po,0) 

where N, a, /3, and e are given by 

For later convenience we note that we can make the simple change of variables here, 
qo ^ qo ~ Pot/'^n-i so that 

^t{p,q) = J dpodqoK{q,p,t\qo,po,0)Wo{po,qo) (7.25) 

K(q,p,t\qo,po,0) = N ex.p {-a(p - pof - 13 (q - qof + e(p - po)(q - qo)) 

W{p,q) = W{p,q-pt/m). (7.26) 

This form makes it clear that the evolution consists of two effects. The first is a shifting 
along the classical trajectories, whilst the second is a spreading in phase space. 

It is useful to consider this process in more detail. In the presence of an environment the 
width of the momentum distribution becomes time dependent, and we have, 

(Ap),2 = {Ap)l + Dt, 

where (Ap)o is the momentum width of the initial state. We recognise an immediate difficulty 
here. Even for initial distributions consisting entirely of left moving momenta, Wt will 
develop support on p > under evolution. This means we cannot strictly regard the current 
as an arrival time distribution, since a typical trajectory will now cross the origin many 
times. Differently put, the arrival time distribution is strictly a measurement of the first 
passage time and this is no longer simply related to the current. (Note that although this 
result is similar to the problems created by backflow, the reasons behind it are very different. 
The spreading of momentum induced by evolution in an environment is a purely classical 



dpodqoK{q,p, t\qo,po, 0)Woipo, qo) 



(7.22) 



A^exp -a{p - pof - P [ q - qo 



Pot 



m 



+e(p -Po){q-qo- 



'Pot_ 
m 



(7.23) 
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effect.) 

However, all is not lost. Although the current is no longer strictly the arrival time distribu- 
tion, it may be a very good approximation to it. This is because the deviation of the current 
from the "true" arrival time distribution will be related to the probability that the state has 
the "wrong" sign momenta. This means that, provided we are in the "near-deterministic" 
limit, 

m't«Pl (7-27) 

where po is the momentum around which the initial state is tightly peaked, the current will 
still be a very good approximation to the true arrival time distribution. The timescale on 
which this analysis breaks down is given by the "stochastic" time. 

After this time we must therefore revert to using the exact expression for the arrival time 
given by Eq.(7.3). We see from the definition of D that, 

for the states we are interested in. This means the stochastic time is much longer than 
the relaxation time 7^^, so that working in the near-deterministic limit does not impose any 
further constraints compared with neglecting dissipation. 

Returning to our arrival probability, in the limit of negligible dissipation Eq.(7.9), becomes, 

p{t2,h) = r dtJit) (7.29) 

J{t) = [ dqdp^-^6{q)Wt{p,q). (7.30) 

This expression is now identical to the unitary case, Eq.(7.2), but with the Wigner function 
evolved under quantum Brownian motion. 

We now turn to the question of whether inclusion of an environment helps ensure the 
positivity of Eq.(7.29). 
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7.4 Properties of the Arrival Time Distribution in 
Quantum Brownian Motion 

In the introduction to this chapter we noted that evolution in an environment typically 
renders Wt positive after a short time. We now wish to examine the effect of this on our 
candidate arrival time distribution Eq.(7.29). To this end we introduce the notation [21] 



(7.31) 



and also the class of Gaussian phase space functions 

^(z; A) = ^-i^ exp (-^z^A-^z) , (7.32) 

where A is a 2 x 2 positive definite matrix, with determinant \A\. g{z;A) will be a Wigner 
function if and only if 

\A\ > -. (7.33) 

A useful result is that 

J d^zg{zi - z; A)g{z - Z2; B) = g{z^ -z^^A + B). (7.34) 
In this notation we can write the propagation of the Wigner function, Eq.(7.22), as 

Wt{z) = j dz'g{z - z';A)Wo{z'), (7.35) 



where 



Since 



after a time 



A^Dt{ ^ l^"^ J. (7.36) 



\A\ ^ ^ (7.37) 
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g[z — z'; A) will be a Wigner function, and thus Wt{z) will be positive because it is equal to 
the convolution of two Wigner functions. Here ti = \J 2mh/ D is the localisation time. 
This is a useful result. Expressing the current Eq.(7.30) in this new notation, 



J{t) ^ dz 



m 



S{zi)Wt{z), 



(7.39) 



we see that after the time given in Eq.(7.38), because Wt > the current Eq.(7.39) will be 
positive if the state consists purely of negative momenta. This means that after this time 
Eq.(7.29) will be a positive arrival time distribution. This holds provided times involved are 
much smaller than Tg, as per the discussion below Eq.(7.26). 

Now wc turn to examining the properties of the current. We wish to find an expression 
for the current in the form Tr(Pp), where P is a projector or POVM. This would allow 
us to express the heuristic arrival time distribution, Eq.(7.1), in the same form as other 
probabilities in quantum theory. Starting from Eq.(7.39) it is useful to write. 



j(t) 



dz 



m 



S{zi)g{z-z';A)Wo{z'), 



(7.40) 



using Eq.(7.35). We can deconvolve the propagator into two gaussians using Eq.(7.34), in 
particular we will let A = Aq + B, where Aq is a minimum uncertainty gaussian, and B is 
the remainder. 



Ao^h 



' 

l/4s^ 



B^Dt 



-Ao. 



(7.41) 
(7.42) 



2 t/m 

Here s is some real number. Using the convolution property Eq.(7.34) we can write the 
current Eq.(7.39) as 



J{t) = j dz" J dz^-^S{z^)g{z - z"; B) J 



dz^ U{z^ + Zot/m)g{z - z"; B) 



m 



dz'Wo{z')g{z" -z';A^) 



(7.43) 



where Q{z) is the Q-function [90], and we have undone the change of variables implied in 
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Eq.(7.26). The Q-function can be written as 

g(z) = - (z|p|z) . 

TT 

We can therefore express the current as, 

J(i)=Tr(Fp), (7.44) 

where 

F = ^ y dz" \z") (z"| (^J dz^-^S{zi + Zot/m)g{z - z"; B)^ 

= [ dzP^^^^S(zi + zot/m), (7.45) 
J rn 



and we have defined the POVM 



P^-\j dz" \z") {z"\ g{z - z"; B), (7.46) 

which is clearly a phase space operator localised around z. F is therefore a smeared version 
of the object used to compute the current classically, —p5{xt)/m. This holds for times 
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Tl<t« Ts. (7.47) 



Assume for a moment that 5 = 0, and so = |z) (z|. The current would then be given 

by 

J(t) = / dz^^^S(zi - Zot/m)Q(z). (7.48) 
J rn 

Since the Q-function is positive by construction the current computed in this way will be 
positive to the extent that the Q-function has vanishing support on p > ((5(z) cannot 
strictly vanish for p > but it can be exponentially small, and this suffices here. See the 
comments below Eq.(7.26).) For S > the Q-function is simply smeared further, and this 
will preserve the property of positivity. In fact one could imagine postulating Eq.(7.48) as 
the definition of the arrival time even in the unitary case, since it clearly satisfies all the 
conditions required. 

The probability of arriving in the interval [ti, t-^ is given by the integral of Eq.(7.44), and 
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because is time dependent (via B) this will not in general have a simple form. However 
there is a separation of time scales here, the time scale of evolution of is seen from 
Eq.(7.37) to be r;. So, if ^2 — << ti we have approximately 



p{tiM)^ / dm{Fp)^ / dzTr pP^ 




S{zi + zot/m) ^Tr{Ep), (7.49) 



where 



dzP49{zi + zoti/m) -9{zi + ZQti/m)] 



(7.50) 



is a POVM representing arrival at x = between ti and t2- Compare this expression with 
Eq.(7.1), which may be written. 



The two expressions are very similar, but crucially Eq.(7.50) is positive for < (ie p < 0). 
In this case therefore, the effect of the environment is simply to smear the unitary result 
over a region of phase space, given by P^ computed at ti. 

Eqs.(7.49) and (7.50) form the first significant result of this chapter. Eq.(7.49) expresses 
the heuristic arrival time distribution, Eq.(7.1), as the trace of an operator times a POVM, 
and thus has the same form as standard expressions for probability in quantum theory. The 
POVM, Eq.(7.50), arises because the environment effectively measures the system. 

We have therefore discovered a range of times for which the current, Eq.(7.39), gives a 
positive arrival time distribution. After a time of order t; interference effects vanish and 
we can regard Eq.(7.39) as the arrival time distribution, which can also be written in the 
form Eq.(7.44). Eventually, however, on a time scale of order the environment causes 
diffusion in momentum to such an extent that the trajectories of the particle are no longer 
sharply peaked around the classical trajectory computed in the absense of an environment. 
Since the particle is still behaving classically after this time there will exist an arrival time 
distribution of the form Eq.(7.3), but our simpler expression Eq.(7.39) will no longer be a 
good approximation to this. It is easy to show that Ts/ti — Eti/U where E is the energy of 
the initial state, and thus these expressions are valid for an large interval if Eti » H. 



p{h,h)^Tr{[e{xt,)-0{xt,)]p). 



(7.51) 
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7.5 The Decoherent Histories Approach to the Arrival 
Time Problem 



We turn now to the question of defining arrival times for open systems via decoherent 
histories. In Chapters 5 and 6 the class operators corresponding to a first crossing of the 
origin in the time interval [tk-i, tk] were claimed to be 

Ck = P{tk)P{tk-i) ■ ■ ■ P{t2)P{ti), (7.52) 

for k > 2, with Ci = P{ti) and where P{t) = 9{—Xt). These clearly describe histories which 
are in a; > at times ti,t2, ■ ■ -tk-i and in a; < at time tk, so, approximately, describe a 
first crossing between and t^. 

In Chapter 6 these class operators were then simplified with the help of the following 
semiclassical approximation 

p{tn)---p{h)p{hm^p{tnm. (7.53) 

Prom this we obtain the class operators for first crossing between tk-i and tk as 

Ck ~ P{tk)P(tk-i) (7.54) 
= P{tk-i) - P{tk)P{h-i) 

^ P{tk-i)-P{tk), (7.55) 

where we use the semiclassical approximation again to arrive at the last line. 

The key step in making contact with the simple heuristic result of Eq.(7.1) was the ap- 
proximation of the class operators Eq.(7.52) by Eq.(7.55). We now want to further explore 
this approximation, and prove that it is in fact closely related to the decoherence condition. 
The aim here is partly to make the the argument in Section 6.5 more precise and partly to 
check whether the inclusion of an environment adds anything new. 

For a free particle, the error in the semiclassical approximation comes from spreading in 
momentum caused by the chopping of the wavepacket. We have seen in Chapter 4 however, 
that for projections spaced by a time greater than h/ E the probability of the measurements 
causing the particle to be refiected is very low. 

We want to find some way of quantifying the error induced by this approximation, for 
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this purpose it is sufficient to consider the approximation as used to go from Eq.(7.54) to 
Eq.(7.55). Writing this in terms of the density matrix, 

Tr{P{tk+,)P{h)pP{h)) = Tr(P(t,+i)p) - Tr(P(tfc+i)P(t,)p) - Tr{pP{h)P{h+i)) 

+Tr{P{h+i)P{tk)pP{tk)) 
^ IV(P(ife+i)p), 

where the last line defines the semi-classical approximation. Noting that, with > tk 

\TT{P{t^)P{tk)p)\' < Tr{P{tra)P{tk)pP{tk)). (7.56) 

A sufficient condition for the validity of the semi-clasical approximation is therefore that the 
object 

^k,m := 1i{P{tm)P{tk)pP{tk)) (7.57) 

is much smaller than 1 for all tm> t^. 

Now assume for a moment that the semiclassical approximation holds, so that the class 
operators are given by Eq.(7.55). The probability of first crossing between and t^+i is 
given by 

Tr(CfcpCj) = Tr(C,p) + 2Tr(pP(tfe+i)) - Tr(pP(t,)P(t,+i)) - Tr(P(t,+i)P(tfc)p). 

However since we are assuming that the semi-classical approximation holds, the last three 
terms cancel and we obtain 

IV(CfcpC;) = TV(Cfcp) = TV(pCS. 

Furthermore, consider a general off-diagonal term in the decoherence functional, where with- 
out loss of generality we take ij^+i < tm 

Tr(CfcpCi) = TY{P{tm)P{tk)p)-T^{P{tm)P{tk+i)p) 

-Tr{P{tm+i)P{tk)p) + Tr{P{tm+i)P{tk+i)p). (7.58) 

We see that this vanishes if the semi-classical approximation holds. 

What we have shown therefore, is that the semi-classical condition and the decoherence 
conditions are closely related. If the object given in Eq.(7.57) is much smaller than 1 then 
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both conditions are satisfied. This means firstly that the class operators are approximately 
given by Eq.(7.55) and secondly that the histories described by these class operators are 
approximately decohcrcnt. 

In general, the probabilities for time of arrival as computed by decoherent histories differ 
from the heuristic ones obtained by the standard quantum mechanical analysis. We don't ex- 
pect agreement since, as discussed in the Introduction to this chapter, the heuristic formula, 
Eq.(7.1), is not in general of the canonical form required for genuine quantum mechanical 
probabilities. It is also the case that decoherent histories only ascribes probabilities to certain 
sets of histories. However if the object Ak^m is small, then we have shown that the decoherent 
histories analysis reproduces the probabilities computed in the standard, heuristic, way. 

As well as clarifying the analysis of [50] and Chapter 6, note that all of the statements 
above apply to the case of a particle coupled to an environment. In this case the trace is to 
be taken over the system and environment, and P{tk) = P{tk)s ® le is a projection onto the 
degrees of freedom of the system, s, tensored with the identity operator on the environmental 
degrees of freedom, e. This is the case of interest in this chapter. 

Since the object A^^m plays a central role in our analysis it is interesting to ask if it has 
any physical interpretation. The answer is that it does. Recall that our class operator for 
crossing between tk-i and tk is given by Eq.(7.55). In Chapter 6 Eq.(6.19) we noted that 

P{h-i) - P{h) = P{h)P{h-i) - P{h)P{h-i). (7.59) 

The approximation leading from Eq.(7.54) to Eq.(7.55) is equivalent to dropping the second 
term on the right hand side of Eq.(7.59). The error in this approximation can be estimated 
by computing the probability associated with the class operator P{tk)P{tk-i)- We see 

Tr{[P{tk)Pitk-iMP{tk)P{tk-i)]'') = Tr{P{tk)Pitk-i)pP{tk-i)) = Ak-i,k. 

However Eq.(7.59) has a simple physical interpretation: it is the decomposition of the total 
current into right and left moving parts. The object Aj. „j is therefore just the probability 
associated with the right moving current. Classically this is small by construction, since 
we will choose to work with wavepackets tightly peaked in negative momentum. Quantum 
mechanically however this term need not be small, indeed the existence of the backflow 
effect [12] shows that this term can sometimes be larger than the left moving current, even 
for wavepackets composed entirely of negative momenta. 

The semiclassical condition we are imposing, and by extension the decoherence condition. 
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is therefore stronger than the condition for the current to be positive. Decoherence requires 
the absence of interference between crossings at different times, whilst the standard heuristic 
analysis includes some interference effects, provided they are not so large as to render the 
arrival time probabilities negative. This is a consequence of the fact that although for a 
suitably semiclassical state histories peaked about the classical trajectories are decoherent, 
the set of sets of decoherent histories contains many more histories than simply the classical 
ones. The significance of this will be explored elsewhere. 

7.6 Decoherence of Arrival Times in Quantum 
Brownian Motion 

7.6.1 General Case 

In the previous section we have shown that an arrival time distribution may be derived from 
decoherent histories and that it agrees with the current provided there is decoherence. We 
therefore turn now to the question of determining for which states and intervals decoherence 
is achieved. Recall that our aim is to show that the inclusion of an environment gives rise 
to decoherence of arrival time probabihties for generic initial states of the system. We will 
start our decoherent histories analysis with a discussion of the general case considered in 
Ref.[54], but most of our detailed results will concern the near-deterministic limit discussed 
in the introduction and in the discussion of the Wigner function above. 
The decoherence functional may be written in path integral form as 

V{a,a') ^ jvx j^Vye^ {^-S[x] - '-S[y] + iW[x,y]^ po{x,y), (7.60) 

where a, a' represent the restriction to paths that are, for example, in x > at and in 
X < Q aX t2. W[x,y\ is the Feynman-Veron influence functional phase which summarises the 
effect of the environment, and is given in the case of negligible dissipation by 



W[x,y]^ 




(7.61) 



It is this phase which is responsible for the suppression of interference between paths x{t) 
and y{t) that differ greatly, and produces decoherence. 

The decoherent histories analysis we are about to perform was first attempted by Halliwell 
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and Zafiris in Ref.[54]. Their conclusions are reasonable, but the analysis leading to them in 
fact contains a small error. This error arose due to a lack of appreciation of the role of the 
Zeno effect, as discussed in the Introduction. We aim to show here how the analysis may 
be modified in line with the treatment of the arrival time problem presented here and in 
Ref.[50]. 

We start from the expression for the density matrix for states that do not cross the origin 
in a time interval [0,t] (Ref.[54], Eq.(4.34)) 



Prr{xf,yf) ^ J^Vx J^Vyexp J ^^i^^ - f) - ^ Jdt{x-yf 



po{xo,yo). (7.62) 



Here the restriction is interpreted as x, y > at times e, 2e.... The probability that the state 
does not cross the origin in this interval is then given by the trace of this expression. There 
are two interesting limits to be taken here. The first is letting e — )■ 0, and this recovers the 
notion of the restricted propagator. The second limit is that of strong decoherence effects, 
so that we can assume that the path integral is tightly peaked in {x — y), this recovers the 
classical limit. The claim in Ref.[54] is that we can take these limits simultaneously. To see 
why this is problematic we write the restricted path integral as a product of propagators, 

POO 

Prr{xf,yf) = Y[ dxkdykJ{xk+i,yk+i,tk+i\xk,yk,tk)po{xo,yo), (7.63) 

where tk — ke, to — 0, tn+i — t. Changing variables to X — x + y, ^ — ^{x — y) gives 

POO pAk / m \ 

Prr{xf,yf) = 11 dXk d^ki^-r-jexp 



k=0 

De 



-^[(Cfc+1 - ik){Xk+\ - Xk)\ 



2[^fe+l + ^k+lik + Cfc] 



3^ 



Po(Xo,eo). (7.64) 



Now in order for the product of propagators to be equal to a restricted propagator, we need 
to take e — > 0. However to recover the classical result we need to replace the limits of 
integration for with ±oo, which requires that De » /l"^ where I is some length scale. 
Clearly these two limits are incompatible for finite D. 

We see the difficulty here immediately, trying to take e — > will lead to the Zeno effect. 
We know from Chapter 4 that for a free particle projections separated by times greater 
than Ti/E will not give rise to the Zeno effect, but the issue is more subtle here, since we 
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have non-unitary evolution. One would expect that the time scale on which projections 
give rise to reflection to be reduced in the present case, since a generic system coupled to 
an environment tends to display classical behavior and the Zeno effect should vanish in the 
classical limit. The claim in Ref.[54] is essentially (although they did not express it in this 
way), that the inclusion of the environment is enough to eliminate the Zeno effect entirely. 
This seems unlikely to be true. This issue will be taken up elsewhere [52] . The more general 
point is that it is not possible to take e — )> and still get a good account of arrival times, 
we should always expect to have to settle for some minimum resolution. For the purposes 
of recovering classical behaviour, however, this is of no significance, since the timescale on 
which a classical "wavepacket" would cross the origin is ~ h/AH and this is much longer 
than h/E. We therefore do not loose anything classically by leaving e small but not zero. 

In any case it seems clear that there exists a short timescale on which we can still approx- 
imate the limits in the ^ integrals above in the manner indicted, but that is short enough to 
give non-trivial crossing probabilities. Following the steps in Ref . [54] then gives 



p{ti,t2)^l-Tr{prr)^ dt f dp! dq^-^d{q)W;{p,q), (7.65) 
where the Wigner function, W[ obeys, 

W[{p, Q) = J dpodqoKr{p, q,t\po, qo, 0) 1^0(^0, qo), (7.66) 
where Kr is the "restricted" propagator, defined by 

Kr{p,q,t\po,qo,0) = Y[ / dpk / dqkK{pk,qk,tk\pk~i,qk-i,tk-i), (7.67) 
k=i 

with tk — ke. However because we are now dealing with an essentially classical system, this 
propagator is well approximated by its continuum limit. Understood in this sense, we see 
that the conclusions of Ref. [54] are in fact correct, even if the argument leading to them is 
not. The interesting subtleties that arise in this argument relate to the classical limit of the 
Zeno effect, this will be explored elsewhere [52]. 

As important as this general case is, it is of interest to examine the simpler case of near- 
deterministic evolution, where the classical limit of the arrival time probability is simply 
expected to be the time integral of the current density. The analysis in this case simplifies 
considerably and we will explicitly exhibit decoherence of histories for suitable intervals. 
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Figure 7.1: The probability that a generic wavepacket arrives between ti and t2 can be ex- 
pressed in terms of the state at ti or, via the propagator, in terms of an initial 
time t — 0. 

7.6.2 The Near-Deterministic Limit 

We have previously seen that, in the near-deterministic limit, if there is decoherence then 
the arrival time probabilities derived from decoherent histories agree with those computed 
from the current. Wc therefore turn now to discussing the conditions under which we have 
decoherence of crossing probabilities. The general picture wc have in mind is illustrated 
in Fig. (7.1), we have an initial wavepacket defined at t = 0, evolved in the presence of an 
environment until time ti, and then we wish to compute the probability of crossing between 
ti and t2. 

Our task is to compute the off-diagonal elements of Eq.(7.60), and this will be rather 
involved in general. However we arc saved from having to do this by the observation in Section 
7.5 about the relation between the scmiclassical condition and the decoherence condition. It 
suffices therefore to compute the quantity A^.m defined in Eq.(7.57), in the presence of an 
environment. We anticipate that this will be small simply from the form of Eq.(7.60). This 
is because the effect of the environment is to cause the density matrix to become tightly 
peaked around the classical path and classically, for p < the probability given by Ak,m is 
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zero. We will see how this works in a specific example below. 

We take /c = 1, m = 2 without loss of generality, and we drop the subscript from now on 

A = Tr{P{t2)P{ti)pP{ti)) ^ jvx jvye^v i^^S[x] - '-S[y\ -^J dt{x - yf^ PtA^i,yi), 

(7.68) 

where the histories a are those that start at Xi, yi < and finish at X2 > 0. In terms of the 
density matrix propagator, Eq.(7.20) 



/■oo pO rO 

A = / dx2 dxi / dyiJ{x2,X2,t2\,Xi,yi,ti)pti{xi,yi) 

Jo J — oo J —oo 



\Trht) 



/•oo /-O pO 

/ dx2 / dxi / dyi ( 1 

^0 J —oo J —oo 

( II \2 \2\ -^(^2 ~ ^l) / \2 \ / 

exp ( ^^^^^3^((X2 -xi) -[x2-yi)) ^ K^x-Vx) \ptr[xx,yx). 



Transforming to new variables 



^ = ^4^, e = ^i-?/i, (7.69) 



and writing the density matrix in terms of the Wigner function via 

p{x,y) = dpe^^(--^)iy(p, ^), (7.70) 

we obtain 



A = - /" / dX I di\ dp( ^ 

Jo J-oo J-X J-oo \ 



2TTh{t2-h] 

We see from this expression that there is a time scale, t/ = ^/2mh/D, set by the environment 
on which localisation effect are important. This timescale is the same as that on which the 
current becomes positive, as we saw earier. 

There are three cases to explore here, the first is where there is no environment, D — 0. 
This is the case covered in Ref.[49]. The second case is the intermediate one, h/ri » 1 
but {t2 — ti)/Ti « 1. This is the most general case in which we can expect to have 
environmentally induced decoherence. The final case is where h/ri, {t2 — ti)/Ti » 1. This 
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is the case of very strong environmental coupling. 

For the free case, D — the integrals over ^ and X2 in Eq.(7.71) may be carried out to 
give, 



where 



/oo 
dpWt,iX,p)f 
-oo 



X / mX 



dy 



h \{t2-h) 



smy 



+ p 



(7.72) 



(7.73) 



This expression for A is now identical to Eq.(5.96) of Chapter 5 and the same conclusions 
apply. We briefly repeat the analysis here however, since it turns out to be relevant for the 
second case. 

We note flrstly that /(+oo) = 0, /(— oo) = 1, /(O) = 1/2. Now we assume that our state 
at is of the form 



W,,{X,p)^-e.^[^ 2^^ ^(P-Po) ), 



2(72 



(7.74) 



and that cr^ is large, so the state is tightly peaked in momentum. We can therefore integrate 
out p, setting p — poio obtain. 



-f 

J — c 



dX 



V27r(72 



exp 



{X -Xq- poti/m) 

2(72 



X f mX 



h \{t2 -h] 



(7.75) 



In Chapter 5 it was noted that this type of integral is dominated by values of X for which 



0<X < 



■4|po|' 



(7.76) 



providied that Po(^2 — ti)/2m » 1, and in this range the exponential term is approximately 
constant. We can approximate this integral then by intergrating from to — 7r/4|po|, taking 
X — in the exponential term and approximating. 



This gives 



1 7/ 

f{u)^--- + 0{u% 

Z TT 



7T h f {Xo+poh/mf 

exp 7^ « 1, 



(7.77) 



(7.78) 



2 8a\po\ \ 2(t2 
there will therefore be decoherence for gaussian wavepackets tightly peaked in momentum. 
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provided their momentum is such that Po(^2 — ti)/2m >> h, ie £^0(^2 >> ^- In Chapter 
5 it was argued that this conclusion also holds for orthogonal superpositions of gaussians. 

Turning to the intermediate case, since (^2 — ti)/Ti << 1 we can set D = in Eq.(7.71) 
and we obtain Eq.(7.72) again. Now however, Wtj^{X,P) is the initial state evolved with 
an environment, and since ti/ri » 1 this should be significant. To proceed we write the 
Wigner function at time ti in terms of the initial state at t = using Eq.(7.22) 



dXodpoWo{Xo,po)F[Xo,po] 



/O /-oo 
dX / dpK{X,p,ti\Xo,po,0)f 
00 J —00 



X ( mX 



(7.79) 



(7.80) 



n \{t2-ti) 

Because the Wigner function propagator is a gaussian the analysis is similar to the first case. 
Recall that we are assuming t « Tg and so p^/Dti >> 1, this means we can integrate out 
p, setting p — pq. This gives us 



dXN, I ^ exp (-^{X-Xo- poh/mf) ) / 



F[Xo,Po] = j 
where, N,a,p are defined in Eq.(7.24). Now, 



X r mX 



+ Po 



pl/Dh » 1 



(7.81) 



(7.82) 



implies that 



1 

Po 



,2 « .2, 



n 



n'/3 \tij ' ^ ^ 

so that for ti > Ti, 1/pq << l/(5Ti^. This means, again like the case of D = 0, that the 
exponential term is roughly constant compared to /. Since we have Po(^2 — ti)/2m >> 
we follow the same procedure as the D = Q case, arriving finally at 



F[Xo,Po] « exp (-'^{X,+pot,/mf 

o\Po\ V 4 \ 4 



(7.84) 



The value of A now depends on the relationship between the width of the initial state a and 
13. However we can obtain an upper bound by ignoring the effects of the exponential term 
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in F, this gives 




A<||;A/7-^,/5f (r)«l. P.85) 



Finally we have the case of strong system-environment coupling. Since (^2 — >> 1 

the integral over ^ in Eq.(7.71) will be peaked around ^ = 0, and we can therefore extend 
the limits of integration to (— oo, oo). This integral and the one over X2 may then be carried 
out and we obtain 



Wt,{X,p), (7.86) 



where Erfc is the complementary error function [2] . For large positive values of the argument 
we have that, 

ErfcM =.£-^(1 + 0(1)). (7.87) 

Therefore, since the Wigner function is peaked around p < 0, A will be very small, provided 
(^2 — ^i) << Ts- This gives us an upper bound on the time interval, [^1,^2], rather than a 
lower one. The lower time scale is provided by the condition t2 — ti»Ti — ^^2mh/D. This 
lower time scale is compatible with the condition that the current be positive. 

Note however that this lower limit is state independent. There will be states for which 
arrival time probabilities decohere on much shorter time scales than this, for example the 
simple cases which decohere in the absence of any environment will continue to do so in the 
presence of an environment, at least until a time ~ r^. 

The key point is that whether or not one can assign arrival time probabilities in decoherent 
histories depends on the form of the state at the time it crosses the origin. Environmen- 
tally induced decoherence produces mixtures of gaussian states from generic initial ones and 
thus, after a short time, arrival time probabilities can be defined whatever the initial state. 
Crucially however it is not necessary for the system to be monitored whilst it is crossing the 
origin. The smallest time interval, 5t, over which we can define a decoherent arrival time 
probability is therefore set by the energy of the system and not the details of the environment 
and we must have E6t » U. This is in agreement with Chapter 5 and Ref.[49], and also 
with the results of earlier works, concerning the accuracy with which a quantum system may 
be used as a clock [88] . 

In conclusion then, for a general state, decoherence of histories requires that we evolve for 
a time much greater than ti — ■\/2mh/D before the first crossing time. This is because this 
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is the time scale on which quantum correlations disappear and our initial state begins to 
resemble a mixture of gaussian states. After this time, we may define arrival time probabil- 
ities to an accuracy St, provided only that E6t » h. States which start as gaussians may 
be assigned arrival time probabilities without this initial period of evolution. This is in line 
with the general result that some coarse-graining is always required to achieve a decoherent 
set of histories in quantum theory. 

7.7 Summary and Discussion 

In this chapter we have been concerned with deriving an arrival time distribution for open 
quantum systems and comparing this with the classical result. We began by discussing the 
generalisation of the current, which is the classical arrival time distribution, to open quantum 
systems and in particular to the case of quantum Brownian motion. We found that in general 
the inclusion of an environment leads to extra terms appearing in the expressions for the 
current, compared with the those valid in the unitary case. However we have shown that in 
the limit of negligible dissipation these correction terms may be dropped. We then explored 
the resulting arrival time distribution and showed that it is non-negative after a time of 
order r/ and that after this time it can be written as the trace of the density matrix times a 
POVM. 

We then turned to the question of deriving this arrival time distribution from the decoher- 
ent histories approach to quantum theory. We extended the decoherent histories analysis of 
the arrival time problem to the case of a particle coupled to an environment. As expected, the 
inclusion of an environment produces decoherence of arrival time probabilities for a generic 
initial state. There are, however, some limitations to the permitted class of histories. For a 
generic state arrival times can only be specified after an initial time t » Ti. Even after this 
time arrival times cannot be specified with arbitrary precision, coarse graining over intervals 
5t » h/E is required to ensure decoherence. We showed that the decoherence condition is 
very closely related to a semiclassical approximation for the evolution of the state, and that 
both conditions are satisfied if the time between projections is sufficiently large. This is a 
specific case of a more general connection between decoherence and classical behavior. 

Our approach has proceeded at two levels. At the heuristic level the simple generalisation 
of Eq.(7.2) to the case of a particle coupled to an environment, Eq.(7.39), is a positive arrival 
time distribution after a time of order t;. This can also be written as the expectation value 
of a POVM, Eq.(7.44), and thus has the same form as other probabilities in quantum theory. 
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On a more fundamental level, these expressions can be derived from the dccoherent histories 
approach to quantum theory, where they are seen to be valid for times much later than r^. 

Although the arrival time probabilities computed from dccoherent histories agree with the 
heuristic ones when we have decoherence, their derivation in this way represents a significant 
advance in our understanding. The great difficulty with regarding the current as the arrival 
time distribution is the arbitrary way in which one accepts these "probabilities" when they 
are positive, but declines to do so when they are not. Because decoherence is an essential part 
of the histories formalism, this arbitrariness is replaced with a consistent set of rules governing 
when probabilities may or may not be assigned. Whilst this may be of no consequence in 
the setting of a laboratory, it may prove hugely important in the analysis of closed systems, 
in particular the study of quantum cosmology [47] . 

Another very interesting result we have presented is that the current becomes strictly 
positive after a finite time whilst assignment of probabilities in decoherent histories is only 
possible asymptotically. In some ways this difference between the heuristic analysis and 
decoherent histories is to be expected. Indeed, the current represents a linearly positive 
history [35], and it is essentially the condition of linear positivity that we have proven 
holds after a time of order r/. It is known that linear positivity is a weaker condition than 
decoherence [60, 46], and thus it is not surprising that demanding decoherence leads to a 
stricter limit on the assignment of probabilities than the heuristic analysis. It would be 
interesting to examine what is gained in this context by imposing decoherence rather than 
linear positivity. 

There are similarities here with the relationship between the current and Kijowski's arrival 
time distribution [68], which may be shown to agree in the classical limit, but not more gen- 
erally. Indeed because it is manifestly positive, one might expect the arrival time distribution 
computed from decoherent histories to be more closely related to Kijowski's distribution on 
shorter time scales. These issues will be discussed elsewhere. 
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8 Quantum Arrival and Dwell Times 
via an Ideal Clock 



But what minutes! Count them by 
sensation, and not by calendars, and each 
moment is a day. 

Benjamin Disraeh 

8.1 Introduction 
8.1.1 Opening Remarks 

So far in this thesis we have been concerned with defining time observables without reference 
to the specific mechanisms by which they might be measured. This is of course reasonable, 
we do not need to specify the means by which the position of a particle is to be measured, for 
example, to be able to state that the probability distribution for finding it at a given point x 
at time t is However since obtaining a satisfactory definition of these observables 

has proven rather tricky, it is worth asking to what extent the ideal quantities we have 
been discussing so far are related to the measured quantities obtained in some particular 
measurement set up. As well as this, we mentioned in Chapter 1 that one way of defining 
time observables was via a particular measurement model. The basic theory is outlined in 
Section 1.5, we look for an arrival or dwell time distribution of the form Eq.(1.27), and then 
try to extract the ideal distribution which is independent of the measurement model. In this 
chapter we use a clock model to define arrival and dwell times and compare the results with 
standard semiclassical expressions. 
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8.1.2 Dwell Time 



Although we have talked at length in this thesis about the arrival time problem, we have 
not so far mentioned the dwell time problem in any detail. Since we will use our model clock 
in this chapter to look at both arrival and dwell times, it is useful to summerise some basic 
facts about dwell time here. For more details see Ref.[83]. 

The dwell time distribution is the probability Il{t)dt that the particle spends a time 
between [t, t + dt] in the interval [— L, L\. One approach to defining this is to use the dwell 
time operator, 

/oo 
dtxixt), (8.1) 
■oo 

where x(a;) is the characteristic function of the region [— L, L] [83]. Since To is self-adjoint 
the distribution Il{t) can then be written as 

n(t) = (V^o|5(t-T^)|V^o). (8.2) 

In the limit \p\L ^ 1, where p is the momentum of the incoming state, the dwell time 
operator has the approximate form To pa 2mL/\p\ so that the expected semiclassical form 
for the dwell time distribution is 

U{t) = {iJo\s(t-^] l^o). (8.3) 



\p\ J 

This is the form against which we shall compare our clock model. 



8.1.3 Clock Model 

In this chapter we will derive arrival and dwell time distributions by coupling the particle to 
a model clock. We denote the particle variables by {x,p) and those of the clock by {y,Py). 
We denote the initial states of the particle and clock by l^'), |0), respectively, and the total 
system state by |^^). We couple this clock to the particle via the interaction Hj — \x{x)Hc- 
The total Hamiltonian of the system plus clock is therefore given by, 

H^Hq + Xx{x)Hc (8.4) 

where Hq is the Hamiltonian of the particle. Here x is the characteristic function of the 
region where we want our clock to run, so that x{^) = ^(^) the arrival time problem and 
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x{x) — 9{x + L)9{L — x) for the dwell time problem. The operator, 

H, = H,{y,py) (8.5) 

describes the details of the dynamics of the clock and we assume that it is such that the 
clock position y is the measured time. The physical situation is depicted in Figure (8.1). For 

t 




Figure 8.1: The arrival time problem defined using a model clock. The clock runs while the 
particle is in x > 0. 

the moment we will assume only that the clock Hamiltonian is self adjoint, so that it may 
be written in the following form, 

J dee\e) (e| (8.6) 

where the |e) form an orthonormal basis for the Hilbert space of the clock. Later on we will 
restrict H^. further by considering the accuracy of the clock. We will also quote some results 
for the special choice of He = Py. whose action is to simply shift the pointer position y of 
the clock in proportion to time. This is the simplest and most frequently used choice for 
the clock Hamiltonian. The physical relevance of this and other clock models is discussed in 
Refs. [71, 88]. 



148 



Our aim, for both arrival and dwell times, is to first solve for the evolution of the combined 
system of particle and clock. We write this as 



{x,y\e-'''^\^o) 

{x, y\ exp {-iHqt - iXx{x)H^T) l^/^o) |0o) 
/ de {y\e) (e|0o) (^^l exp{-iHoT - iXx{x)eT) \ipo) 



(8.7) 



and we then solve for the propagator in the integrand using path integral methods. We will 
take the total time r to be sufficiently large that the wave packet has left the region defined 
by x(^)- We will then compute the final distribution of the pointer variable y, which is 



Our main aim is to show that the predictions of the clock model Eq.(8.8) reduce, in certain 
hmits, to the standard forms described above. Our main tools in solving Eq.(8.7) will be 
those we have used throughout this thesis, namely the path decomposition expansion and 
the semiclassical approximation it suggests. These results are covered in Chapter 2. 

8.1.4 Connections to Earlier Work 

Clock models of the type Eq.(8.4) for arrival and dwell times have been studied by numerous 
previous authors, including Peres [88] , Aharanov et al [4] , Hartle [58] and Mayato et al [71] . 
These studies are largely focused on the characteristics of clocks. Refs.[88, 4] are the works 
perhaps most closely related to the present work. They concentrate on the case of a clock 
Hamiltonian linear in momentum, with some elaborations on this basic model in the case 
of Ref. [4]. Here, we focus on a different issue not addressed by these works, namely the 
dependence of the distribution Eq.(8.8) on the initial state of the particle, for reasonably 
general clock Hamiltonians. In particular, we determine the extent to which the standard 
semiclassical forms derived above are obtained for general initial states of the particle. We 
also use path integral methods to perform the calculations, in contrast to the scattering 
methods used in most of the previous works. 




(8.8) 
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8.1.5 This Chapter 



The rest of this chapter is arranged as follows. In Section 8.2 we compute Eqs.(8.7), (8.8) 
for the arrival time problem, and similarly in Section 8.3 for the dwell time problem. We 
conclude in Section 8.4. This chapter is based on Ref.[101]. 

8.2 Arrival Time Distribution from an Idealised Clock 

We now turn to the calculation of the arrival time distribution, Eq.(8.8), recorded by our 
model clock. Using the path decomposition expansion in the form Eq.(2.13) the state of the 
system Eq.(8.7) can be written as. 



^0 

We can simplify this expression in two different regimes, the weak coupling regime of ^ Ae, 
and the strong couphng regime of <^ Ae. 

8.2.1 Weak Coupling Regime 

In the limit E ^ Xewe can make use of the semiclassical approximation to the PDX formula. 
This yields. 



*(x,y,T) 





x/ dt{x\exp{-iHo{T -t j)6{x)pcxp{-z{Ho + Xe)t)\xpo) . (8.10) 



This means that the arrival time distribution is 




X / dtdt' {ipo\ ex.p{i{Ho + Xe')t')p5{x) ex.p{-iHo{t' - t))5{x) 



xpexp(-i(ifo + Xe)t) \ijjo) 



(8.11) 
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To proceed, we first note that for any operator A, we liave 

S{x)A6{x) = 6{x){0\A\0) . (8.12) 

Using tliis in Eq.(8.11) gives, 

n{y) = —J dede' (0o|e') (e'ly) {y\e) (e|0o) 



X / dtdt' (V'ol exp(i(ifo + Xe')t')pS{x)pexp{-i{Ho + Xe)t) \ijjo) 
Jo 

X {0\exp{-iHo{t' -t))\Q) . (8.13) 

We see here the appearance of the combination p5{x)p, and the main challenge is to show 
how this turns into the current operator, 5{x)p + p5{x). 
Next we rewrite the integrals using 

pT pT pT pT pT 

/ dtdt' = / dt j dt'+ dt' / dt. (8.14) 
Jo Jo Jt Jo Jf 

In the first term we set u — t, v — t' — t, and in the second we set u — t', v — t — t' to obtain, 

n(y) ^ I dede'(0o|e')(e»(?/|e)(e|0o) f du f ^ dv 
^ J Jo Jo 

X {{tpol exp{i{Ho + \e)u)p6{x)p ex-p{-i{Ho + Xe){u + v)) [tpo) (0| cxp{iHov) |0) 
+ (V^ol exp(i(i/o + Xe'){u + v))pd{x)pexp{-i{Ho + Ae)^) iV'o) (0| exp{-iHov) |0)} 

(8.15) 

Since wc take the time r to be large, wc can extend the upper limits of the integrals to 
infinity. The integral over v can then be carried out, to give. 



U{y) = I dede' {(Pole') {e'\y) {y\e) {h\<Po) 

2^ du{i;u\e'^'^{pS{x) + S{x)p)e-'''^^\i;u) 



n(y) = du\^y,u)\'{^u\{p5(x) + 5(x)p)\^|;u) 

dtmy,t)\^J{t) (8.16) 
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where 




(8.17) 



is the wavefunction of the clock, and J{t) is the current, Eq.(1.22). 

This form shows that, in the weak couphng hmit, our arrival time probability distribution 
yields the current, but smeared with a function depending on the clock state. We thus get 
agreement with the expected result, Eq.(1.27). Note that the physical quantity measured, 
the current, is not affected by the form of the clock Hamiltonian. 

Although the form Eq.(8.16) holds for a wide class of clock Hamiltonians, not all choices 
make for equally good clocks. To further restrict the coupling He we require that different 
arrival times may be distinguished up to some accuracy St. For this to be the case we re- 
quire that the clock wavefunctions corresponding to different arrival times are approximately 
orthogonal, so that 



where St = t — t'. Clearly this expression is equal to 1 ii St = 0. Suppose now that |0o(£)P 
is peaked around some value eo with width Ug. This integral will approximately vanish if 



and so the resolution of the clock is given by l/Xa^. The relationship between t and the 
pointer variable y will depend on the specific model. It is easily seen that a clock with 
good characteristics may be obtained using, for example, a free particle with a Gaussian 
initial state. But clocks with more general Hamiltonians can also be useful if they evolve 
an initial Gaussian along an approximately classical path (as many Hamiltonians do). See 
Refs.[88, 4, 58, 71] for further discussion of clock characteristics. 

For the special case He — Py, \e) — \py), the expression for the arrival time distribution 




(8.18) 



We easily see that. 




(8.19) 



Xa^St » 1, 



(8.20) 
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simplifies, since 



(8.21) 



The time is related to y hj t = y/X and the expected form Eq.(1.27) then becomes a simple 
convolution. 

8.2.2 Strong Coupling Regime 
Special Case: He — Py 

We now turn to the limit of strong coupling between the particle and clock. The analysis 
of the case of general clock Hamiltonian is rather subtle, so before we tackle this we first 
examine the special case where the clock Hamiltonian is linear in the momentum. That is, 
we have. 



We start from Eq.(8.9), and insert a complete set of momentum states for the particle, p to 
obtain. 



where (poipy) is the initial momentum space wavefunction of the clock, and E — p^ /2m is 
the kinetic energy of the particle. Note the appearance of the momentum p in the integrand. 
The expression involving the integral over t has been computed previously using the final 
crossing PDX [41, 49]. In the hmit r — > oo it is given by. 




(8.22) 




(8.23) 




dt {x\ exp{-i{Ho - Xpye{-x))t) |0) exp{iEt) 
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We can now write our probability distribution for y. Carrying out the x integral we obtain, 



n(?/) 



dpydp'ydpdp'(t)l{p'y)(t)Q{py) {p'y\y) {y\py) ex.p{-iX{py - Py)T) 



X— W) (pI^o) exp{-t{E - E')t) 




=S U2m{E + Xpy) - ^2m{E' + Ap^) 



(8.24) 



Using the formula S{f{x)) — 5{x)/f{0), we can carry out the p'y integral to give, 



n(2/) 




dpydpdp\(t)Q{py)\^— (V^oIpO {p\^o) exp I -i 



exp i-ffov p5(£)pexp -iHq- iV'o) 



(8.26) 



where A is some constant whose explicit form is not required, and we have used the fact 
that E « Xpy. We see therefore that in this limit the probabihty of finding the clock at 
a position y is given by the kinetic energy density of the system at the time t — y/X, in 
agreement with Eqs.(1.25) and (1.26). 

Note that there is no response function involved in this case, as one might have expected 
from the general form Eq.(1.27). (A similar feature was noted in the complex potential 
model of Ref.[41]). It seems hkely that this is because the strong measurement prevents the 
particle from leaving x > until the last moment, so that the response function R{t, s) is 
effectively a delta-function concentrated around the latest time. 

General Case 

As well as the approximations valid for E « Xpy, the key to the analysis in the special 
case presented above is that the position space eigenfunction of the clock Hamiltonian with 
eigenvalue Py takes the simple form. 



{y\py) 



1 



exp{iypy). 



(8.27) 



V2n 
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This greatly simplifies the resulting calculation. For the case of a more general clock Hamil- 
tonian, the eigenstates will not have this simple form. Instead we make a standard WKB 
approximation for the eigenstates of the clock, 

(y|e) = C(|/,e)exp(z5(y,e)) (8.28) 

where S{y,e) is the Hamilton Jacobi function of the clock at fixed energy. This means 
Eq.(8.25) becomes, 

U{y) ^ J dedpdp'l (6|0o) {Mp') (pI^o) ^6 + ^^^ly^ {y\e) 
- Jde\{e\<^o)\'^,/^\Ciy,e)\^ 

X {^o\ exp (^zHo^^^^^^ p5{x)pe^p (^-iH,^^!^^^ \^,) (8.29) 

where we have used, 

(e+ (^1.) m « ICfe.ejPexp (-,(^^) , (8.30) 

which is vahd for £■ — £"<< Ae. 

We now suppose that the clock state is a simple Gaussian in y, or equivalently in py. 
It follows that it will be peaked in e about some value eo- This means that the integral 
over eo in Eq.(8.29) may be carried out. The result for Il{t) will again be proportional to 
the kinetic energy density, of the form Eq.(1.25), where the relationship between t and the 
pointer variable y is defined by the equation 

A de 

as one might expect from Hamilton- Jacobi theory [34]. Hence the arrival time distribution 
has the expected general form, Eq.(1.25) (and therefore Eq.(1.26) also holds), but the precise 
definition of the time variable depends on the properties of the clock. 
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(8.31) 



8.3 Dwell Time Distribution from an Idealised Clock 



We now turn to the related issue of dwell times. Here the aim is to measure the time spent 
by the particle in a given region of space which, for simplicity, we take to be the region 
[— L,L]. This is portrayed in Figure (8.2). In this section we will work exclusively in the 
weak coupling regime where E » Xe. 



t 




Figure 8.2: The dwell time problem, defined using a model clock. The clock runs while the 
particle is in the region [—L,L\. 

The starting point is the final state of the particle plus clock, Eq.(8.7), which we write as 

y, r) = {x,y\ e-*(^o+^^^x(£))r |^^^ _ (8 32) 

where xi^) = 6{x + L)6{L — x). Wc wish to rc-cxprcss this using the path decomposition in 
a similar way to Eq.(8.9). For this case wc need a PDX which is more general than the one 
used for the arrival time, since there are now crossings of two surfaces. One way to proceed 
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is to use the path integral expression for the first crossing of x — L and x — — L, which is 

'^{x,y,r) = — / de {e\(f)o) {y\e) / ds dt {x\ex.p{-i{Ho + Xx{x)e){T -t - s))\-L) 
J Jo J-oo 

{-L\pex.p(-i{Ho + Xx(x)e)s) \L) {L\pexp(-iHot) |V'o) , 

This is shown in Figure (8.3). However there are other choices we could make. We could 

t 



Xi 




Figure 8.3: The PDX used for the dwell time problem: Paths from xq > Q to x\ < Q have a 
first crossing of x — L at a time t and a first crossing of x — —L at a time t + s. 

consider the first crossing oi x — L and the last crossing oi x — — L for example. In the 
semiclassical limit these choices lead to equivalent expressions for the dwell time. It would 
be interesting to explore what differences do arise in other regimes. This will be addressed 
elsewhere. 
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It will prove more useful to work with the wavefunction in position space for the clock and 
momentum space for the particle. Changing to this representation, and making use of the 
semiclassical approximation we obtain, 

{p,y\-^r) ~ de{y\e){e\(t)Q) ds dt {p\e^v{-iHG{T - 1 - s))\-L) 

TT^ J Jo J-oo 

X {-L\pex.p{-i{Ho + Ae)s) \L) {L\pexp{-iHot) \ipo) (8.33) 

Now we make the standard scattering approximation of letting the upper limits of the inte- 
grals go to infinity. This means we can carry out the t and s integrals to obtain, 

(p, y\^r) ^ Jde {y\e) (e|0o) exp(-i£;T) exp(-i2LmAe/|p|) (p|Vo) , (8.34) 

where we have used the standard integral [49], 

poo 

iHos)p\0)e'^' - mex-p{i\x\V2mE) (8.35) 

Jo 

and 

f°° dt 

/ —{x\peM-KHo-E)t)\^o) = {x\p5{Ho-E)\^o)=m{x\5{p-p)\%), (8.36) 

Here, we have neglected the term involving 5{p + p) since this corresponds to reflection, 
which will be negligible in this semiclassical limit. We have also use the fact that E 3> Ae. 
We therefore obtain the distribution for y as, 

Ii{y) = J dp\{p,y\^r)\' 

= I dede' (0o|e') (e'|y) {y\e) (e|0o) J dp\Mp)\' (^^(^' - ^)) 

= J dp\Mp)\'My,'2Lm/\p\)\' (8.37) 

where, 

—i2LmX 



^y,2Lm/\p\) = J de {y\e) {e\(f)o) exp 



\P\ 



-e 



(y|exp( i^c)|0o) (8.38) 
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is the clock wavefunction. We may rewrite this as 

U{y) ^ Jdt my,t)\'{i^o\S (^t - ^) iV^o) (8.39) 

It is therefore of precisely the desired form, Eqs.(8.3), (1.27), with |$(y,i)p playing the role 
of the response function. The discussion of clock characteristics is then exactly the same as 
the arrival time case discussed in Section 3. 

8.4 Conclusion 

We have studied the arrival and dwell time problems defined using a model clock with a rea- 
sonably general Hamiltonian. We found that in the limit of weak particle-clock coupling, the 
time of arrival probability distribution is given by the probability current density Eq.(1.22), 
smeared with some function depending on the initial clock wave function, Eq.(1.27). This is 
expected semiclassically, agrees with previous studies and is independent the precise form of 
the clock Hamiltonian. 

In the regime of strong coupling, we found that the arrival time distribution is proportional 
to the kinetic energy density of the particle, in agreement with earlier approaches using a 
complex potential. The fact that two very different models give the same result in this 
regime suggests that the form Eq.(1.26) is the generic result in this regime, independent of 
the method of measurement. It would be of interest to develop a general argument to prove 
this. (See Ref.[23] for further discussion of this regime). 

For the case of dwell time, we have shown that the dwell time distribution measured by 
our model clock may be written in terms of the dwell time operator in semiclassical form, 
smeared with some convolution function, Eqs.(8.3), (1.27). 

In all of these cases, the precise form of the clock Hamiltonian and clock initial state 
determine the relationship between time t and the pointer variable y and they determine the 
form of the response function R in the general form Eq.(1.27). These are particularly simple 
for the special case He = Py explored previously. However, what is important is that, once 
the definition of the time variable is fixed, the clock characteristics do not effect the form of 
the underlying distributions - the n(s) in Eq.(1.27). The n(s) are always one of the general 
forms Eqs.(1.22), (1.25) and (8.3), no matter what the clock characteristics are. This means 
that these general forms will always play central role, irrespective of how they are measured. 



159 



9 Summary and Some Open 
Questions 



This is the first thing 
I have understood: 
Time is the echo of an axe 
Within a wood. 

Phihp Larkin 

9.1 Summary 

The main result presented in this thesis is an analysis of the arrival time problem in the 
framework of decoherent histories. This analysis has given us an axiomatic way of under- 
standing how probabilities may be assigned to this observable, as well as the conditions 
under which these probabilities are meaningful. We have carried out this analysis both for 
a free particle, where decoherence is achieved for a limit class of states, and for a particle 
coupled to a environment, where decoherence is achieved for a general class of states after a 
given time. The main tool which allowed us to proceed was the equivalence between pulsed 
and continuous measurements set out in Chapter 4. This result allows us to construct the 
appropriate class operators either by writing them in terms of complex potentials, or by 
giving us the confidence to use semiclassical approximations provided the energies and times 
involved satisfy the conditions laid out in Chapters 4 and 5. The major advance achieved 
compared to previous work is a decoherent histories based framework that respects the Zcno 
limit. In addition we have also shown how the semiclassical, heuristic formulae may be de- 
rived from framework. In Chapter 8 we also demonstrated that the probabilities obtained 
via decoherent histories agree with those computed from model clocks. This general scheme 
of decoherent histories, using path integral methods and complex potentials to define class 
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operators, seems to be general enough to encompass any time observable in quantum theory. 
We encourage the reader to test this for themselves. 

9.2 Open Questions and Future Work 

We list here, with some short comments, a few obvious questions which arise from the work 
presented in this thesis, but which time constraints have prevented the author from studying 
in any depth. The majority of these questions would form the basis of interesting further 
study. 

Decoherent Histories Analysis of the Dwell Time Problem 

The most obvious question one might ask of the work presented here is whether it can be 
extended to cover other time observables. The most obvious such extension would be to 
the definition of dwell times. Dwell times were discussed very briefly in Chapter 8 where 
they were defined using a model clock, but no attempt has been made here to perform a 
decoherent histories analysis. The key step in carrying this out lies in the construction of 
the relevant class operators. One might imagine doing this either via the introduction of a 
suitable complex potential or potentials in the manner of Chapter 5, or via some semiclassical 
approximation in the manner of Chapter 6. What is particularly interesting about the dwell 
time problem is that there exists a relatively straightforward operator formalism, sec Ref.[83], 
and it is not clear at the moment how the decoherent histories analysis fits in with this. 

Quantum Zeno Effect for Open Systems 

Another topic which we touched on briefly in Chapter 7 was the occurrence of the quantum 
Zeno effect in an open system. Whilst we studied the evolution of a particle undergoing 
continuous measurement in detail in Chapter 4, it seems at first sight a difficult task to extend 
the analysis to a particle coupled to some environment. As we pointed out in Chapter 7, 
this is an important question since it is the Zeno effect that limits the accuracy with which 
we can define time observables in quantum theory. On the face of it one might expect 
that the Zeno effect should somehow vanish in the classical limit, so that inclusion of an 
environment should reduce the "Zeno time" from 1/E. However it is not at all obvious 
how this might come about, since one generally assumes that all time scales associated with 
the environment are much longer than 1/E. In part, this is related to the issues that for a 
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genuine step potential the reflection and transmission coefficients, Eqs.(2.22) and (2.24), do 
not depend on h and thus have no naive classical limit (see Chapter 2). 

Time Observables and the Backflow Effect 

Another question which we have encountered but not discussed in any detail in this thesis 
concerns the role of the backflow effect in the deflnition of time observables. The backflow 
effect causes the current Eq.(1.22) to be negative even when classically it should be positive, 
and thus prevents decoherence of histories, see Section 5.5.4. However unlike the quantum 
Zeno effect, we do not have a good grasp of the time scale on which this effect occurs. Whilst 
it is true that one can always flnd a state which displays backflow for an arbitrarily long 
period of time [12], one suspects that for a given state this time interval is determined by 
some property of the state such as the moments of the Hamiltonian. Since it it the backflow 
effect and the Zeno effect that together determine the accuracy with which time observables 
may be defined, further investigation of these effects would be useful. Work on this has 
already begun [102] . 
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Appendix: Some Properties of the 
Current 



We have derived the expression 



p{0,T) = [ dt J{t) 



(9.1) 



as the approximate probabihty for crossing the origin during the time interval [0, T], where 
J{t) is the usual quantum-mechanical current. The current itself is not necessarily positive 
due to backfiow. Here we explore the possibility that averaging it over time might improve 
the situation. On the one hand, the results of Bracken and Melloy [12] show that there is 
always some state for which p(0, T) defined above is negative, for any T. On the other hand, 
for a given state, one might hope that p(0, T) will be positive for sufficiently large T. Here 
we give a brief argument for this, which also makes contact with the negativity of the Wigner 
function. 

The current can be written in terms of the Wigner function W{p, x) as 



The Wigner function evolves freely according to W{p, x, t) — W{p, x — pt/m, 0). We assume 
it has support only on negative momentum states, with average momentum po < and 
momentum width (Tp. 

Consider a time interval < t < T over which backfiow occurs. It is clear that in order 
for this to occur the Wigner function must be negative for at least some of this interval. We 




(9.2) 



163 



can write, 



p{0,T) 




\p\T/m 

dx W{p,x,0) 




(9.3) 



So p{0, T) is given by the average of the Wigner function over a region of phase space. We 
now recall a standard property of the Wigner function which is that, broadly speaking, it 
will tend to be positive when averaged over a region of phase space of size greater than 
order 1 (in the units used here where % = 1). This region is of size of order \pq\T /m in the 
x-direction but infinite in the p-direction. However, the Wigner function has momentum 
spread Cp, so the effective size averaged over is apP^T/m which is approximately the same 
as Ai7 T. This means that we expect that p(0, T) will be positive as long as 



Hence, as expected, the integrated current will be positive for T sufficiently large and the 
key timescale is the Zeno time. This heuristic argument will be revisited in more detail 
elsewhere. 
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